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I. THE NEUTRAL ATMOSPHERE 


1 . Atmospheric Nomenclature 

The scientific study of the upper atmosphere has become known as 
aeronomy. As generally happens in the study of a complex natural system, 
nomenclature has been developed to describe the different parts of the 
atmosphere (Chapman, I960). The description may be based on chemical 
composition or on temperature or on the dominant physical processes. For 
each of these alternatives there exists a series of words terminating in 
"-sphere, " each using some property to characterize the atmosphere within 
a certain range of altitude. The upper boundary of any "-sphere" may be 
denoted by a similar word ending in "-pause"; thus, the "tropopause" is 
the upper bound of the troposphere. Sometimes these levels can be defined 
to within a few kilometers, but in other cases, the "-pause" nomenclature 
is inappropriate because the divisions can only be specified within tens 
or hundreds of kilometers. Figure 1 contains all the "-sphere" terms that 
occur in this review. 



FIG. 1. REGIONS OF THE ATMOSPHERE, SHOWING 
CONVENTIONAL NAMES DESCRIPTIVE OF LEVELS, 
PHYSICAL REGIMES, AND CHARACTERISTIC CON- 
STITUENTS, The temperature profile is 
taken from the U.S, Standard Atmosphere 
and the electron density profile represents 
average daytime conditions for middle 
latitudes, high solar activity. 
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The lowest layer is the troposphere , in which the composition is 
uniform and the temperature generally decreases upward. The emission and 
absorption of infrared radiation by molecules such as water vapor, carbon 
dioxide and ozone provides the most efficient transfer of heat between 
different levels in this region. If we solve simplified equations for 
radiative transfer, we find that the equilibrium temperature in the lower 
atmosphere decreases approximately linearly with the partial pressure of 
these constituents. At low heights, this would lead to a "lapse rate" 
(negative temperature gradient) exceeding the "dry adiabatic lapse rate," 
which is about 10 °K/km. Such a situation is unstable and the lapse rate 
is maintained by turbulence and water vapor condensation at a somewhat 
smaller value, about 7 °K/km, throughout the troposphere. Under certain 
conditions, particularly at night, "inversions" may be set up in which 
the temperature gradient near the ground is positive. 

Until the turn of the century, it was supposed that the temperature 
continues to decrease upward and that the atmosphere terminates at about 
50 km, there merging into cold interplanetary space. However, experiments 
with kiteborne thermometers revealed a nearly isothermal region beginning 
at an altitude of about 11 km in mid-latitudes at a temperature of about 
220 °K which has since been named the stratosphere . Its base, named the 
tropopause, is somewhat higher, and the stratospheric temperature is about 
20 °K cooler at the equator than in the polar regions. Uniformity of 
temperature is maintained principally by radiation exchange in the infrared 
part of the spectrum between the minor constituents of the atmosphere. 

The existence of a temperature inversion or positive temperature 
gradient above the stratosphere was suggested by observations of sound 
propagation over distances of 100 km or more, which seemed to result from 
the refraction of sound waves in the upper atmosphere. We now know that 
such an inversion results from the presence of a trace of ozone in the 
atmosphere. Although it i,s a very minor constituent, only contributing 
a few millionths of the total ground level pressure, the presence of ozone 
is very important. It absorbs all solar ultraviolet radiation of wave- 
length less than 2900 X, and partially absorbs wavelengths between 2900 X 
and 3600 X. The atmosphere is heated by this absorption and the earth's 
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surface is shielded from the otherwise lethal radiation. This region of 
elevated temperature is known as the mesosphere . It does not seem to have 
a well-defined lower boundary ( stratopause) but its upper boundary (meso- 
pause) lies at 80 to 85 km, and is the coldest level in the entire atmo- 
sphere (about 180 °K). Heat flows toward this level by conduction, from 
above and below, and can only be removed by radiation which is primarily 
in the infrared, such as the 63p line of atomic oxygen. The existence of 
the temperature minimum arises from the lack of any strong heating mech- 
anism at this height. Ozone cannot exist in appreciable quantities at 
this level or higher, being rapidly destroyed by photochemical reactions. 

The shorter UV radiations are absorbed at greater heights, in the 
thermosphere, and are responsible for the high temperatures existing there. 
Most of the heat liberated in the thermosphere is removed by downward con- 
duction, so that the temperature increases upward. Finally, the heat con- 
ductivity becomes so good that the region of the upper thermosphere is 
maintained nearly isothermal at a relatively high temperature (about 
1000 °K ) . In the exosphere , collisions between molecules are so infre- 
quent that neutral particles move in ballistic orbits subject only to 
gravity, whereas the ionized particles are constrained by the magnetic 
field. No more will be said about this region now, since a later section 
(l-6) is devoted to it. 

The foregoing description is closely linked to the vertical tempera- 
ture profile, which is shown in Fig. 2. Parallel to this, a classification 
in terms of composition can be made. The term, ozonosphere , may be prin- 
cipally associated with the mesosphere, but is not too well defined 
because appreciable concentrations of ozone exist in the stratosphere. 

The vertical distribution of ozone extends, roughly, between 5 and 60 km. 
The peak of concentration lies at about 30 km, though on occasions two 
peaks have been detected, the lower at 10 - 15 km. The T 'median of the 
ozone distribution is at about 25 km. The maximum rate of heat production 
occurs higher up, at about 45 km, near the level of unit optical depth for 
the incident solar radiation. 

The ionosphere may be defined as the part of the earth's upper atmo- 
sphere where ions and electrons are present in quantities sufficient to 
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affect the propagation of radio waves. It extends down to perhaps 50 km 
and thus overlaps the ozonosphere, It has no well-defined upper boundary, 
but merges into (or may be extended to include) the heliosphere , where 
neutral and ionized helium are important constituents, and the protono- 
sphere , which is composed principally of ionized hydrogen. The latter 
"-spheres" are still poorly defined in extent. 



30952 

FIG. 2. U.S. STANDARD ATMOSPHERE 
(Government Printing Office, 

Washington, D.C., 1962). Verti- 
cal distribution of pressure p, 
density p, temperature T and 
mean molecular mass M to 250 
km. The composition is assumed 
constant up to 100 km. 

The remaining terms used in Fig. 1 refer to the physical regimes at 
different levels in the atmosphere. Up to about 100 km, the atmosphere 
is well mixed by turbulence. The relative abundances of major constitu- 
ents may be assumed independent of height, except at the highest heights. 
The abundances of chemically active constituents, such as ozone, are of 
course subject to variation, as are constituents such as water vapor and 
contaminants in the troposphere. In the lower thermosphere, the composi- 
tion can be modified by photochemical reactions, such as dissociation of 
molecular gases. 
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The level at which turbulence ceases may be called the turbopause ; it 
is rather sharply defined, and lies at about 100 km. Probably this depends 
on the increase of. viscosity with height. At greater heights, the lack of 
turbulence enables a condition of diffusive separation to be established, 
in which the vertical distribution of each neutral gas depends on its 
molecular weight. The distribution of chemically active gases may however 
be influenced to some extent by photochemical reactions as well as by dif 
fusion. This is especially true of the ionization, which attains a dif- 
fusely controlled distribution only at heights above 400 km. 

Finally, there is the term magnetosphere , denoting the region in which 
the earth's magnetic field controls the dynamics of the atmosphere. It is 
difficult to define a lower limit since the movement of ionization is geo- 
magnetically controlled at all heights above about 150 km (or even less), 
but the magnetosphere certainly includes the whole atmosphere above the 
level at which ionized constituents predominate over neutral, probably at 
about 1500 km. Because of this magnetic control, the earth's atmosphere 
may be said to terminate at the magnetopause , the boundary of the geomag- 
netic field which lies at 10 to 20 earth radii (see Sec. VII ). 

2. Evolution of the Earth's Atmosphere 

It would be very surprising, indeed, if the atmosphere as we know it 
today had remained relatively unchanged in the several thousand million 
years in which our planet has been in existence. There has naturally been 
much study directed toward understanding the sequence of events which have 
combined to produce our nitrogen-oxygen environment and a variety of 
theories have been developed to explain it. The origin of the atmosphere 
must be closely linked to the development of the solid earth, but geologists 
are not in agreement concerning the earliest history of the planet. Most 
theories may be considered within a "molten globe hypothesis" or an 
"accretion hypothesis" (Chamberlin, 1949). In the former, the globe was 
originally molten and developed a crust as it cooled, while in the latter, 
the planet was accumulated by the gradual collection of small masses of 
solid material. 
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In either event, the relative abundance of the various elements in 
the solid earth and in the atmosphere may be compared with the ’’cosmic 
abundance” of the elements as estimated from meteorites and solar studies. 
When this comparison is made (cf. Kuiper, 1949), the lighter elements and 
particularly the rare gases of the earth are found to be less abundant by 
five to ten orders of magnitude. Oxygen is the first exception which 
appears with near normal cosmic abundance. This would suggest that any 
gases present early in our planet ! s history must have escaped and that in 
some way the oxygen was retained, perhaps in the rocks. 

The present atmosphere may then be one primarily established by vol- 
canic emissions. Water vapor and carbon dioxide are the main gases lib- 
erated while nitrogen concentrations are also important. Photosynthetic 
action by plants would seem to provide an adequate explanation for the 
conversion of CO^ to oxygen (Chamberlin, 1949) accompanied by the deposi- 
tion of carbon in plant material. 

Berkner and Marshall ( 1963) argue very forcefully that the primordial 
atmosphere must have been one in which molecular oxygen was very scarce, 
at no more than 0.1 percent of its present atmospheric level. The oxygen 
produced by the photodissociation of water vapor would lead to production 
of 0 3 (thereby shielding the water vapor from further dissociation ) and 
would also react with surface materials to form oxides. Concentrations 
of 0 2 substantially greater than 0.1 percent would appear to require the 
generation of oxygen from biological sources. 

The work of Hoering and Abelson (l96l) places the beginnings of photo- 
synthesis at more than 2.7 x 10 years ago and from this point the equili- 
brium levels of oxygen abundance may have begun their slow increase. As 
a convincing illustration of the importance of photosynthesis and respira- 
tion to our present oxygen levels ( Rabinowitch, 1951 ), it is found that 
all the carbon dioxide in the atmosphere is exchanged within about 400 
years, all the oxygen in about 2000 years and all the waters in the oceans 
may have participated in the biological process within about two million 
years. 

However, it was only in the much more recent Cambrian period (600 
million years ago) that higher forms of life, leaving fossil remains, 


SEL-64-111 


6 



came into existence. With oxygen levels at perhaps 1 percent of their 
present amount, it was possible for organisms to exist near the surface 
of the oceans without being destroyed by the intense ultraviolet radiation 
from the sun. Over the course of the next two hundred million years, 

Berkner and Marshall suggest that the oxygen content rose to about 10 per- 
cent of present levels, at which" time _ the ozone shielding was adequate to 
permit oxygen-breathing organisms to live on land. They identify these 
two "critical levels" at which life could first exist near the water’s 
surface and then on the land, with the interrelated rise of oxygen abundance 
and ultraviolet absorption by ozone. It is interesting to observe that 
current models of the Martian atmosphere, in which nitrogen and carbon 
dioxide are considered to be the principal constituents, may be character- 
istic of the earth's primitive atmosphere. 


3 . Structure of the Atmosphere 

The distribution with height h of a neutral atmospheric gas may be 
assumed to be subject to the perfect gas law 

p = nkT ( I -1 ) 


and the hydrostatic or barometric equation 


dp 

-= nmg 


( 1 - 2 ) 


The symbols p, p, and n are used to denote pressure, density, and con- 
centration. If it is necessary to specify any particular gas, the chemical 
symbol is appended in square brackets, e.g., nfN^l* T denotes absolute 
temperature and g the acceleration due to gravity. Molecular mass is 

denoted by M (in atomic units) and m (in conventional units) so that 

23 

p = nm, The ratio M/m is denoted by N, which is 6 X 10 per gram-mole 

20 

(Avogadro's number), but takes the value 6 x 10 per kilogram-mole in 

mks units. If R stands for the molar gas constant and k for Boltzmann's 

constant, then 
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« = “ - VT 

k m aN 


(1-3) 


By combining these equations, we find that 

(1-4 

These equations provide two alternative definitions of the scale height 
H (more precisely, the pressure scale height). For each individual gas, 
these definitions are consistent at a level where diffusive separation 
exists and where the distribution is determined solely by the balance of 
the partial pressure gradient against gravity. They are also valid for 
a gas in a completely mixed atmosphere, in which processes such as tur- 
bulence are strong enough to cause each gas to conform to the scale 
height of the mixture as a whole. The masses M and m then refer not 
to the individual gas, but to the "composite" gas, air. 

It is sometimes convenient to measure altitude in units of the scale 
height H, from some reference height h . The reduced height z, mea- 
sured in this way, is defined by the equations 


1 dp Mg = 1 
p dh " RT " H 



o 


The second equality applies if H is independent of height, but in 
any case, the vertical pressure distribution is given by p/p^ - e Z . 

This relation is valid for the partial pressure of any constituent which 
obeys the hydrostatic equation, as well as for the whole atmosphere, pro- 
vided that the proper scale height is used in the equation for z. 

Also, we see that if the whole atmosphere above any arbitrary level 
z = 0 were compressed to a uniform pressure p , its vertical extent 
would be the scale height at z = 0, wherever this level is chosen. 

The total number of particles above this level is n H . as can be deduced 

o o' 

from the gas law and the definitions of H and z. Since n = p/kT = p/mgH, 
we have 
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nH dz 


n H 
o o 


( 1 - 6 ) 



n dh 





-z 


e 


dz = 


We emphasize that these equations assume m and g to be independent 
of height, and thus apply to a single gas in hydrostatic equilibrium. In 
this case, gradients of H arise only from gradients of T, so that 
dH/H = dT/T. By taking logarithms and differentiating the gas law, we 
obtain the useful relations 


__ _ dh _ dp _ dii dH 
H ~ p ~ n + H 

Although many authors do not use the concept of reduced height, we 
shall find it useful in subsequent discussions. For example, let us 
examine the changes which occur when a column of air is heated. Each 
small volume of thickness 6h will be called a "cell," If there is no 
horizontal flow of air, the weight of air above any cell does not change; 
therefore, the pressure within the cell remains constant and so does its 
reduced height z. If we assume that the column of air is initially at 
a uniform temperature T, and is then heated to a uniform temperature 
tT, the gas concentration with each cell is reduced from n to n /t and 
its thickness expands to T&h. Since the scale height has increased from 
H to TH, the real altitude, h, measured from the ground ( z = 0 ) is 
increased in the same ratio: that is, from zH to zTH, If n is the 

initial concentration at z = 0, then the concentration at a distance h 
above z = 0 is r^e Z , or n^ exp (-h/H). After the heating, the 
vertical distribution is given by 


(1-7) 


n(h) 




( 1 - 8 ) 


H being the original scale height. For small values of h, this equation 
shows that the concentration n is decreased by the heating. But at 
greater heights n is increased because the reduction of the exponent in 
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the equation more than compensates for the factor i/t. At some inter- 
mediate level, therefore, n is unchanged, and for small changes of 
temperature (t » l) this may be shown to be just one scale height above 
z = 0. This is known in meteorology as the isopycnic level. Obviously, 
the total integrated concentration throughout the atmosphere, / n dh, is 
unaffected by the heating, in the absence of horizontal flow. 

Although the above discussion relates to a case of an isothermal air 
column in contact with the ground, it can easily be generalized to apply 
to the diurnal temperature variations in the thermosphere. There are 
three important differences, which complicate the analysis but do not 
affect the principle of the discussion. 

First, the part of the atmosphere subject to the heating does not 
extend to the ground, but is bounded by the temperature minimum (mesopause) 
at about 80 km. Secondly, the temperature in the thermosphere is not 
uniform but varies with height. Thirdly, the atmosphere at this height 
is subject to diffusive separation. Consequently, the hydrostatic equa- 
tion must be solved for each constituent separately with the appropriate 
scale height being inserted for each. When the composition at the lower 
boundary of the heated region is kept constant, it is found that at any 
fixed pressure-level the relative concentrations of all the gases are 
unaffected by the heating (Garriott and Rishbeth, 1963). 

In detailed calculations of atmospheric parameters, it must be remem- 
bered that the acceleration due to gravity decreases upward, in accordance 
with the inverse square law 


g(h) 


gm e 

2 


r 


g(o) 


r 


2 

E 


♦ -> 2 


(1-9) 


where G is the gravitational constant, M and r are the earth’s mass 

E E 

and radius, r = r + h is radial geocentric distance and g(o) is the 
sea-level value of g. This variation may be included in the equation for 
the scale height, Eq. (l-4). Sometimes, a "geopotential height" is used, 
which may be defined in various ways. The simplest definition is that in 
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which the work done in raising unit mass to height h, e.g., / g(h) dh, 
is equated to the geopotential height h* multiplied by the sea-level 
value, g(o). In addition, centrifugal force may become appreciable in 
problems involving the exosphere. 

Various "standard” model atmospheres have been proposed. They present 
detailed distributions of pressure, density and temperature which are 
based upon observations and computed with the aid of the equations given 
above. At heights where the mean molecular weight is independent of 
height, knowledge of the vertical distribution of one of the basic param- 
eters, plus boundary assumptions, enables the others to be computed from 
Eqs. (i-l) and (l-2). This is not true in the thermosphere above 200 km 
where the molecular weight varies with height. Here the best determined 
parameter is density, as determined from satellite drag data. But even 
if the complete density distribution p(h) is available, further informa- 
tion is required in order to compute the temperature profile and gas con- 
centrations. The calculations which have been published generally start 
from an assumed chemical composition at a lower boundary level (variously 
taken as 100 to 150 km). If diffusive separation is assumed to exist 
above this level, the variation of mean molecular weight with reduced 
height, z, can be found at once by applying the equation p / p q = e~ Z to 
the partial pressure of each individual gas. But conversion of reduced 
height z to real height h involves knowledge of the temperature pro- 
file T(h). In some published models, this is found by a trial -and -error 
fitting of the model to the observed density profile p(h) (e.g,, Yonezawa, 
I960) or by detailed solution of the heat -conduction equations (e.g., 
Nicolet, 1959). In the very detailed models of Harris and Priester 
( 1962a, b) , the time-varying conduction equations are solved for the tem- 
perature distribution T(h,t), for different heat inputs corresponding 
to several different levels of solar activity. Moreover, none of the 
models so far published takes account of possible variations of composi- 
tion at the lower boundary, nor of horizontal transport of air, which 
might have an important influence on thermospheric structure. 

Figure 2 gives the vertical distributi on of temperature, pressure, 
and density from 0 to 250 km, as tabulated in the U.S. Standard Atmosphere 
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(Champion and Minzner, 1963). Some typical data for higher heights is 
given in Fig. 3. In a later section, we consider methods of obtaining 
the experimental data on which such models are based. 


DENSITY ( kg/m 3 ) 



FIG. 3. THEORETICAL MODELS OF HARRIS AND PRIESTER 
(l962a,b). Vertical distribution of density and 
temperature for high solar activity (10 cm solar 
flux, S = 250 ) at noon (l) and midnight (2); and 
for low solar activity (s = 70) at noon (3) and 
midnight (4). These models assume the pressure, 
composition, and temperature at 120 km to be 
invariant . 


4. Dissociation and Diffusive Separation 

Above the mesopause the composition of the neutral atmosphere varies 
with height. This is associated with the dissociation of molecular oxy- 
gen and the occurrence of diffusive separation above the turbopause. 

The dissociation of oxygen has been studied in detail by Nicolet and 
Mange (1963). They point out that even though 0 molecules are dissoci- 
ated by radiation in the Schumann -Runge continuum (wavelength < 1750 S), 
the lifetime of an 0 molecule above 100 km is of order days, the actual 

Ct 

lifetime being very dependent on the intensity of the solar radiation. 
However, the 0 concentration in the thermosphere is replenished by 

Gi 
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diffusion in a time shorter than this, so that the vertical distribution 
of should more closely correspond to the hydrostatic equation (l-2) 
than to a condition of photochemical equilibrium. 

The N molecule, unlike 0 , has a very small cross-section for photo- 
dissociation and most of the nitrogen in the upper atmosphere is thought 
to be in molecular form. Such atomic nitrogen as does exist is produced 
by other photochemical processes, such as those involved in the production 
of ionization. Its abundance is linked with that of oxides of nitrogen. 

We may note in passing that nitric oxide, NO, is thought to exist in small 
quantities in the mesosphere and lower thermosphere. Probably at no level 
does NO represent more than one part in 10^ of the total concentration, 
yet its presence may be significant for the formation of the D region 
(Nicolet and Aikin, I960). 

The question of diffusive separation of atmospheric gases has been of 
interest ever since the composite nature of air was recognized in the 
eighteenth century. It was once thought that diffusion separation would 
take effect just above the ground, but it is now clear that major con- 
stituents are well mixed up to at least 100 km. Evidence is accumulating 
that diffusive separation begins near this altitude. 

From observations with rocket-borne mass spectrometers, Meadows and 
Townsend (1958) found molecular nitrogen and argon to be diffusively 
separated at heights above 115 km. The situation concerning the gases 
N , 0 , and 0 remained obscure, partly because Ahe chemical activity of 
atomic oxygen creates difficulties in mass-spectrometer sampling. The 
loss of atomic oxygen by chemical reactions within the apparatus and con- 
tamination by gases from the vehicle have been major problems. Some 
measurements have been published which show only a slow upward decrease 
of the [Ngl/fo] ratio at heights up to 210 km (Pokhunkov, 1963). However, 
the results of Schaefer and Brown (1964) have found a rapid upward decrease 
of the [0 ]/[0] ratio above 100 km, which indicates the existence of 
diffusive equilibrium even at this altitude. Recent measurements by Nier 
et al (1964), shown in Fig. 4, also indicate that diffusive separation of 
O, O and N occurs as low as 110 km. 
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FIG. 4. RATIOS OF THE CONCENTRATIONS OF THE MAJOR 
NEUTRAL CONSTITUENTS OF THE UPPER ATMOSPHERE 
(after Nier et al , 1964). 


We note that according to these data the level at which the concentra- 
tions of O and N are equal is near 200 km. Further measurements such as 
these should eventually resolve the discrepancies between existing models. 
For instance, Harris and Priester (1962) place this level at 190 - 300 km 
(depending on temperature) although Norton, Van Zandt and Denison (1963) 
place it at 130 km. 

The neutral composition can also be determined from observations of 
the solar spectrum at a series of heights, such as those of Hinteregger 
and Watanabe (1962). It is necessary to know the ultraviolet absorption 
cross-sections, as a function of wavelength, for each gas. This is essen- 
tially a trial -and -error method and has been applied by Norton, Van Zandt 
and Denison (1963) although their ratio of [n ]/[o] is considerably 
smaller than that measured by Nier et al , as we have noted. 

The question of diffusive separation is linked to that of turbulence, 
which is believed to be absent at heights above about 100 km. Blamont and 
de Jager (l96l) observed the rate of expansion of sodium vapor trails 
released from rockets. Below 102 km the trails were rapidly distorted by 
turbulence, but above this level they expanded more slowly, at a rate 
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consistent with molecular diffusion. A similar result was found by 
Rosenberg (1963) whose results place the turbopause at 110 to 115 km. 

The theory of mixing and diffusive separation is incomplete. A number 
of different approaches have been compared by Mange (1957). Most of these 
define a "criterion" for estimating the time taken for a minor constituent 
to attain a diffusive equilibrium, starting from an initial nonequilibrium 
(or mixed) condition. This time is in general of order (scale height) / 
(diffusion coefficient). If it is less than some time-constant character- 
istic of the process responsible for mixing, then mixing will not be 
maintained. Naturally, the problem is more complicated if photochemical 
processes are operative as well as diffusion. 

More recently, Hines (1963) discusses the problem of mixing of major 
constituents, with special reference to turbulence. If mixing were sud- 
denly to cease, the gases would separate with a relative velocity which 
is found by balancing the resistance to motion, due to collisions, against 
the "buoyancy" forces which tend to separate the gases. The power 
dissipation of the collisional forces can than be found: it is of order 

10~ 9 watts m 9 for a mixture of 0^ and N^. Since Hines estimates the 
power available from turbulence to be a few times larger than this, between 
90 and 100 km, he concludes that turbulence can just maintain mixing of 

0 and N up to 100 km, but not at greater heights. The precise level of 
2 2 

the turbopause presumably also depends on the strength of mechanism causing 
the turbulence, which may be gravity waves. Hence, changes in atmospheric 
motions below the turbopause might result in changes in composition at 
higher levels. A more complete discussion of turbulence, and of the per- 
tinence of parameters such as Reynolds and Richardson numbers, is given 
by Hines in the review cited above. 

Unless some other mechanism for preventing separation is effective at 
greater heights, it seems that diffusive separation of atmospheric gases 
must be established near the turbopause, or at least w 7 ithin one or two 
scale heights above it. As previously mentioned, the most recent avail- 
able experimental data do indeed indicate that diffusive separation begins 
at heights near 110 km. 
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5, Thermal Balance 


The heat balance in the thermosphere depends on the following pro- 
cesses: 

Production (q): (a) photochemical processes, such as photoionization, 

dissociation , recombination, attachment, (b) corpuscular ionization. 

(c) atmospheric absorption of gravity waves, hydromagnetic waves. 

(d) joule heating by electric currents. 

Loss (l): Radiation at wavelengths for which the thermosphere is 

optically thin. These include the visible airglow lines which are 
transmitted to the ground, and the infrared line of atomic oxygen 
at 63^. 

Conductive Transport: This is written as the divergence of a heat flux 

$ and thus resembles diffusion. 

These processes may be included in an equation of continuity. Fol- 
lowing Nicolet ( 1961 ) , we find it most convenient to write the equation 
in terms of temperature T, and make each item represent a change of 

energy per unit volume per unit time. Let c be the specific heat per 

P 

unit mass at constant pressure. Then, if only vertical transport is con- 
sidered , 


pc P 



" 5h 


(i-io) 


Before considering the form of the terms in this continuity equation, we 
note that, at the lower boundary of the thermosphere, the temperature 
minimum at the mesopause at 80 to 85 km, acts as a barrier to conduction. 
Thus, the thermosphere and the mesosphere can exchange heat only by radia- 
tion (or by forced convection driven by large-scale wind systems). The 
thermal regime in the mesosphere depends on the presence of ozone, which 
absorbs a wide band of solar ultraviolet radiation not significantly 
absorbed above the mesopause. 

Solar ultraviolet heating is thought to make a major contribution to 
the production term Q. The equations are akin to those describing the 
rate of photoionization which we will investigate later. We write 


Q(h) 


J l(A,h) £ e i o' i (A)n i (h) 


dA 


( 1 - 11 ) 
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in which the summation refers to different atmospheric gases and the 

formula is integrated with respect to ultraviolet wavelength A. 0\ 

_ _ . th , ,, 

represents an absorption cross-section for the 1 gas, and e the 
amount of heat liberated per unit of radiation absorbed. At any height 
the intensity of radiation l(A,h) is related to optical depth T(A,h) 
and the intensity l(A,°°) incident at the top of the atmosphere by 

ln[l(A,h)/l(A,°°) ] = T = sec X ^ H 1 (l-12) 

where % is the solar zenith angle, assumed not near 90 deg. At heights 

above about 250 km, t « 1 for all wavelengths. In this case, the total 
flux, I = / I ( A ,oo)dA, together with a total gas concentration n(h) and 
mean values F and a, can be used to replace the integration in (i-ll) 
which then reduces to 

Q(h) = I w ? a n(h) e T (l-13) 

The rate of heat input Q(h) is greatest at altitudes between 100 km 

and 200 km, where t(A) ^ 1 for the wavelengths involved. Under these 

circumstances the full equation (i-ll) must be. used, as the approximation 
(l-13) is inaccurate. Even the more detailed equation neglects two further 
complications connected with photoionization heating. One is that ejected 
pho t oel ec t rori s may travel considerable distances before giving up their 
kinetic energy, so that some of the energy may appear in remote locations. 
The other is that part of the energy is stored in the ionization, and 
does not appear as heat until the ionization recombines. This implies a 
time lag which may be of the order of hours at 300 km. Other sources of 
heat, previously mentioned, may be important, especially at times of 
magnetic disturbance. 

The "radiative loss" term, which is given by Bates (1952), is gener- 
ally regarded as having little effect above 150 km altitude and is omitted 
in many published solutions of the” continuity equation. 
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If only vertical variations are considered in the conduction term, 
the heat flux $ in the continuity equation (i-io) may be written 


0 = 


-AT 


1/2 dr 
5h 


(1-14) 


The product AT represents the thermal conductivity, A being a con- 

stant given by kinetic theory. According to Nicolet (l96l), its value 
in (joules m sec (°k) ' ) is 0.21 for atomic hydrogen, 0.036 for 

atomic oxygen and 0.018 for N 2 or 0 2 . 

The conditions at great heights in the thermosphere are largely 
governed by the conduction term, since the three other terms in the con- 
tinuity equation all decrease upward as the atmospheric density decreases 
If, as is generally supposed, the flux $ tends to zero at great heights 
then 5r/5h 0 by (l-14) so that the upper thermosphere is isothermal. 

However, there exist complications (to be discussed later in Sec. II -5) 
concerning the thermal disequilibrium between ions, electrons and neutral 
particles . 

To justify the assumption that $ -> 0 at great heights, we note that 
the outermost atmosphere is fully ionized and its thermal conductivity is 
large parallel to the magnetic field but extremely small across it. The 
flow of heat across the magnetopause (which effectively marks the termi- 
nation of the atmosphere), is then very small except when hot plasma from 
the sun actually pentrates into the geomagnetic field. Hence, the sug- 
gestion that the earth’s outermost atmosphere is at the temperature of 

5 

the surrounding solar corona, perhaps 10 °K, is not generally accepted 

(Chapman, 1959). We shall use the term ’’limiting (thermospheric) tempera 

ture" for the temperature in the isothermal region, although other terms 

such as ’’thermopause temperature" and "exospheric temperature" have been 

used. By solving the equation for equilibrium conditions, (such as 

approximately hold in the afternoon), Hunt and Van Zandt (l96l) obtain 

-2 -1 

an estimate of ~1 erg cm sec for the heat produced by ionizing radia- 
tion in the ionospheric F region, which corresponds to the upper part 
of the thermosphere. 
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If time variations are considered, it must be noted that the properties 
of the atmosphere at any given height h will vary on account of thermal 
expansion or contraction. Hence, it is convenient to "follow the motion’ 1 
of each "cell" of air, thereby applying the principles described in Sec, 

1-3 and formulate the equations in terms of reduced height z or pressure 
p, instead of real height h. Both Cummack (1962) and Harris and Priester 
(1962) use an approach of this type. The actual height and time varia- 
tions of thermospheric temperature are not known in sufficient detail to 
determine all the parameters of the models, but observations of satellite 
drag, in particular, enable good estimates of the limiting thermosphere 
temperature to be made (jacchia, I960). The data indicate that the limit- 
ing temperature is greatest about 14 hours local time, whereas the com- 
puted models give a later maximum, at about 17 hours. This discrepancy 
led Harris and Priester to postulate that a substantial fraction of the 
heating is due to a "corpuscular" source, which peaks at about 09 hours 
local time. This question must still be regarded as open. After sunset, 
photoionization heating is absent and conduction appears to equalize the 
neutral ion and electron temperatures above 200 km within two hours or so. 

Some graphs of scale height (or temperature) against height, obtained 
by combining data from several satellites, often show changes of slope 
around 200 km. These have been interpreted as "temperature inversion, " 
or layers in which T decreases upward (Priester and Martin, I960), and 
attributed to the fact that the rate of absorption of solar ultraviolet 
energy has a maximum in the lower F region, below 200 km. However, the 
rate of temperature increase due to heating depends on the heat production 
per particle which, by Eq . (l-13), is Q/n = I € cF e T and is therefore 
greatest at the top of the atmosphere where T = 0. Even if some other 
heating process should operate so as to produce a greater rate of tempera- 
ture rise below 200 km, the high thermal conductivity at 200 km and above 
would oppose the formation of an inversion, and doubts have been expressed 
as to its reality. 
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6. The Exosphere 


The exosphere is, by definition, the region in which the mean free 
path exceeds the scale height. Using data contained in the U.S. Standard 
Atmosphere, we may estimate this level to be at about 590 km. A molecule 
moving upward at the base of the exosphere is unlikely to make any colli- 
sion above this level, and thus moves in a ballistic orbit under the 
influence of gravity. This orbit will be elliptical if the upward velocity 

of the molecule is less than the velocity of escape v from the earth's 

E 

gravitational field, which is about 11.4 km/s. Should a molecule acquire 
a velocity exceeding 11.4 km/s as the result of collisions, its orbit will 
be hyperbolic and it will escape from the atmosphere unless it suffers 
further collisions. Some molecules will possess elliptical "satellite” 
orbits which do not intersect the base of the exosphere, in which the 
probability of colliding with other molecules is small. The computation 
of the radial distribution of gas concentration n(r) is extremely com- 
plicated, and has been the subject of some controversy (see Chamberlain, 
1963, and references contained therein). 

If we define an "escape temperature" T by the simple equation 

9 E 
l/2 mv^ = 3/2 kT , we find T, to be 5200 °K, 21,000 °K, and 84,000 °K 
EE E 

for H, He, and 0 atoms respectively. At sunspot maximum, the temperature 
of the exosphere may exceed 2000 °K and there will be a steady loss of 
hydrogen atoms at the upper end of the energy distribution. Helium will 
also be lost, but much more slowly. This gas is liberated at a steady 
rate by radioactive decay processes in the earth’s crust, at a rate which 
can be estimated from geological data. Since the rate of escape depends 
on the temperature of the upper atmosphere, it was thought that this 
temperature could be estimated from the observed abundance. However, it 
seems that the problem is too complex to yield useful information on this 
subject (Chamberlain, 1963). 

Hydrogen is generated by dissociation of water vapor and methane in 
the upper atmosphere, and then escapes from the top of the exosphere. 

In a detailed discussion of the hydrogen and helium content of the atmo- 
sphere, Kockarts and Nicolet (1962) estimate the upward flux of hydrogen 

7 -2 -1 

atoms at 100 km to be 2.5 X 10 cm sec . It may then be deduced that 
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the average lifetime for hydrogen atoms in the atmosphere is a few days, 
as compared to a few years for helium atoms. 

Great interest has been aroused by spectrometric observations which 
are connected with the hydrogen content of the atmosphere. Purcell and 
Tousey (i960) observed the solar Lyman a line from a rocket above 130 km. 
They found it to be broad (~1 X) with a narrow absorption core, which 
could be attributed to a column of hydrogen atoms, 3 x 10 (m ). The 
width of the absorption core corresponded to a temperature of order 
1000 °K, more likely to be associated with the terrestrial atmosphere 
than the interplanetary medium. Moreover, the absence of a Doppler shift 
suggested that the hydrogen is stationary with respect to the earth. 
Diffuse Lyman 0i radiation has also been detected in the night sky by 
rockets and it is believed to be solar in origin, and scattered by hydro- 
gen in the outer atmosphere. However, the amount of hydrogen required 
to account for the observed intensity is several times greater than the 
content quoted above, deduced from daytime observations. This difficulty 
has been well discussed by Donahue and Thomas (l963a,b) who suggest that 
it could be explained by the presence of a distant "cloud" of hydrogen 
("geocoma") associated with the earth, or by a large day-night variation 
of the hydrogen concentration. However, Hanson and Patterson (1963) con- 
sider that a lateral flow of hydrogen would be efficient in preventing 
such a diurnal asymmetry. 

7 . Experimental Techniques 

There are a number of experimental techniques available for the study 
of the neutral atmosphere. Below about 30 km it is possible to measure 
pressure, temperature and density with relatively inexpensive balloons. 

A "searchlight" probing technique has permitted ground-based measurements 
as high as 67 km (Elterman, 1951, 1954). In this method, a powerful 
searchlight beam which may be modulated to discriminate against the steady 
background radiation from the night sky, is directed upward. Another 
remote photo-sensitive detector then scans the searchlight beam to deter- 
mine the intensity of the scattered light as a function of altitude. When 
only Rayleigh scattering is observed, as appears to be the case above 
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about 10 km, the received intensity is proportional to the concentration 
of air molecules. Furthermore, when the equation of hydrostatic equili- 
brium is introduced, the temperature and pressure profiles may be deduced 
from the measured concentration profile. 

On one occasion, Elterman has obtained temperature profiles from rocket, 
sound propagation and balloon measurements nearly simultaneously with the 
searchlight experiment and all were in satisfactory agreement. His mea- 
surements have found that below 45 km there is little diurnal or seasonal 
temperature variation, but near the mesopeak at about 50 km and above, 
day-to-day fluctuations as well as diurnal and seasonal variations exist 
of 20 to 40 °K. The invention of the laser, since these experiments were 
performed with its attendant narrow bandwidth, offers the possibility of 
extending the measurements to even greater altitudes. 

Another very fruitful method for study of the atmosphere below about 
90 km has been the measurement of sound delay times from grenade explosions 
at high altitudes. As a rocket climbs to 100 km or more, the grenades are 
ejected and exploded every few kilometers. The delay times are recorded 
at a network of ground stations and from the records it is possible to 
deduce not only the temperature profile, accurate to within ±3 °K (Stroud 
etal, I960), but also the horizontal winds in the mesosphere. A sim- 
plified presentation of the theory and number of other references are 
given in the paper of Stroud et al . Winds in the mesosphere show little 
or no diurnal variation but there is a strong seasonal variation in the 
arctic. Strong westerly winds (100 to 150 m/sec) in the winter--the 
so-called "polar jet streams" --change to more moderate easterly winds in 
the summer. 

Most of the other methods require rockets to place the instruments 
within the environment it is intended to measure. For example, an experi- 
ment has been designed in which a collapsed sphere is launched to a high 
altitude and then ejected and inflated after the rocket power has been 
exhausted and the atmospheric drag has reached a tolerable magnitude 
(Bartman et al , 1956). The large area/mass ratio allows the drag force 
to exert measurable accelerations even in the very tenuous upper atmo- 
sphere. The drag force may be written as 
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(1-15) 


» = 5 PV A C D 

where p is the air density, V is the sphere velocity, A is the cross- 
sectional area and is a dimensionless drag coefficient. Although 

the value of is a rather complicated function of how the air molecules 

arc reflected from the surface of the sphere, Faucher et al ( 1963) believe 
that it can be calculated to within a few percent. 

After the sphere has been erected and vibrations have damped out, a 
sensitive accelerometer is used to measure the drag force and then the 
density is calculated from Eq^ (l-15). In a recent rocket flight reach- 
ing an altitude of about 250 km, Faucher et al detected measurable atmo- 
spheric drag below 135 km an d have computed density values below this 
level. The results are not only useful in themselves, but may be com- 
pared with rocket data based on pressure gauge measurements which extend 
to even higher altitudes. 

Lagow et al (1959) have compared the density and pressure profiles 
obtained in the arctic and in mid-latitudes from a number of rocket 
flights. A variety of pressure gauges were used in the measurements and 
from these, both the ambient pressure and the gas density were calculated 
up to an altitude of 210 km. Their results show rather large differences 
in scale height (and therefore temperature) between White Sands, New 
Mexico and Fort Churchill, Canada, on summer days. At high latitudes the 
density and scale height values appear too large to match satellite drag 
data. Sterne (1958) indicates reservations about the validity of the 
pressure measurements and also Johnson (i960) has pointed out the diffi- 
culty in maintaining large latitudinal pressure and temperature differences 
above 200 km. 

Following a suggestion by Bates (1950), a number of experiments have 
been performed in which sodium and other contaminants have been released 
at high altitudes. When the gas is released all along the rocket trajec- 
tory, a trail develops which becomes highly distorted as time proceeds 
owing to wind shears and molecular diffusion in the upper atmosphere. The 
sodium trail observed at twilight is believed to be due to resonant scat- 
tering of the 5580 $ sunlight and it is rapidly extinguished following 
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sunset at the appropriate altitude (Bedinger et al , 1958). However, 
even at night a brief trail is observed with the energy for the photo- 
emission apparently derived from chemical reactions with atmospheric 
constituents . 

At altitudes above 200 km, different methods must be employed for 
the measurement of atmospheric density. A "ribbon microphone" in which 
the incident flux of neutral molecules is modulated by a rotating shutter 
has been described by Sharp, Hanson and McKibbin (1962). The resulting 
ribbon displacements of only a few Angstroms produce an alternating 
voltage which is amplified with sufficient accuracy to obtain density 
measurements (proportional to the pressure oscillations) at an altitude 
of 550 km. Their results are consistent with values deduced from satel- 
lite drag observations and have the outstanding advantage of being 
nearly instantaneous measurements. 

Since 1957, it has been possible to measure atmospheric density from 
a careful measurement of the rate of decrease of the orbital period of 
a satellite. Very important advances have been made in our knowledge of 
the upper atmosphere from these observations. Although the gas density 
at 500 km is about 10 10 of that at the ground, appreciable drag is 
experienced by most satellites at this altitude and useful results have 
been obtained between 200 and 1100 km (jacchia, I960). The method of 
analysis is particularly simple for circular orbits, as will be shown 
below, but such an orbit is very difficult to achieve in practice, and 
the gas density is not spherically symmetric in any case. The accurate 
equations applicable to elliptical orbits have been derived by Sterne 
(1958) and others, and King-Hele (1963) has given a relatively simple 
and useful derivation of the most important terms. 

By equating the centrifugal and centripetal forces acting on a satel- 
lite in a circular orbit about the earth, we may obtain the relation 

2 -11 / 2 / 2 \ 
GM = rv , where G - gravitational constant, 6.67 X 10 (newton m /kg 

E +24 

M = mass of earth, 5.97 X 10 (kg), r = geocentric radius of the orbit 

E 

and v = velocity of the satellite. Then, we find that the kinetic energy 
is given by 
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(1-16) 


KE = 77 m v 2 = GM m /2r 
2 s E s 

where is the satellite mass. From first principles we know that 

the potential energy is 

PE = - GM m / r 
E s' 

We should note that the magnitude of the PE is just twice that of the 
KE at any radius and that their derivatives with respect to r are also 
in the ratio of 2:1. 

When atmospheric drag is included, the rate at which satellite energy 
UAost is Dv, where D is the drag force given in (l-15), This may 
be equated to the rate change of total energy: 

_Dv = ii (KE + PE) = £ ( RE + PE ) if ( l - 18 ) 

Solving for the rate of change of altitude, 


(1-17) 
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(1-19) 


For small values of p, the satellite will gradually spiral deeper into 
the atmosphere, gaining kinetic energy while losing potential energy 
twice as fast, in what may be termed a ’’quasi ”Circular ,f orbit. 

We observe the paradoxical situation in which the drag force acts in 
the direction opposite to the velocity vector, yet the velocity of the 
satellite continually increases. This may be understood as follows: the 

impact of a molecule at one point in the orbit does reduce the KE and 
thus the satellite velocity by a very small increment. This puts the 
satellite in a slightly elliptical orbit with the point of impact at 
apogee. When the satellite reaches perigee half an orbit later, the 
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velocity will have increased to a value greater than prior to impact at 
the expense of a reduction in PE. For the quasi -circular orbit, the 
variations of KE, PE and r are all continuous, of course. 

The measurements of atmospheric density rely upon the observation of 
the average rate change of satellite period, usually obtained from radio 
tracking stations. From Kepler r s third law, 


T 


2 



( 1 - 20 ) 


we obtain 



Again using 


GM = rv 


2 


and Eq . ( I -19) , 


dT 

dt 



( 1 - 21 ) 


( 1 - 22 ) 


The value of is usually considered to be about 2.2 for satellite 

applications. Equation (l-22) can then be used to calculate the density 

at the satellite altitude. The analogous equations for elliptical orbits 

l/2 

will not be derived, but it develops that the quantity pH is the 

parameter now determined by drag analysis. The values of p deduced 
when H is assumed are most accurate about one-half scale height above 
the satellite perigee. 
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II. measurement of ionospheric parameters 


1 . Introduction 

In the preceding section, the vertical distribution of the nonionized 
molecules and atoms in the earth’s atmosphere has been discussed. This 
is a very important preliminary to the discussion of the ionized portion 
of the atmosphere because the neutral molecules are obviously the parents 
of the charged particles to be discussed here. Our knowledge of the 
charge distribution has tended to precede that of the neutral composition 
principally because more simple measurement techniques are available.- 
These measurements, in turn, have contributed to the study of the neutral 
composition. 

Without discussing the theory of the formation of the ionosphere for 
the moment, in this section we will want to consider the various means 
available for the measurement of the pertinent ionospheric parameters. 
Principally, these are the electron ’’density” (or more precisely, con- 
centration), the ion concentration and composition, the ion and electron 
temperatures and the movement of these charges. 

The presence of the ionosphere was detected long before modern radio 
methods were developed. Over eighty years ago Balfour Stewart suggested 
the existence of a conducting layer in the upper atmosphere within which 
large currents were generated by the tidal oscillations of the atmosphere. 
These currents explained the observed magnetic deflections at ground level 
which wore found to vary with solar activity and geographic latitude. 

Then, after Marconi had demonstrated the feasibility of trans-Atlantic 
radio communication in 1901, Kennelly and Heaviside suggested independently 
that a reflecting or refracting layer must exist such that the radio waves 
were bent” around the curve of the earth. Direct radio investigations 
of the ionosphere may be dated from the pioneering work of Tuve and Breit 
(l925) and that of Appleton and his colleagues (1925). 

2. Determination of Electron Density by Sounding 

After forty years, the pulse sounding technique of Breit and Tuve is 
still a basic tool of ionospheric research. A sounder is a type of radar 
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which is capable of obtaining echoes from the ionosphere over a wide 
range of operating frequencies. The model C-4 ionosonde developed by 
the National Bureau of Standards has been in use at several stations for 
a number of years and is representative of the many designs now in use. 

In a typical mode of operation, the sounder is swept from 1 - 25 Mc/s in 
13.5 sec, using a pulse repetition frequency of about 100 per second and 
a peak power of up to 30 kilowatts. Newer types of ionosonde have been 
developed, which provide better resolution and greater versatility in 
operation , 

In a conventional recording system, the echo received from the iono- 
sphere is used to modulate the intensity of a spot of light on an elec- 
tronic time base. Distance along the time base represents the "time of 
flight" of the radio pulse which, if divided by the free space velocity, 
gives the equivalent path length, and this is twice the "virtual height 
h’ of the reflection point in the ionosphere. In the recorder, a photo- 
graphic film is moved at right angles to the time base as the frequency 
is varied, so that the spot of light traces a graph of virtual height h' 
against radio frequency. Range and frequency calibration markers are 
usually inserted automatically. These recordings are commonly known as 
ionograms or h'(f) curves. 

The pioneer ionospheric observatory is situated at Slough in England, 
where an ionosonde has been operated regularly since 1932. Other stations 
were opened before the war at Washington, D.C. , Watheroo in Western 
Australia, and at Huancayo on the magnetic equator in Peru. Since World 
War II many other stations have been operating and during the International 
Geophysical Year, in 1957-58, the total number was about 150. Normally, 
each station makes a sounding once an hour (sometimes once every 15 min 
or less) and the total number of ionograms now in existence amounts to 
some tens of millions. 

To understand the form of the h'(f) curves, it is necessary to con- 
sider briefly the equations for the propagation of waves in the ionosphere. 
Let N be the electron concentration (electron density), m and -e the 
mass and electronic charge, and £ q the permittivity of free space. In 
the absence of a magnetic field, the phase refractive index p and group 
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refractive index p 1 , for waves of frequency f are related by the 
equation pp ! = 1, and 


2 , Ne 2 

^ = 1 ' ~ 2 

4tc me f 
o 


(n-i) 


Note that p- — 1 — M- ' (which is usually true even in the presence of a 
magnetic field). The wave is reflected from a level where p = 0 , or in 
other words, where the frequency f is equal to the local value of the 
plasma frequency f given (in c/s for N in electrons/m 3 ) by 

f N = ~ 9 /N (lI-2) 


In the presence of the magnetic field, the ionosphere is a doubly 
refracting medium, and two modes of propagation exist, for which the 
names "ordinary" and "extraordinary" are taken from crystal optics. The 
complete Appleton-Hartree equation replaces the simplified equation 
(il-l) and it is no longer true that pp' = 1 (Ratcliffe, 1959). The 
equation contains the frequency of gyration of electrons in the geomagnetic 
field B,' given by f = eB/2jtm, and the magnetic dip angle. It thus 
depends on the location of the observing station. For a more accurate 
treatment, the effect of collisions between electrons and heavy particles 
should also be included. An idealized h*(f) curve is sketched in Fig. 

5; its most obvious features are the "cusps'* which mark the "critical 
frequencies " of the E and F layers. 

Both "ordinary" and "extraordinary" traces are sketched in the figure. 
By day, the F layer is often divided and shows an "FI cusp" at a fre- 
quency lower than the ordinary mode penetration frequency which is nor- 
mally known as the "F2 critical frequency" and denoted by f F2. Waves 

o 

of frequencies exceeding f^F2 cannot be reflected at vertical incidence, 

but penetrate the ionosphere completely. The symbols f FI and f FI 

o x 

are used for the FI cusp frequencies, if present, and f E and f E 

o x 

for the E layer cusp frequencies. At times additional stratifications 
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FIG. 5. AN IDEALIZED IONOGRAM, SHOWING VIRTUAL HEIGHT VS 
FREQUENCY. The solid and broken curves apply to the 
ordinary and extraordinary modes, respectively. The 
dotted curve represents the plasma frequency -real 
height profile necessary to produce the virtual height 
curves . 

are observed; in fact, ionograms reveal a wide variety of phenomena, some 
of which are touched upon in subsequent sections. Also shown as the 
dotted curve in Fig. 5 is the plasma frequency profile which could have 
produced the virtual height curves. 

The critical frequency f Q F2 (some authors omit the subscript "o") 
is the most widely studied property of the ionosphere. It is related to 
the greatest electron density N^ within each layer by the standard 
equation (lI-2). But we note that (except for NJF2) these values may 
not represent actual peaks (maxima) of the N(h) distribution, but may 
just be points of inflection. The virtual heights at these points are 
large, and are actually infinite if dN/dh = 0. This, of course, is 
always the case at the F2 peak, but the E and FI cusps may correspond 
to quite small "ledges" in the N(h) profile, such as is sketched in 
Fig. 5. 
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We note that the separation of the ordinary and extraordinary cusps 
can be computed from the Appleton-Hartree equation. It is given by the 
equations 


f 

o 



X 


f f 
x H 


1/2. 


f 

X 


i f H 


(11-3) 


The approximation is good if the critical frequencies are much greater 

than f , as is generally the case for the F layer. (A typical mid- 
H 

latitude value of f is 1.4 Mc/s). Under certain circumstances, a 

H . . 

third trace is seen on ionograms , notably at high latitudes. The critical 
frequency f of this so-called z-trace is approximately ^f^ - 2 fjj) * 

A discussion of this phenomenon is given by Ratcliff e (1959). 

The virtual heights observed at and near the critical frequencies 
have little physical meaning, but the real heights h^E , h^Fl and h^F2 
are important quantities. The virtual height (or "group height") h 
and the real height h at which a given frequency f is reflected are 
connected by an equation that depends on the distribution of electron 
density below the reflection point; it is 



Since p f ^ 1, the group height h* ^ h. Much effort has been devoted 
to solving this equation to obtain N(h). This is accomplished by writing 
it as an integral with respect to plasma frequency instead of height 
(which assumes that N and f vary monotonically with height) 

h'(f) = f n'(f,N) df + h(0) (II-4) 

0 N 

where h(o) is the height of the base of the ionosphere, below which 
f N = o - 

A good review of the ways of inverting this equation has been given 
by Thomas (1959). The availability of electronic computers has made 
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obsolete the earlier methods in which a special form, such as a parabola 
was assumed for the N(h) profile, though these are of academic interest. 

The present-day methods of solving the integral equation (II-4) fall 
into two main classes, "lamination" and "polynomial." The former class, 
represented by the "matrix method" of Budden (l955) replaces the inte- 
gration in (lI-4) by a summation over a number of thin slabs, each cor- 
responding to a discrete interval of plasma frequency. For this purpose, 
some simplifying assumption is made, the commonest being that the gradient 
dh/df^ is constant within each slab so that the function h ( f N ) is 
represented by a series of linear segments. 

Suppose there are n slabs, and that h’(f.) is the virtual height 

^ t h 

for the frequency f reflected at the top of the i slab, at a real 

height h . We can then form a set of n + 1 equations 


b'(0) = h(0) 


h(f ) - h(o) A 

h’(V = h(0) + — i — — J ^'(f r f N ) df N 


h'(f 2 ) = h(0) 


h(f ) - h( 0 ) A h(f ) - h(f ) , f 2 

+ f, - 0 L A f 2’ f N^ df N + f„ - f J 


) df v 


2 1 f 


(II-5) 


h-cfj) = h(o) + y 

i=l 


h(f.) - h(f i _ 1 ) . i 


f. - f. . 
1 i-l 




i-1 


Additional equations are written up to j = n (note that we take 
f^_^ = 0 for i = l). By using the abbreviation 


M . . = (f. 
Ji 1 



“'('j'V d, N 


(II -6) 
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we can write the equations compactly in the form 


j 

h'(fj) - h'(0) = [h(f.) - hCf.^)] M ( II-7a) 

i=l 


or 

j 

= 2 (Ah i ) “ji (II -7b) 

i=l 


Note that Ah ! . is the virtual height at frequency fj measured from the 

bottom of the layer, Ah, is the real width of the i ^ slab, while M 

! ^ ’ ji 


is the average value of 
f 


th 


p’ witnm tne i slab for waves of frequency 

This equation can obviously be written in matrix form, Ah 1 = M • Ah , 

J -1 — 
whose solution is Ah = M • Ah' . From the width of these slabs, the 

real height profile is obtained as 


h j - h '(°) * J *»l 

i=l 


It is common practice to perform both of the last two steps at the same 

time with an additional matrix multiplication (Budden, 1955), 

It may be advantageous to use an unequally spaced set of frequencies 

f^. These can be so chosen that the summation reduces to a simple 

averaging of a number of values of h l , which can be carried out manually 

(Kelso, 1952). Schmerling (1958) has shown that this method can be 

adapted to include the magnetic field. Some authors do not explicitly 

include the term h(o) in their equations. For example, Thomas (1959) 

imposes an additional condition on the matrix elements which takes 

account of this term, but neglects the effect of ionization in the first 

lamination (i.e., below the height h(f^) at ^ he first frequency 

f is reflected), 
o / 
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The assumption that dh/d± N = constant within each slab may be varied. 
In a comparison of different methods, King (1957) advocates the assumption, 
dh/d(ln f ) = constant, whereas others have used a parabolic variation 
within each slab. 

In the "polynomial" method of Titheridge (l96l), it is assumed that 
the function h ( f N ) can be represented by a polynomial in f N> of the 
form 

n 

h - h(o) = h* = ^ a. f^ +1 (II-8) 

j=l 

First a set of sampling frequencies f , f f is chosen. Then 

X ^ XI 

the heights h* at which the plasma frequency takes the values ^ can 

be included in a matrix equation, h* = A • a, in which each element of 

the column matrix a is the coefficient appearing in Eq. (lI-8), and 

the element a is equal to f? +1 . Now the set of virtual heights at 
1J 1 

frequencies f can be written as h ' = B • a, in which the elements 
of B contain integrals with respect to frequency of the group refrac- 
tive index p' . Hence h* = ( AB ) • h*' . The elements of the matrix 
AB _1 depend on the frequencies f and the parameters of the geomagnetic 
field, and can be computed for a given location and a given set of sam- 
pling frequencies. Thus sets of values of h' and h are connected 
by a matrix equation, as in the lamination method. 

Whereas the lamination method represents the function h ( f N ) a 

series of linear segments, the polynomial method instead assumes that a 
smooth polynomial of degree (n + l) in f N can be drawn to pass through 
the n sampling points although the coefficients of this poly- 

nomial are not required to be known. For a given number of sampling 
points, the polynomial method gives greater accuracy than the lamination 
method. Titheridge (l96l) also describes a "modified polynomial method" 
of greater complexity, which is suitable for reducing ionograms containing 
a number of cusps. 


SEL-64-111 


34 



There are several difficulties to be faced in. the s reduction of h ! (f ) 
curves to N(h) profiles. One of these is the so-called ’Valley 
ambiguity. ,f All the methods of analysis that we have described give 
unambiguous results only if electron density N increases monotonically 
with height. If the N(h) profile is re-entrant so that the same value 
of N occurs at more than one height, then the reduction procedures do 
hot give a unique result. For instance, it is possible that N decreases 
with height immediately above the E layer to give a "valley" between 
the E and F layers. Indeed, this was generally thought to be the 
case until rocket observations proved otherwise. It is also possible, 
though less likely, that there could be a valley between the FI and 
F2 layers. 

Another difficulty arises from the fact that conventional lonosondes 
have a lower frequency limit, usually about 0.5 - 1 Mc/s, and record no 
echoes below this limit. Consequently, they obtain no information about 
the electron distribution below the E layer in the daytime or below 
the F layer at night. Special sounders have been built which operate 
down to 50 kilocycles (Watts, 1957), but this equipment is not in general 
use. The lack of data about these low frequencies can lead to errors of 
10 - 20 kilometers in the height corresponding to higher plasma frequencies. 

- Titheri dge (l959) has shown that if data from the extraordinary ray 
virtual height are included in the analysis, it is possible to correct 
for the low level ionization not included in the ionograms, and also to 
make some estimate of the ionization contained in the valleys. Although 
the method does not give a detailed distribution of electrons in the 
valley, Titheridge was able to show that on some occasions, only a very 
small valley existed between the E layer and F layer at night. Paul 
and Wright (1963) have given some results obtained with a very refined, 
method of ionogram reduction. It seems, however, that we shall have to 
wait for the improved records which should be obtained from the most 
advanced ionosondes before the full potentialities of this method can be 
realized. It should then be possible to calculate the real height at 
which a given electron density occurs to within about 1 kilometer, whereas 
with the older methods of reduction, in which the difficulties described 
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above have not been fully overcome, the errors might be of the order of 
ten or even twenty kilometers (Thomas, Haselgrove and Robbins, 1958). 

The successful launching of the Canadian satellite, "Alouette, 11 in 
September 1962, inaugurated the era of the "topside sounder." This 
satellite is in a nearly circular orbit at a height of just over one 
thousand kilometers, and carries a compact sweep-frequency ionosonde, 
operating on radio frequencies of 0.5 - 11.5 Mc/s, Whenever the satellite 
is activated, soundings are made at intervals of 20 sec, during which 
time the satellite travels about 100 km. The interpretation of the iono- 
grams obtained from the satellite and the reduction to N(h) profiles 
presents some interesting problems. The analysis is complicated by the 
fact that neither the electron density N nor the gradient dN/dh is 
zero at the sounder, as is the case with a conventional ionosonde on the 
ground. On the other hand, there is generally less structure in the top- 
side than in the bottomside of the ionosphere, so that there is less dif- 
ficulty with valleys. Nelms (1963) has developed a reduction program which 
is based on the original Budden matrix method, whereas Fitzenreiter and 
Blumle (1964) use a lamination method which takes dh/d(ln f ) = constant 
within each slab. Thomas, Long and Westover (1963) have used a polynomial 
method in which each real height is expressed as a polynomial in f - f , 
f being the plasma frequency at the vehicle. The topside records show 
a number of interesting magneto-ionic and plasma phenomena, but we shall 
not describe these in detail here (Calvert and Goe, 1963). 

3. Propagation Methods 

In contrast to the relative ease with which electron density measure- 

3 

ments are made above 90 km where the concentrations are greater than 10 
per cc, the region between 50 and 90 km presents formidable experimental 
difficulties. One useful method first described by Gardner and Pawsey 
(1953) relies on partial reflections from small discontinuities or 
gradients in the index of refraction in the lower ionosphere. The tech- 
nique involves the transmission of short rf pulses (2 to 3 Mc/s) and the 
reception of the reflected energy as a function height (time). The 
received energy must be separated into its two characteristic modes, which 
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will be right hand and left hand circular except near the equator, as it 
is the ratio of the amplitudes of the modes that is important. For either 
mode, the received amplitude is given by 



A oc R exp 



(n-9) 


where R is the power reflection coefficient and K is the attenuation 
per unit length. The upper limit of the integral is the reflection 
.height. The power ratio of the extraordinary mode to the ordinary mode 
is therefore 


A 

x 

A 

o 


R 

x 

r exp 

o 


-2 J <K X - K 0 ) 


dh 


( II -10 ) 


From the appropriate expressions for the index of refraction, it may 

be found that the ratio (r /r ) is a function of the electron collision 

v x o' 

frequency v, but not the electron concentration, while the differential 
absorption (j^ - & ) depends on a different function of N and v 
(Holt, 1963). If the collision frequency is determined independently, 
the appropriate functions can be calculated and Eq. (iI-lO) can be 
revised to give 


f 

ln(A /A ) - ln( R /R ) = -2 / (& - K ) 

K x o x' o' X o' 


dh 


f h 

= / f(v,h) • N(h) dh (il-ll) 

J 0 


In this equation (A /a ) is measured while (R /r ) and f(v.h) are 

v x o' ^ v x o' \ / 

calculated from the assumed collision Frequency profile. Then, the equa- 
tion may be differentiated with respect to height and solved for N(h). 


I . ±. r -_± 
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An example of these calculations and observations is shown in Fig. 6 
taken from Holt (1963). The ratio of the reflection coefficients rises 
with altitude owing to the increasing collision frequency. However, at 
the higher altitudes, the relative absorption of the extraordinary wave 
increases even faster, producing an amplitude variation which first 
increases, then drops off rapidly. Figure 6b shows the electron density 
profile calculated from these observations. At the lowest altitudes, 
the absorption is relatively small and the right hand side of (il-ll) 
is approximately zero. Here, (A^/A^) « (r^/r ) and the collision fre- 
quency may be calculated. 
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FIG. 6. (a.) OBSERVED AMPLITUDE RATIOS OF PARTIALLY 

REFLECTED EXTRAORDINARY AND ORDINARY WAVES (OPEN 
CIRCLES); THE BARS SHOW UPPER AND LOWER QUARTILES. 
The reflection coefficient ratio is calculated from 
an assumed collision frequency profile. (b.) THE 
DEDUCED ELECTRON DENSITY PROFILE. The measurements 
were made in August 1962, in Norway (after Holt, 
1963). 


An alternative technique first described by Fejer (1955) is based on 
radio wave interaction and permits the calculation of both the collision 
frequency and the electron density as a function of altitude. This 
method requires one transmitter (usually at a frequency of two to three 
Mc/s) to emit short pulses which are reflected from the ionosphere and 
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received at the ground. The received amplitude of this "wanted” signal 

is carefully measured. When a short interval (t ) exists prior to the 

o 

reception of every alternate wanted pulse, a second "disturbing” trans- 
mitter emits another short pulse at a different frequency and with as 
much power as possible. This situation is sketched in Fig. 7, In reality 
T , T^ and R may all be at one location, but they are shown separated 
here for clarity. 


LOG v (sec" 1 ) 

6 7 8 9 
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FIG. 7. ILLUSTRATING THE MEASUREMENT OF N AND v BY THE RADIO- 
WAVE INTERACTION TECHNIQUE. A pulse from the "disturbing” 
transmitter T^ interacts with alternate pulses from the 
"wanted” transmitter at the altitude h Q . A region of 

slightly enhanced collision frequency follows in the "wake" 
of the disturbing pulse. The N(h) and v(h) profile 
obtained by Barrington and Thrane (1962) is shown at the 
right. 


At the instant shown, the downcoming "wanted” wave is just meeting 

the upgoing "disturbing" wave at the height h « ct . The electro- 

o o 

magnetic wave comprising the pulse from T forces the electrons in the 
ionosphere to oscillate at its wave frequency, and this energy is con- 
verted to heat (e.g., random velocity) by collisions with the neutral 
particles. If a fraction G of the excess electron energy is lost at 
each collision, the electron will come back to its original temperature 
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in a time T ^ (Gv) 1 . Thus, the small excess electron temperature (AT) 
'will diminish behind the pulse in a distance Ah ~ cT as shown, where 
c is the velocity of light. We may expect v to depend on temperature, 
and in this way, the disturbing pulse generates a region of slightly higher 
collision frequency as it passes through the ionosphere. 

Now we wish to see what effect the increased v produces in the 
amplitude of the wanted wave. When the wave frequency exceeds v, an 
increase in v will increase the absorption and reduce the received 
amplitude. However, at the lowest altitudes, below about 70 km, a col- 
lision is likely to occur before the electromagnetic wave has fully 
accelerated the electron ( v > <jo) , and in this case the increased v 
reduces the energy loss and the absorption. The received amplitude of 
the wanted pulse which interacts with the disturbing pulse at altitudes 

below h is then compared with the next wanted pulse, for which no 
o 

disturbing pulse was sent. The difference in these amplitudes may then 

be related to v at h . As the time difference t is altered, the 

o o 

collision frequency at other altitudes is measured. 

After the v(h) profile has been determined, it is possible in 
principle to determine N(h) from the height attenuation of the wanted 
wave. The calculation tends to be rather inaccurate and "best-fitting” 
profiles are sometimes employed (Barrington and Thrane, 1962). 

In the last few years, the absorption of cosmic radio noise has led 
to another method of studying the lower ionosphere. An instrument called 
a riometer (for relative ionospheric opacity meter ) scans across a small 
portion of the radio frequency spectrum, perhaps a hundred kilocycles, 
and determines the minimum noise level observed (Little and Leinbach, 

1958, 1959), The scanning procedure is necessary to avoid strong inter- 
ference at various (usually numerous) points across the band. Variations 
in the minimum noise level with local time are first corrected for the 
variation of incident radio flux with sidereal time and then the remaining 
changes are usually attributed to ionospheric absorption (Lusignan, I960). 
Under normal circumstances at mid-latitudes, the D and T regions may 
each contribute about 1 db to the total absorption at wave frequencies 
between 20 and 30 Mc/s, 
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When the D region absorption can be separated or when it predominates 
during "disturbed ionospheric conditions/’ the absorption on several fre- 
quencies can be used to estimate the D region electron density profile 
(Parthasarathy et al , 1963). Since the absorption is proportional to 
the integral of (Nv) with respect to height, it is first necessary to 
assume the collision frequency profile. (Most recent work assumes that 
V oc T, from which an "effective collision frequency" may be defined.) 

Then various electron density profiles may be tried and adjusted to best 
fit the absorption data. 

One of the most widely used methods for obtaining electron density 
profiles is the cw propagation technique pioneered by Seddon (l953). The 
experiment is performed by firing a rocket equipped with two coherent 
transmitters in a trajectory as nearly vertical as possible. By Eq. (il-l) 
the index of refraction for radio waves is |i « 1 - in which 

the magnetic field has been neglected and a sufficiently high wave fre- 
quency is assumed. The constant = 40.3 in inks units. The two radio 

transmissions are received on the ground and their precise frequencies 
(or phases) compared in the following manner. The phase velocity of 

either frequency is v = c/u and the wavelength is A = v /f = A /u. 

p m p' o' ^ 

These quantities will both change with height as the electron concentration 
varies. The "electrical length" of the paths is given by 


L = 



( 11 - 12 ) 


The quantity L is simply the number of electrical wavelengths between 
the rocket and the receivers, and the time rate of change, dL/dt = L, 
is the Doppler shift. For simplicity let us imagine that two frequencies 
in the ratio f = 2f ^ , were used. If the electrical length at the higher 
frequency is divided by two and then compared with the lower frequency, 
we obtain 



2 





dh ( 11-13) 
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Since 


in which the integral extends along the ray path, assumed vertical, 

between the ground receiver and transmitter located at height h^. 

?\. = 77 A , , this becomes 

o 2 2 ol 



The time derivative of (ll-14) may be called the "differential Doppler 
f requ ency 


5 




_ K N dh 

3 JL_ v 

4 cf ][ dt 


(11-15) 


From this equation, the electron density profile N(h) may be deduced 
from the measured differential Doppler frequency. In practice, there are 
additional complexities introduced by the magnetic field and the horizontal 
component of the rocket velocity. However, the effect of these correction 
factors may be reduced by making the vertical velocity as large as possible. 

An analogous calculation can be made for the case of coherent trans- 
mitters in a satellite in a circular orbit above a spherically symmetric 
ionosphere; however, in this case only height-integrated electron density, 
or the electron content is obtained (Garriott and Bracewell, 196l). As 
before 


r S 

AL = — — — I N sec Y dh { 1 1 - 1 6 ) 

4 Cl l J o 

The obliquity of 
where )( is the 
the satellite in 


the ray path has been included by letting ds = sec dh, 
vertical angle of the ray. Now, the geometric range of 
wavelengths is r/A^ = 1 A ol / sec X dh . thus 
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(H-17) 


3 1 . 

■t —7— sec Y J N dh 

AL 4 cf l F 

R Aol " A oi sec X h s 

In the numerator, sec X should be evaluated at the height of the centroid 
of the electron density profile in the F region, while in the denominator 
it should be evaluated at about (h^/ 2 ), But since sec X is nearly con- 
stant with altitude, except at very large zenith angles, the values of 
sec X m ay be cancelled, revealing that both sides of (ll-17) should be 
constant with time during a satellite passage. Furthermore, their time 
derivatives 


5 

r/a 

o 


3K “s 

~2~ J N 


4f\ 

1 s 


dh 


(n- 18 ) 


arc in a constant ratio. Since 6 is the most frequently measured 
parameter and R and h g are known from the satellite orbit, the elec- 
tron content, f N dh, may be calculated. Again, more complicated pro- 
cedures are available to account for the lack of a symmetrical electron 
distribution (Garriott and de Mendon^a, 1963), but this simple derivation 
contains the essential elements. 

Electron content measurements have also been made based on an entirely 
different principle which requires the presence of the earth's magnetic 
field in the ionosphere. When an electromagnetic wave passes through an 
ionized medium in the presence of a magnetic field (as it does in the 
ionosphere) the polarization of the wave is observed to rotate as long 
as propagation is not too close to perpendicular to the magnetic field. 
This phenomenon is called Faraday rotation. The angle through which the 
polarization is rotated is 


- mi 


NH cos 0 sec )( dh 


(11-19) 
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where is a constant equal to 0.0298 in mks units, H = magnetic field 

intensity, 0 = angle between magnetic field and ray direction and the 
rotation angle fi is expressed in radians. Usually, the terms 
(H cos 6 sec x) are taken outside the integral sign as an approximation, 
since they vary little over the portion of the path in the ionosphere 
where N is significantly large. A measurement of fi in addition to 
knowledge of the pertinent geometric and magnetic field terms will 
therefore permit the electron content to be calculated. 

This method was applied to the polarization of moon echoes by Evans 
(1957) at Jodrell Bank and the diurnal variation of electron content 
calculated. The most satisfactory approach involved the transmission 
of two closely spaced frequencies from which the "differential rotation 
angle" could be observed. 

Even more extensive studies based upon the Faraday effect have been 
accomplished with satellite transmissions (Lawrence et al , 1963). There 
is a variety of analysis methods available usually based on either the 
change in rotation angle in a certain time interval or the rate of 
rotation during the passage. It is possible to improve the accuracy 
of the calculations by combining both Faraday and Doppler observations 
in a "hybrid" method (de Mendon^a and Garriott, 1963). From these 
investigations \ the variation of electron content with latitude, local 
time, season and solar activity has been investigated. 

One of the principal purposes for continuing the study of electron 
content is as an indicator of ionospheric temperature. When the shape 
of the layer is assumed, the product of N m and scale height (or tem- 
perature) is proportional to the electron content. Since can be 

determined readily, measured values of electron content may be related 
to the plasma temperature. 

4. Direct Measurements 

One of the more useful direct measurement instruments is the dumbbell- 
shaped electrostatic probe devised by Spencer and his associates (Spencer 
et al, 1962), The instrument is launched to a high altitude and then 
ejected from the rocket to avoid any local perturbations. One end of 


SEL-64-111 


44 


the device is forced to a negative potential by a sawtooth wave form and 
the resulting plasma currents are measured as in a Langmuir probe. The 
measurements are telemetered back to the ground where they may be inter- 
preted in terms of the local ion concentration and the electron tempera- 
ture. The results of five flights are compared by Brace et al (1963) 
and permit the comparison of high and middle latitudes, quiet and dis- 
turbed conditions and diurnal effects. 

The ion concentration profiles follow quite closely the results of 
ionogram analysis and the profiles deduced from simultaneous two-frequency 
radio propagation experiments which were discussed, in the previous section. 
However, the daytime electron temperatures are considerably larger than 
the temperatures deduced for the neutral atmosphere from satellite drag. 
This is not unexpected, since the electrons released by photoionization 
will have an excess energy of some 10 ev and a portion of this energy is 
eventually shared with the rest of the ambient electrons. The results 
of their electron temperature measurements are summarized in Fig, 8. 

The theoretical results of Hanson (1962) are very close to the quiet 
day profile shown here. At night the electrons have time to drop back 
to the ion and neutral temperature but under disturbed conditions or at 
high latitudes the electron temperatures remain above 2000 °K. 

This temperature nonequilibrium has been further explored by Nagy 
et al ( 1963) by measurements made with a spherical ion trap and a thin 
cylindrical probe. A rocket was launched near noon from a site in 
Florida, U.S.A., and both ion and electron temperatures were calculated. 
Between 200 km and 300 km the ion temperature was about 1400 °K while 
the electron temperature was about 2600 °K, An extensive survey of 
electron temperature, as a function of position and local time, has been 
obtained from a probe experiment on the UK-1 satellite "Ariel" (Willmore, 
Henderson, Boyd and Bowen, in press, 1964). The data refer to altitudes 
between 400 and 1200 km. 

Other direct measurement devices include the Gerdien condenser (use- 
ful in the lower ionosphere), rf impedance measurements at frequencies 
above the local plasma frequency, resonance probes operating near the 
plasma frequency, intermediate frequency measurements at VLF, and quasi- 
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D.c. instruments similar to Langmuir probes and ion traps. Most of these 
instruments have been operated on a number of occasions in the ionosphere 
and their results intercompared . The paper by Crawford and Mlodnosky 
(1964) reviews their operation and indicates the required assumptions in 
regard to electron and ion mobility in each frequency range. 
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FIG. 8. ELECTRON TEMPERATURE PROFILES 
OBTAINED ON FIVE ROCKET FLIGHTS FROM 
WALLOPS ISLAND, VA. (Lat. 38 °N) AND 
FORT CHURCHILL, CANADA (55 0 N). 
Variations with latitude, local time 
and magnetic activity all appear to 
be evident (after Brace et al , 1963) 


5 . Incoherent Scatter 

Conventional ionospheric sounding depends on the reflection of radio 
waves by an ionized gas, as described before. This reflection process 
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depends on the collective behavior of electrons, which can be described 
(at least approximately) in terms of the refractive indices of the 
Appleton -Hart ree equation. A second kind of echo, much weaker than the 
first, arises from scattering from irregularities or sharp gradients of 
ionization when the frequency of the radio waves exceeds the local plasma 
frequency. A third kind of echo, even weaker than the others, is due to 
the classical Thomson scattering of waves from individual electrons. 

This has become known as "incoherent scatter"; the power returned is 
directly proportional to the electron density. 

In his initial paper on this subject, Gordon (1958) discussed the 
conditions under which ionospheric electrons scatter incoherently, so 
that the return is of the third kind. He found that incoherent scatter- 
ing occurs if the electronic mean free path exceeds the radio wavelength 
and exceeds the scale size of irregularities and concluded that these 
conditions hold at meter wavelengths at heights above 100 km. Given a 
sufficiently powerful radar, the method could be used to measure electron 
densities throughout the F region and beyond, to a distance of one 
earth radius or more. 

Experiments first carried out by Bowles (l96l), with a 41 Mc/s radar, 
at 4 - 6 Mw peak power and a 100 psec pulse length, yielded "scatter 
profiles" of electron density, from about 100 to 700 km altitude. A 
high powered scatter radar has now been built at Jicamarca, on the 
magnetic equator in Peru. It operates at 50 Mc/s and has a maximum peak 
power of 5 Mw feeding a large array of dipoles over 20 acres in size. 

Such powerful equipment can be used for purposes other than incoherent 
scatter, such as investigation of ionospheric irregularities by coherent 
scatter, and even planetary radar. Examples of electron density profiles 
to 7000 km altitude have been given by Bowles (1963). Another very large 
installation is operated by Cornell University in Arecibo, Puerto Rico. 

Following the first observation, a number of theoretical papers on 
incoherent scatter were published (Fejer, 1960; Salpeter, 1960; Dougherty 
and Farley, I960). These were largely concerned with the broadening of 
the spectrum of the returned signal due to thermal motions of the scat- 
tering particles. It is found that the breadth of the spectrum depends 


47 


SEL-64-111 


on the ionic (rather than the electronic ) thermal velocity, and that the 
spectrum does not have a simple Gaussian form. Applications of this 
theory give values of temperature broadly consistent with other estimates. 

In fact, the scattering is not fully incoherent, but depends on weak 
statistical fluctuations of electron density; in this case, the scattering 

cross-section corresponds to half of the classical value, the latter 

4 / 2 2 4 , - 28 2 

being e /4^e me = 1 X 10 m . 

More recently, efforts have been made to extend the theory to include 
complications existing in the actual ionosphere. These are the inequality 
of electron and ion temperatures in the F region, the existence of light 
ions (He + and H + ) in the upper F region, and the geomagnetic field. 

The situation appears to be complicated because it may be difficult to 
distinguish between the effects of the first two factors, though magnetic 
field effects do not seem likely to be of importance in practice. Some 
of these matters are dealt with by Fejer (l96l), Hagfors (l96l), Farley, 
Dougherty and Barron (l96l), Salpeter (1963) and Moorcroft (1964), 

There are three ways of obtaining data about ionospheric temperatures 
from incoherent scatter data. First, the slope of a topside N(h) profile 
gives the sum of electron and ion temperatures, provided the mean molecu- 
lar mass M of the ions is known, because as we shall see in Sec. Ill 5 

d( In N)/dh = M + g/R(T e + T.j (II-20) 

This equation is not, however, applicable near the magnetic equator. Sec- 
ondly, the ratio (T /t. ) can be deduced from the echo spectrum, again 
provided the ion species is known. Thirdly, the scattering cross-section, 
which depends on T /t , can be deduced from the intensity of the echo if 
the electron density is known independently from ionosonde data. Obviously, 
this method requires accurate absolute calibration of the radar system. 

It is scarcely surprising that the results so far published involve 
certain assumptions, mostly about the ionic composition or height depen- 
dence of T and T . Evans (1962) and Pineo and Hynek (1962) find 
e i 

that T /t =1 at night but is greater by day; Evans reports a maximum 
e i 

ratio of about 1.6 just before noon. Both of these investigations were 
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carried out with about 2 Mw peak power at 440 Mc/s and cover a range of 
altitudes between 200 and 700 km. Greenhow, Sutcliffe and Watkins (l963), 
using 0.1 Mw peak power at 301 Mc/s, measured the diurnal variation of 
scattering cross-section (and hence of T g /T ) near the F2 peak. They 
found that = 3 during the day, but gradually decreases at night 

to a value of unity after midnight. They also find that by day, T /t. s 4 
near 200 km, but decreases upward. Bowles, Ochs and Green (1962), observ- 
ing at low latitudes, also find increases of T /t near sunrise 

e i 

Another proposed application of high-power radar equipment may be 
mentioned here. This is the artificial heating of electrons in the F 
region. Farley (1963) shows that small but probably detectable changes 
in the electron profile might be produced by the 50 Mc/s radar at Jicamarca. 
If a lower frequency, near f Q F2, were employed, it should be possible to 
produce increases of T^ of several hundred degrees and appreciable 
reductions of N^F2 with average powers as low as 100 kw. 
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Ill . PROCESSES IN THE IONOSPHERE 


1 . The Balance of Ionization 

We now turn to consider the physical processes which control the 
ionosphere. In this section, we present only the simplified theory 
against which the behavior of the actual ionospheric layers can be com- 
pared in Sec. IV. We deal here only with the large-scale structure of 
the ionospheric layers, and completely neglect small-scale irregularities 

which are considered in Sec. IV-6. 

The processes can be divided into two broad categories: those that 

result in production or destruction of ionization, and those that result 
in movement of ionization. The terms "photochemical" and "transport" 
serve as convenient, though not ideal, labels for these two categories. 
Later, we shall see that the relative importance of these categories varies 
with height; photochemical processes dominate the lower ionosphere 
(D and E regions), and we shall in fact suggest that the F2 layer 
represents a transition from a "photochemical" to a "transport" regime. 

It is usual to form an equation of "continuity" or "balance," whose 
terms represent the effects of the various processes which alter the 
electron density N. Within a cell of unit volume, we have 

[Rate of change of electron density n] = [Gain by production] 

- [Loss by destruction] - [Loss by transport] 

Continuity equations can also be written for the positive and negative 
ions, or indeed for any constituent whose concentration is subject to 

change. 

If the transport processes result in a net drift velocity v, then 
the loss due to transport is the divergence of the flux Nv. Using sym- 
bols q and L to represent production and loss, we have 


dN/dt = q - L(N) - div(Nv) 


(HI-l) 
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This is not unlike the heat conduction equation (l-io) of Sec. 1-5. 

Before we attempt to solve this equation, or even evaluate the terms, it 
is worth discussing its nature. As it stands, the equation contains 
derivatives with respect to space and time. But, except in a few special 
circumstances (as perhaps near sunrise) horizontal gradients of N and 
V are likely to be much smaller than vertical gradients. Horizontal 
variations generally involve scale distances of hundreds or thousands of 
kilometers, but vertical scales are a few tens of kilometers. So we can 
often retain just the vertical contribution to the transport term, 
3(Nw)/dh, where w is the vertical velocity. 

The principal production process for the creation of ion-electron 
pairs is generally accepted to be the absorption of solar UV and X ray 
radiation, at least in low and middle latitudes. Photons with energies 
greater than about 12 ev can ionize one or more of the major atmospheric 
constituents. This process not only produces the ionization but it also 
provides the heat input which is necessary to maintain the high tempera- 
tures which are found above the E region. The excess energy of the 
photon is transformed into the kinetic energy of the ion-electron pair 
and then the remainder of the ionization energy is transformed eventually 
to heat upon recombination. 

At high latitudes and during magnetic storms (and perhaps at other 
times), there also appears to be a significant production of ions and 
electrons from collisions between high-energy charged particles precipi- 
tated in the atmosphere and the neutral molecules. This "corpuscular 
ionization" might be produced by the leakage of charge particles stored 
in the Van Allen belts or perhaps by charged particles in the solar wind. 
The importance of corpuscular ionization is not yet well established and 
is a topic of considerable current interest. Finally, if negative ions 
should be formed in the lower ionosphere by an attachment process, the 
electrons can be released by photodetachment, which provides another 
mechanism for electron production. 

The important loss processes may be summarized as atomic ion and 
electron (radiative) recombination; molecular ion and electron (dissocia- 
tive) recombination; and, in the lower ionosphere, the attachment of an 
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electron to a neutral molecule. The production and loss processes which 
are believed to be appropriate to the ionosphere are discussed in some 
detail in Secs. III-3 and III-4. 

The transport term includes the effects of several processes. For 

many purposes, we may think of the plasma as a gas which represents a 

minor constituent of the atmosphere. Like the other constituents, it is 

acted on by gravity and by a force arising from any gradient in its own 

partial pressure. Unlike the other constituents, the charged particles 

* 

are also acted upon by electric and magnetic forces. The plasma will 
diffuse through the neutral air if the forces acting on it are not in 
equilibrium. Electrons and ions diffuse together because, if they did 
not, their separation would lead to a large electric field which would 
quickly bring them together. This so-called "ambipolar " or "plasma" 
diffusion is opposed by collisions between the charged particles (prin- 
cipally the ions) and the neutral gas. It proceeds rapidly in the F 
region, but slowly in the lower ionosphere, where collisions are frequent. 
We develop the equations relating to plasma diffusion in Secs. III-5 and 
I I I -6 , 

The plasma tends also to be set in motion by movements of the neutral 
air, which may be due to large-scale wind systems or to temperature 
changes. But, particularly in the F region, the plasma tends to move 
parallel to the geomagnetic field lines. It can also be moved across the 
field lines by large-scale electric fields. Under certain circumstances, 
especially in the E region, the ions and electrons may move with dif- 
ferent velocities so that an electric current flows. We find it convenient 
to postpone discussion of this topic until Sec. VI, which deals with the 
relationship between geomagnetism and the ionosphere. 

A particular form of plasma diffusion is the interchange of ionization 
between the ionosphere and the magnetosphere, where the positive ions are 
mainly protons. This requires a charge exchange process which has been 
discussed by Hanson and Ortenburger (l96l). Possibly a downwards diffu- 
sion of plasma from the exosphere could help to maintain the F2 layer 
at night, though it is not clear that this process would be sufficient. 


SEL-64-111 


52 



In the lower ionosphere, two simplifications can often be made. 

Transport processes are not very important, and if they are neglected 
entirely, a "photochemical Tt equation results, containing only the one 
derivative Sn/cH, Furthermore, the "time constant" associated with the 
loss term L(n) may be so short that /bt is much smaller than the 
other terms, so the "photochemical equilibrium" equation q = L(n) is 
adequate. This is generally the case in the D, E and FI layers by 
day, except for rapidly-varying phenomena such as eclipse effects. Trans- 
port can be then included as a small perturbation if required. 

2 . Chapman f s Theory 

A surprisingly good description of many of the features of the earth's 
ionosphere can be obtained from a theoretical study of the simplest model 
of the atmosphere and the incident flux. The theory was developed by 
Chapman ( 1931a, b) and has proven such a useful description that it is 
identified with his name. We will therefore consider an atmosphere com- 
posed of a single constituent with a monoenergetic UV flux incident 
upon it. A "flat earth" will be assumed at first but this can be general- 
ized to a spherically stratified atmosphere with little difficulty. The 
temperature profile will remain arbitrary, but we will later specialize 
our discussion to the case of a linear temperature variation with altitude. 
Chapman’s original work considered the isothermal case, which is readily 
obtained from the equations below. 

The intensity of the incident solar radiation beyond the atmosphere 

2 

will be denoted by I (photons/cm /sec). There is a certain probability 
that each molecule will be ionized in unit time when subjected to unit 
flux density, and this probability is usually expressed as a "cross-section" 
which will vary with wavelength and the molecular species. For our simple 
model, a single cross-section a is assigned. As the radiation passes 
through the atmosphere, the incremental reduction in intensity is 

dl = -I d n ds (lII-2) 
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in which ds is an element of path length in the direction of power 
flow, (Many derivations replace (cr n) with (Ap), in which A is the 
"mass absorption coefficient" and p is the gas density.) Since 
-ds = sec x dh, where \ is the zenith angle above a flat earth 

(dl/l) = cr n sec x dh (lII-3) 

When this equation is integrated, we obtain 


!( h )/l = 


-o* sec x n dh 
J h 


Using the result of Eq . (l-6), namely 


( 1 1 1 -4 ) 



dh 


n(h) H(h) 


0 - 6 ) 


we obtain 


l(h)/l = exp {-cr n H sec x) = (ill -5) 

which serves to define the "optical depth," T. The height at which 
T ~ 1, corresponding to a reduction of (l/e) in the incident radiation, 
provides a convenient measure of the depth of penetration of the radiation 
We will assume that the absorption of one photon produces one ion- 
electron pair, which is true for most of the UV spectrum. Then, the 
production rate q (pairs/cm /sec) is obtained from 

q(h) = dl/ds = cos X (dl/dh) (lII-6) 

q(h) = -I M cos X e -T (dx/dh) = -I OT cr e' T -jt (n H) (lII-7) 
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From Eq. (l-6) we see that d(n H)/dh = -n, so that 


q(h) 


I cr n( h ) 

CO V / 


-T 


e 


( IH-8) 


To locate the production peak, we set d(ln q)/dh = 0, giving 


1 dn dT d / y— \ 

— tt = "TT = O' sec Y — n H 
n dh dh ^ dh 


( I II -9 ) 


We see that q is greatest at the level where the downward increase 
of gas concentration just compensates for the rate of attenuation of the 
radiation. From the perfect gas law p = nkT we deduce, as in Eq. (l-7), 


1 dn 1 / dH\ 
n dh " " H T + dh/ 


(III-10) 


When Eqs. (ill -9) and (ill-io) are combined, we observe that 


n cr sec \ E = 


1 + 


dH 

dh 


T 


(in-11) 


at the peak' For positive scale height gradients, the peak of production 
lies just below the level of unit optical depth. 

This equation suggests two convenient levels at which our reference 
altitude, h = h Q , might be established. We could set (a) <r = 

1 + (dH/dh)^, that is at the altidue of maximum production for an over- 
head sun, or (b) n cr H =1, the altitude of unit optical depth for 
overhead sun. Although the choice is quite arbitrary, we will select 
the former definition. In the simple case in which H is independent 
of height, these levels naturally coincide; otherwise level (a) lies T 
scale heights below level (b), where = (dH/dh)^. 

We can now obtain the maximum production rate for overhead sun from 
Eq. ( I II -8 ) : 


q( h 0 ) 


I 

00 


-a+r,) 

cr n e 
o 


i (i + r ) 

oo v O 

H exp ( 1 + T ) “ Q o 
o o y 


(in-12) 
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We would like to put these equations in terms of the reduced height 

parameter, z. It is therefore convenient to specialize our equations to 

the case of a linear temperature and scale height variation. We may 

write H = H + r(h - h ) and then proceed as follows: 
o o 


z 



dh 

H + r(h - h ) 

o x o 



( III-13) 


obtaining 


H = H 

o 


+rz 

e 


(in-14) 


+Tz 

and there is a similar equation for T = T^e . The expression for the 
neutral concentration is 


n 


P 

kT 


p o 6 -z(i+r) 

+Tz = n o e 

k T e 
o 


(ill -15 ) 


When Eqs. (III-14) and (III-15) are inserted in (lII-8), we find 

q(z) = q Q exp [(l + r)(l - z - e" Z sec x) ] (III-16) 


(A rather different expression is obtained when h is identified with 

o 

the height of unit optical depth. ) 

To find the reduced height of the production maximum, we solve Eq. 

(Ill-ll) for z ; 
v ' m 

-z (i+r) +rz 

cr sec Vne He = 1 + T ' (III-17) 

v o o v / 


Therefore, 


z = In (sec y) 
m x 


( III -18 ) 
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which, of course, is zero for overhead sun. When Eq. (III-18) is inserted 
in Eq. (lll-16) we find the maximum production rate to be 

q = q (cos X) 1 ^ (III-19) 

m o' 

For the case of an isothermal layer (F = 0), Eq. (III-16) reduces 
to the simple Chapman formula 

q = q^ exp [l - z - e Z sec )(0 (ill- 20 ) 

This function has the Interesting property that as X varies, its shape 
is unchanged, its peak is shifted to the level = In (sec x) anc * its 

amplitude is scaled by the factor cos X- This can be shown be writing 
the above equation in the form 


q = (q^ cos x) ex P 


1 - (z - z m ) 


z -z 
m 
e 


(in-21) 


In Fig. 9, the ratio (q/q Q ) is plotted for several values of X • We 
see that above the peak, q depends very little on X- This is also 
illustrated in Fig. 10, which shows q(x) curves for several fixed values 
of z. At large positive values of z, q(x) is almost independent of X 
until X approaches 90 deg and sec X > e • As z decreases, the 
function q(x) becomes increasingly sensitive to X* * n the ionosphere, 
this means that the maximum of q, which occurs at noon when X is least, 
becomes increasingly sharp with decreasing z. Some of the properties of 
this function have been discussed by Cummack (l96l). 

In the actual ionosphere, the production formula is considerably more 
complicated. First, there may be several atmospheric gases, differently 
distributed, in which case the product (cr n) in Eq . (lII-2) must be 
replaced by a summation, Z o\n_^ . Some equations for a two-gas mixture 

were given by Rishbeth and Setty (l96l). Also, the ionizing radiation 
is not monochromatic, but consists of a range of wavelengths, and the 
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An approximation to the Chapman function can be obtained in the 
following manner. In Fig. 11, we are shown a point a distance h Q above 
the earth at which the sun's zenith angle is X 0 and for which we wil1 
evaluate the Chapman function. The minimum altitude of the ray is 
h = (R + h ) sin X " R , which we may call the "grazing height." 

gr ' O O 7 O 



FIG. 11. GEOMETRY USED IN THE 
CALCULATION OF Ch (X Q ). 


From the two preceding equations, we obtain 

J n sec X dh 
Ch (Xq) = H n(h Q ) 


(HI-24) 


although this is not in the same form as originally expressed (Chapman, 
1953). Expanding these terms, 


00 

f n(h ) sec X ex P ' 
J h e 

S' 

h - h 
£ 

1 H 

dh ± 

n 

r ° 

j n(h ) sec X ex P 

h g 

_i_ — 

h - h 
H | 

l ‘ 

dh 

H n(h ) exp 
g 

h - h 
g ° 

H 



(HI-25) 


SEL-64-111 


60 



In this equation, the upper sign is used for X Q > 90 deg, while the 
lower sign is applied when X Q < 90 deg. 

If we now measure the reduced height from the grazing level, so that 
T = (h - h^)/H, the last equation becomes 


Ch (x) = e ° 

CO 

f -T J , 

1 / e sec X dT ± 

T 

r O 

f ~T 

/ e sec x dT 


J 

[ 0 

0 


(III-26) 


Some additional work will show that these integrals may be approxi- 
mated by 

i e ' T sec X dT ~ (III-27) 

which is sufficiently accurate as long as H « R^, so that exp ( -T ) -► 0 
well before the height variable approaches the value R^ . Also, it should 
not be applied for small X> where the optical depth is proportional to 
sec X anyway. However, it is quite a good approximation near X “ 90 deg 
where the Chapman function is required, and is quite useful in almost all 
cases of ionospheric interest. The integrals may then be evaluated in 
terms of the error function, giving 


HR 

Ch (X 0 ) ~ /- 2 h exp (t q ) (1 ± erf (III-28) 

Since the error function is tabulated in many places, Eq. (ill -28) 

may be more conveniently applied than Ch (X Q )> anc * ma Y a lso be more 

readily evaluated in a computer. Ch (x ) has been tabulated by Wilkes 

(l954) as a function 'of Y and the ratio of the ionosphere radius to 
„ v o 

scale height. Ch (x ) departs markedly from sec X for X ~ 80 de S 
and for large H, as shown in Fig. 12. The approximation in Eq . (lII-28) 
does not depart significantly from the true values of Ch (X Q ) shown 
below. 
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FIG. 12. SEC ^ AND Ch (*<,) VS X©- 
PARAMETRIC IN (r/h). R and H are 
the geocentric radius and scale height 
of the atmospheric layer. 

When considering problems at large zenith angles, it is necessary to 

replace sec X with Ch (X Q ) wherever it occurs in the preceding equa- 
tions. From Fig. 12, it is clear that the optical depth does not become 

infinite at y = 90 deg and that significant production may occur for 

zenith angles as large as 95 deg or 100 deg. This is particularly true 
at the greater altitudes, because the ray paths from the sun do not 
penetrate deeply into the atmosphere even for \ ~ 100 deg. 

Having found the production function q(z), it is a simple matter to 
find the photochemical equilibrium electron density profile when the 

transport term in Eq. (ill-l) is neglected. When electron loss is pro- 

2 2 
portional to N , the continuity equation becomes cSl/ht = 0 = q - CtN , 

from which 
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N = 


Vq/a 


(hi - 29 ) 


and therefore the peak electron density and the critical frequency 

of a Chapman layer are connected by 


f 

o 


(cos X) 1 / 4 . 


(in-30) 


If there is a vertical temperature gradient, we may expect the exponent 
of cos X to t> e different. We can see from Eq . (HI-19) that if 
T = dH/dh is independent of height, 


f <x 

o 


(cos x) 


0.25(l+r) 


(in-31 ) 


3. Production and Loss 

The production of ionization is principally complicated by the wide 
spectrum of the radiation generated in the chromosphere and the corona of 
the sun. Wavelengths from a few Angstroms to almost 2000 % are all 
important in the ionosphere. The spectral lines characteristic of hydrogen, 
helium, carbon, silicon and many other elements are all found in the 
radiation spectrum. 

We should also distinguish between an absorption cross-section which 
expresses the probability that a given molecule will absorb a photon and 
an ionization cross-section which is the probability that an ion-electron 
pair will be produced. In the previous section on Chapman layers, they 
were assumed identical. More precisely, cr should be interpreted as the 
absorption cross-section and then each equation for production multiplied 
by the ratio of the ionization to the absorption cross-section. Since 
these cross-sections are dependent upon wavelength and molecular species, 
we may expect the incident radiation to be absorbed over a wide range of 
altitudes as implied by Eq . (ill-ll). 

The best current values for flux intensity vs wavelength come from 
the rocket measurements of Hinteregger and Watanabe (1962) and the mea- 
surements made in the Orbiting Solar Observatory (Lindsay, 1963). Fig. 13 

[ ■ - SEL-^4-111 


shows the EUV flux vs wavelength and is derived principally from the 
references just mentioned, but it includes some results of other experi- 
menters. From Eq. (III-12) we can calculate the maximum production rate 
at any wavelength, but the absorption cross-section is needed to compute 



FIG, 13. EXTREME -ULTRAVIOLET (EUV) FLUX VS 
WAVELENGTH, BASED ON THE DATA OF WATANABE 
AND HINTEREGGER (1962). Several of the 
more prominent line emissions are shown, 
along with the distributed radiation. The 
5000 °K black body emission is shown as a 
broken line, and it may be seen that the 
solar radiation exceeds this value below 
about 1100 %. 
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the altitude at which it occurs. Based on Hinteregger 1 s cross-sections 
and the model atmosphere used in Sec. I, Fig. 14 shows, as a function of 
wavelength, the height of unit optical depth, at which vertically inci- 
dent radiation is attenuated to a fraction l/e of its intensity above 
the atmosphere. We can now compute the production profiles as well, and 
these are shown in Fig. 15 for the case of overhead sun. They refer to 
a fairly low level of solar activity, current at the time of their rocket 
experiment in 1961. We see in Fig. 15a that q[o + ] is greatest at about 
170 km. 



FIG. 14. THE ALTITUDE OF UNIT OPTICAL DEPTH VS 
WAVELENGTH (after Friedman, I960). 


Figure 15b shows that q [ N^ ] also peaks at this level, though as 
mentioned later it is uncertain how much the ionization of N^ contributes 
to the observed electron density. According to Fig. 15c, a major part of 
the E region ionization arises from the wavelength band 911 - 1027 2, 
which cannot ionize 0 or N , and which includes the very strong solar 
Lyman p line at 1026 X. This greatly exceeds the contribution made by 
X rays in the range 10 - 170 X, although Norton, Van Zandt and Denison 
(1963) consider the latter to be more important in the E region. 
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3)282 PRODUCTION RATE {elec /cmvsec) 


FIG. 15d. ILLUSTRATING THE ZENITH ANGLE DEPENDENCE 
OF THE PRODUCTION yffiFFOR' 0+, Ng and o£. 

The same models were u sed as in Figs. 15a, b,c. 

Three curves are drawn for each zenith angle 
X = 0, 90, and 98 deg. The curve with the small- 
est production rate in each set is q[o + ]. The 
difference between this curve and the middle 
curve in each set is q^Nj], The difference 
between the middle curve and the total production 
rate is qfOrj], as is shown for \ = 0 deg. 

Ionization can be produced not only by solar photon radiation, but 
also by energetic particles which can enter the atmosphere most easily 
at high magnetic latitudes. The enhancements of ionization observed at 
high latitudes during magnetic disturbance are certainly due to fast 
electrons or protons. The mechanisms by which these particles acquire 
their energy are part of the general problem of high-latitude and mag- 
Tibtic storm phenomena. We note, however, that the depth to which parti- 
cles penetrate depends on their energy. Some calculations of penetration 
depths have been given by Bailey (1959) and Rees (1963). It has been 
suggested that corpuscular ionization is responsible for much of the ioni 
zation and heating in the F2 layer (Antonova and Ivanov -Kholodnyy, 1961 
Harris and Priester, 1962a). The only radiation likely to be absorbed at 
these heights would be soft electrons of not more than a few hundred ev . 
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The ionization process would load to emission of the 6300 8 airglow line 
of atomic oxygen, and this can be used to set an upper limit to its rate. 

The luminosity vs height profile for a monoenergetic source has the 
same shape as the ionization profile; it is only necessary to alter its 
magnitude. For the 3914 8 radiation from excited N*, it is estimated 
that about 1 photon is emitted for every 50 ion-electron pairs produced, 
and each pair requires about 35 ev for ionization. If the luminosity 
profile is measured with a photometer, the production profile may be 
calculated from which the energy and intensity of the incident flux may 
be estimated (Rees, 1963). 

The kinds of photochemical reaction which are thought to control the 
ionospheric electron density are considered next. The idealized chemical 
equations are given although the participating ions and neutral particles 
will not be identified in most cases until the experimental data are dis- 
cussed in the next section. The mechanisms for loss of ionization have 
been discussed in two important articles by Bates and Massey (1947); the 
pertinent reactions are listed below, 

(a) Ion-ion recombination (coefficient ) 

X + + Y~ -» X + Y 

(b) Electron-ion recombination (coefficient ct ) 

Three-body: X + e+ M-*X + M 

Radiative: X + + e -> X* X + hv 

Dissociative: XY + + e -* X* + Y* 

The asterisks indicate that the atoms may be left in excited states, and 

subsequently lose this energy by radiation or during collisions with 

other particles. The symbol M denotes a neutral particle which exchanges 

energy and momentum but does not take part in the chemical reaction. 

Three-body processes can occur in the lower D region, but are so rare 

at greater heights as to be quite unimportant. Apart from this means, 

recombination of electrons and atomic ions can take place only by the 

-12 3 

very slow radiative process (coefficient ~ 10 cm /s) or possible by 
other equally slow reactions. Only in the uppermost F region is 
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Hiative recomblnatioji likely to be the fastest loss process, and »e 
"all see that at such heights transport processes (especially diffusion) 
are so completely dominant that loss coefficients are irrelevant. Else- 
where in the F and E regions, the dissociative recombination process, 


_____ -7 -8 

with coefficients of order 10 or 10 


(cm'Vs) , are more important. 


large p rgjjdfftio n of the ions are originally atomic, dissocia- 
wiFfej- preceded^y reactions involving formation of 

molecular ions, namely: 


(c) Ion-atom interchange (rate coefficient y) 

pllgfeo - x + + YZ - XY + + z 

While it is generally accepted that ion-atom interchange followed by 
dissociative recombination is the principal loss process in the E and 
F regions, there is considerable controversy as to precisely which 
reactions are important. Although charge-exchange reactions can also 
lead to formation of molecular ions, the ion-atom interchange process is 
thought to be more rapid (Bates, 1955). 

(d) Attachment (coefficient a) and collisional detachment 
(coefficient d^) 

X + e + M ft X +M 


(e) Associative detachment 

x” + Y -* XY + e 


(f) Photodetachment (coefficient d p ) 

x" + hv -» X + e 


A summary of the conclusions of this section is provided in Table I with 
additional comments regarding the importance of the various processes. 


4 The D, E and FI Photochemical Regime 

In the lower ionosphere, transport may be neglected and only photo- 
chemical terms appear in the continuity equations for the ion concentra 
tions N , N and the electron concentration N e - For the present, we 
add the + subscript "e" for the sake of clarity. These equations are 


therefore written as 
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processes. Hence, by day, A decreases upward, and is probably small 

above 90 km. At night, however, d p « 0 and A depend on the ratio 

a/d Above 100 km, ionic recombination would limit the equilibrium 
' c ' , 

nighttime value of N_, and at greater heights where 1/a exceeds a 

few hours, equilibrium conditions are scarcely applicable. 

When, as in the daytime D region, attachment and detachment are 

so rapid that A is constant, we can substitute Eq. (III-35) into Eq. 

( III -32) and obtain 

2 

(1 + A) dN e /dt = q - (1 + A)(a e + ta.)N e 

= q - CLN 2 ( HI “37 ) 

E e 

The so-called If ef f ective ,f coefficient as defined in these equa- 
tions, rapidly approaches with increasing altitude because of the 

decrease of A. 

So far , we have described the chemical processes in rather general 
terms. In reality, discussion of the D region is hampered by the fact 
that the ions present have not been identified. Since nitrogen does not 
form stable negative ions, it has generally been assumed that 0 or 
0^ is the most important. There is, however, some controversy about 
the electron affinity of atomic oxygen, which is important because it 
determines whether visible light or UV is required for photodetachment, 
and thus influences our interpretation of D region phenomena, especi- 
ally at sunrise. Other suggestions, such as 0 g and NOg, as to the 
dominant negative ion have been advanced (e.g., Reid, 196l). Neither 
have the positive ions been identified, but N0 + (which, as we shall see, 
is the dominant ion in the E layer) may be the principal ion. It can 
be formed either by photoionization of NO or by various chemical reactions 

Presumably 0* and N+ are also present. Further consideration of D and 

E region photochemical processes may be found in Nicolet and Aikin 
(i960), Crain (l96l), Poppoff and Whitten (1962, 1964), Pierce (1963), 
and Webber (1962). 
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Let us now examine the situation at greater heights, in the E and 
FI regions, where the reactions involving negative ions are unimportant. 
Under these circumstances, the only important recombination process is 
dissociative recombination of electrons and molecular positive ions. 

Since a large proportion of the neutral atmosphere is atomic, especially 
in the F region, atomic ions are produced by photoionization, but these 
do not recombine with electrons directly, except by the very slow radia- 
tive process. Instead, they undergo an ion-atom interchange reaction and 
the molecular ions thus formed combine with electrons. We shall consider 
the electron distributions resulting from this situation in a general 
way, before discussing the actual reactions which occur. 

In writing the continuity equations for the concentrations of elec- 
trons N, atomic ions N and molecular ions N__, we assume that the 

-a M 

ion-atom interchange reaction involves molecules "m, " but that we can 
neglect the direct production of molecular ions by photoionization. This 
could be a realistic situation, since it is possible that the dissociative 
recombination rate of N 2 is quite large (perhaps 10 7 or more) in which 
case the ion would recombine so rapidly that it would not contribute to 
the "observable ionization," Proceeding on this assumption for the moment, 
we write 


dN/dt = q - aNN M 


dN A /dt = q - 7n(m)N A (III-38) 

dN M / dt = 7n(m)N A - aNN^ 


Charge neutrality again requires N = N. + N . If we write 6 = m(m) 
and assume equilibrium, we can eliminate the ionic concentrations N 

A 

and V and obtain a quadratic equation in N (Hirsh, 1959): 

aPN 2 - aqN - Pq = 0 (ill- 39) 
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of which the positive root is 


This reduces to 


N = (q/2p) [l + (1 + 4 P 2 /aq) l/2 ] 


(ill -40 ) 


N = N a = / q/a if 4p 2 » aq 


N = Np = q/p 


if 4p « aq 


(ill -41 ) 


or N 2 = N + NN„ 

a P 


in general 


Alternatively, the quadratic can be written as 


1 _ JL_ _1_ 

q ~ PN + 2 

aN 


( III -42) 


which reduces to 


q = aN 2 if p » aN 

q = PN if P « aN 


( III -43) 


The conditions "»" and "«" occurring in Eq. (III-43) are, of course, 
roughly equivalent to those in Eq. (lll-4l). 

These equations demonstrate how the two-stage loss process proposed 
by Bates and Massey (1947) gives rise to a transition between a "quadratic" 

n 

loss law (aN ) and a "linear" loss law (PN). Returning to the balance 

equations (III-38), we see that N./N^ = aN/p at equilibrium, so that 

when the rate of electron loss is determined by the dissociative recom- 

2 

bination reaction, the ions are mainly molecular and the "ON law" 
applies. The coefficient a may depend on temperature, but is otherwise 
not height dependent, unlike P which varies as the molecular concentration 
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n(m) and therefore decreases rapidly upward. So at greater heights, we 
expect the condition P « aN to apply, in which case the ion-atom inter- 
change reaction controls the rate of loss and the ions are mostly atomic. 

As we shall see in a moment, rocket data on the ionic composition con- 

2 

firm that this transition between Tf CfcN " and M (3N M loss laws occurs in 
the lower" F region, at about 160 - 200 km. This happens to coincide with 
the level at which the F region production rate q is greatest, and 
Ratcliff e (1956) showed that this might account for the splitting of the 
F2 layer into FI and F2 components. To demonstrate this, we investi- 
gate the shape of the equilibrium electron distribution by evaluating the 
formula (ill -40 ) using the Chapman production function q(z,x) for an 
isothermal layer, Eq. (lll-2l), and a height -independent recombination 
coefficient a. If we introduce a new quantity k = IF/h^, the ratio of 

the scale heights of the ionizable gas and the linear loss coefficient, 

— kz 

we can write p e 

Since P = 7n(m), k depends on the molecular mass of the gas partici- 
pating in the ion-atom interchange reaction, which may be 0^ or N^. Since 
the principal ionizable gas is 0, we have k = 32/l6 = 2 for 0^ and 
k = 28/16 =1.75 for provided these gases are diffusively separated, 
which is found to be the case above the E region. If the coefficients 
y were very temperature dependent, the value of k might be modified, 
but probably is not very different from the above-mentioned values in the 
F region. 

We found in the previous section that the shape of the Chapman produc- 
tion function does not change as the solar zenith angle X varies, but 
that the level of the peak varies according to the relation z = In sec X- 
The form of Eq . (lII-40) suggests that a single parameter, (3 /aq, deter- 
mines the shape of the electron density distribution N(z). Let r denote 
the value of P/y^Qq at the level of peak production. Provided r (and 
k) are kept constant, changes of q and p can alter the magnitude of 
the function n(z) and displace it with respect to the z-axis, but do 
not affect its shape. We note that from Eq. (lll-4l), 


r = p/y aq = N /n 


a' P 


(ill -44 ) 
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The graph is drawn 


and in Fig. 16 plot n(z) for r = 1, r = 2, r = 3. 
in such a way that the dotted curve, which represents N 
z = 0, is the same in each case, but different values of 
so as to give three dashed curves N^ , and three solid curves for N. We 
see that if r = 1, the n(z) profile is almost smooth, but that a 
"ledge" appears for r = 2 and is more prominent for r = 3. On an 
h'(f) curve, such a ledge would produce a prominent cusp, as shown by 
Hirsh (1959), and we may identify this with the splitting of the F layer 
into FI and F2 components. 


= yJ~^/cx for 
P are chosen 



FIG. 16. ELECTRON DENSITY PROFILES FOR THE "TRANSITION 
REGION," ASSUMING A CHAPMAN PRODUCTION FUNCTION q(z) 
WITH PEAK AT z = 0 AND A RECOMBINATION COEFFICIENT 
a, INDEPENDENT OF HEIGHT. The linear loss coefficient 
is (3 0 exp ( -1 . 75z ) , and three values of P Q are 
used, such that r = $ Q /yJ aq Q takes the values 1, 2, 
3. The N(z) profiles correspond to: 

....... N a = y/lja 


Np = q/p, 


for three values of B 


N computed from full equation (III-43) 
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The relevance of this analysis of the F region lies in the fact 
that, since r is evaluated at the level of peak q, it varies with 
solar zenith angle Remembering that the level of peak production 

is given by z = In (sec )(), that the peak value of q is (q cos x) 
and that p>(z) = p^e , we have 

r (x) = p( z )/v/c«r = p 0 ( c ° s x) k /y c ° s x 

= r(o) (cos x ) k ' l/2 (HI-45) 

Thus, if r(o) is constant, as would be expected on a first-order theory 
of the atmosphere, we see that as increases, r()() decreases and the 

"splitting" of the layer becomes less pronounced, or even disappears. 

This is consistent with the observations that, at mid-latitude stations, 
the FI layer is most prominent around noon, and is more commonly observed 
on summer days than in winter. Since r(o) oc l/^q” and varies with 

the solar cycle, we can also account for the observation that the FI 
layer is more prominent at sunspot minimum than at sunspot maximum. More- 
over, the FI layer sometimes appears during a solar eclipse at times 
when it would not normally be seen, and this may arise from the reduction 

of q and consequent increase of r. 
o 

The above discussion does not depend critically on the identity of 
the reactions, but we should now be more specific about the ionic con- 
stituents. Examples of the experimental results obtained from rocket - 
borne ion mass spectrometers are shown in Fig. 17 (Johnson, Meadows, and 
Holmes, 1958). Below 150 km molecular ions (NO + and 0*) are predominant, 
but atomic ions (0 + ) dominate above 200 km. This is broadly consistent 
with the "transition" theory outlined above. Other ions, such as 
and N + are present as very minor constituents. We have previously shown 
evidence that the major neutral constituents in the thermosphere are 0 
and N , with some 0 . Thus the ions produced by photoionization are 0 , 
and since 0^ is a minor constituent in most of the thermosphere, 

we expect most of the 0* ions to be produced by "transfer" reactions 
(i.e., ion-atom interchange or charge exchange), and this is certainly 
true of the N0^~ ions. 
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FIG, 17. RELATIVE ION CONCENTRATION VS ALTITUDE FOR THREE ROCKET 
FLIGHTS (after Johnson et al , 1958). 

In Table II we give a list of reactions which is not exhaustive, but 
includes those reactions which seem likely to be of greatest importance. 
Reactions which have been omitted include those which are endothermic and 
some which are probably unimportant because they involve only "trace" 
constituents (such as atomic nitrogen) or for other reasons. Orders of 
magnitude of the coefficients J and a have been given in Table 1, 

The formation and decay scheme sketched above is the simplest which 
could, in principle, account for the observed ion compositions. The 
equations relevant to this scheme have been given by Yonezawa and Taka- 
hashi (i960). It should be noted that the production of atomic nitrogen 
by these reactions does not necessarily lead to an appreciable concen- 
tration of this gas in the atmosphere; this will depend on the removal 
mechanisms which we do not discuss here. 

Since N* ions are so scarce in the E and F regions, this scheme 
would imply that they are destroyed so rapidly by the recombination (R3) 
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TABLE II. PHOTOCHEMICAL REACTIONS 


Photoionization 


0 + hv -> 0 + e 


+ hv + e 

+ h v -> Og + e 


Transfer or Interchange 
0 + + 0 2 - 0 ^ +0 
0 + + N 2 -> N0 + + N 

N* + 0 -» N0 + + N 

N 2 + 0 -* 0 + Ng 

Ng + ° 2 -* N0 + + NO 


Dissociative Recombination 


°2 + 

e 

-*• 0 * 

+ 

0 ** 

N0 + + 

e 

- N* 

+ 

0 * 

N 2 + 

e 

- N* 

+ 

N** 


(rate q) 

(Qi) 

(Q2) 

(Q3) 

(rate coefficient 7 ) 
(Tl) 
(T2 ) 
(T3) 
(T4) 
(T5 ) 

(rate coefficient a) 

(Rl) 

(R2) 

(R3) 


Simple Scheme for Positive-Ion Formation and Decay 


o-^o + 


Q2 


Rl 

< 1 > 2 °’° 



v’X + R2 

NO — » N , 0 


N 2 R3 > N , N 
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that they do not contribute to the observed ionization, though this does 
contribute to the input of heat. However, it now seems that this scheme 
is too simple, and that the coefficient of (R3) is not large enough to 
account for the scarcity of N* ions, so that additional reactions [such 
as (T3), (T4), (T5 ) or others] must be operative. There is considerable 

controversy even as to which of these is possible, let alone which is 
most important. Further discussion of possible reactions and data on the 
results of laboratory determinations of reaction rates are contained in 
the comprehensive reviews by Nawrocki and Papa (l96l) and Whitten and 
Poppoff (1964) and in the paper by Nicolet and Swider (1963). 

Although measurements have been made of the dissociative recombination 

coefficients for 0 +, N0 + and N+, and of the rate coefficients for the 

transfer reactions (Tl) and (T2), none of these parameters can be 

regarded as well known. In some cases, different results are at variance 

by an order of magnitude. Thus, although it is easy to write down the 

continuity equations for the various positive ion concentrations and the 

electron concentration, and to devise computer methods for their solution, 

it is not so easy to reach definite conclusions about which particular 

reactions are most important, even though the ionic composition is known 

as a function of height. However, it is possible to make some rather 

general statements, such as the necessity for some removal process for 

N + ions besides dissociative recombination; and to find ranges of values 
2 

within which the various coefficients night lie. 

There is also a long-standing difficulty concerning the rates of 
reactions (Tl) and (T2), which almost certainly make the major con- 
tribution to the loss coefficient (3 in the F2 layer (e.g.. Bates and 
Nicolet, I960). We have 

P(h) = 7 T1 n[0 2 ] + 7 T2 n[n 2 ] (III-46) 

Most available data on all the quantities on the right-hand side of this 
equation give values of p which seem to be too large for the F2 layer, 
although some recent values quoted by Langstroth and Hasted (1962) for 

7 and 7 may be small enough to remove the discrepancy. 

Tl T2 
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The products of some reactions, notably the dissociative recombina- 
tions ( R1 ) - (R3) in Table II, are formed in an excited state and may 
subsequently emit radiation. The energy-level diagrams for the atoms are 
described by Chamberlain (1959). Sufficient energy is liberated in 
reaction (Rl) to raise at least one atom to a state from which it 
cascades to the ground level, with successive emission photons of 5577 X 
green line and the 6300/6364 X red lines. The intermediate state has 
such a long lifetime (~1 min ), however , that at E region heights the 
atom is more likely to suffer collisional deactivation than to emit the 
red line. However, reaction (R2) liberates only enough energy to 
excite the red lines, 6300/6364 X. Nighttime rocket experiments have 
shown that nearly all of the red emission originates above 150 km, and 
nearly all of the green line originates below this height (Heppner and 
Meredith, 1958; Tousey, 1958). The absence of green emission from the 
F2 layer can be used as evidence that reactions (T2) + (R2) are more 
important in the F2 layer than (Tl) + (Rl). Since the red line is 
emitted from the F region, one would expect its intensity to be related 
to the integrated rate of loss of ionization, J PN dh. Correlations 
between its intensity at night and the critical frequency ^ ave 

been found at low latitudes (Barbier, 196l) but not at middle latitudes 
(Duncan, I960). Since it is now possible to observe the 6300 8 line 
even by day (Noxon and Goody, 1962), this technique should provide a 
valuable check on loss coefficients in the F region. 

5 . Plasma Diffusion 

The foregoing description of F region photochemical processes fits 
the hypothesis of Bradbury (1938), concerning the formation of the F2 
layer at about 300 km. No mechanism seemed to be capable of causing a 
peak of production at such a height, and so Bradbury supposed the produc- 
tion peak to lie at a lower height (now known to be near the FI layer) 
and attributed the upward increase of electron density to a rapid upward 
decrease of loss coefficient. 

This hypothesis does not of itself account for the location of the 
F2 peak. If the linear loss coefficient (3 varies as e (the ratio 
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k = IT/h^, as defined previously), then the electron density well above 
the production peak is approximately q/{3> « exp [ (k-1 )z ] , If k > 1, 
as seems nearly certain, then q/f3 increases indefinitely upward, and 
we must find some explanation for the existence of a peak of N. The 
possibilities include: 

(i) Failure of the Chapman formula for q at heights where the 
ionized/neutral concentration ratio (N/n) is not small, 

(ii) Existence of an additional loss process, such as radiative 
recombination, which might dominate at great heights. 

(iii) Lack of equilibrium, such that N never approaches the 
limiting value q/|3. 

(iv) Action of some transport process (such as diffusion) to 
limit N at great heights. 

Of these possibilities, (iv) and perhaps also (iii) are most likely 
to control the F2 layer; in particular, the process of plasma or 
ambipolar diffusion is thought to limit the upward increase of electron 
density and lead to the formation of the F2 peak. Diffusion arises 
from the tendency of the electron-ion gas or plasma to assume a hydro- 
static distribution under gravity. We recall that the neutral atmospheric 
gases are in hydrostatic equilibrium above 100 km, but the great chemical 
activity of the ionization causes its distribution to be determined by 
photochemical processes right up to the F2 peak. 

The importance of diffusion was suggested by Hulburt (1928) and the 
mathematical expressions were derived by Ferraro (1945). Early solutions 
of the continuity equation, with diffusion included (Mariani, 1956), 
treated it as a small perturbation only. But Yonezawa (1956) discussed 
in detail the problem of the formation of the F2 peak and showed that 
diffusion could provide an explanation. 

To formulate the diffusion equation for the ionization, we balance 
the partial ion pressure gradient against the total force per unit volume 
acting on the ions. These are gravity, the electrostatic force due to 
the macroscopic electric field E which we assume to be present, and 
frictional forces due to collisions with ions and electrons. Considering 
only vertical forces and velocities, the equations for ions and electrons 
are: 
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( II I -47 ) 


S(N.kT.)/Sh = -N.m.g + N. eE - N. m. v. (w - w ) 
x i i 7/ 11 i l i in i n 

S(N kT )/Sh - -N m g - N eE - N m v (w - w ) (ill -48) 

v ee ee e eeen e n 


in which we symbolize concentration N, Boltzmann's constant k, tem- 
perature T, gravitational acceleration g, mass m, electronic charge 
-e, collision frequency v (with neutral particles) and vertical drift 
velocity w. Suffixes i, e, n refer to ions, electrons, and neutrals. 
The frictional terms can be derived from kinetic theory. We simplify 
the equation by the following assumptions: itu » = N^, w_^ = - 

w , w =0 (neutral air at rest), m. v. » m v (collisions with 
neutral particles important for ions but not for electrons). On adding 
the equations, the terms in E vanish (and, because of Newton's Third 
Law, so do the ion-electron collision terms which we omitted from the 
equations). This means that the electrons interact with the ions only 
via the electric field. Rearranging to solve for the drift velocity of 
the plasma 



1 


m. v 
i in 


1 b 

N <5h 


[Nk(T. 


+ T ) ] + nug 


( III -49 ) 


In the F region, the electron, ion and neutral temperatures may all be 
different. Except at great heights, it is probable that T = T (the 
suffix n has been dropped). Then we introduce the neutral scale 
height H = kT/m g and set u == m./ 2m and T = T /t. = T /t. The 
ion-neutral diffusion coefficient is defined as D = kT/nL so that 

the^ equation becomes 


-w_ 


/, I“l 1 St 2 \i bnc/bh 

= D m (1 + t) [n 3h + T 3h + (1 + T)H + 1 + T. 


( II I -50 ) 


and 


if T e = T^ = T at all heights (i.e., T = l) 


-w_ = 


nT , r i c^n i u] 

2D - - [n 5h + T 3h + 5 ] 


D in 


(ill -51 ) 


- 83 - 


SEL-64-111 


Since D - 2D. = 

m 

gas concentration and 


2kT/m^v. , and v is proportional to both the 
/T. we may write, to a first approximation, 



n 


( III -52 ) 


Ferraro (1945) uses a more complicated temperature dependence, so that 
b is a slowly-varying function of T; this is further discussed in a 
later note (Ferraro, 1957). Values of b can be calculated from theory 
(Chapman and Cowling, 1952) but the factor of 2, discussed above, was 
omitted by Ferraro. This was pointed out by Johnson and Hulburt (1950), 
who treated the whole diffusion problem from a different standpoint. 

Ferraro assumed a mean molecular mass of 25 a.m.u. for the plasma, whereas 
it is now believed that in the F2 layer the neutral gas is mainly atomic 
oxygen and the ions mostly 0 + . Thus, numerical values of b must be 
treated with caution, especially since Dalgarno (1958) finds that, for 
ions diffusing through their parent gas, b is reduced by charge-exchange 
between ions and neutrals, possibly by a factor of four. There are also 
many variations of notation in the literature; for instance, some authors 
write D = b/n and so include the dependence on T within the quantity b. 

The contribution of diffusion to the continuity equation can be written 
as 


-d(Nw D )/dh = DJ9N (III-53) 

where J9 is a differential operator. The evaluation of J9 involves 
fairly complicated algebra, but expressions have been given by Shimazaki 
(1957) for the case in which the scale height gradient dH/dh is inde- 
pendent of height. We shall restrict our derivation to the isothermal 
case, dH/dh =0. We also take p = 1/2 (corresponding to m n = m ^) in 
Eq, (lll-5l) and express J9 in terms of reduced height z by taking 
H • b/bh s b/bz . For this simple model, the neutral particle concentra- 
tion n ce e Z so that Dec e +Z . Using these relations, Eq. (lll-5l) 
reduces to 
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ID 1 3R 


’ W D = H N ^ + 2 


and by Eq. (ill -53) the diffusion term in the continuity equation takes 
the form derived by Ferraro (1945) 


n ,0 XT D ON 3 dN N 


(III-55) 


The effect of gravity and the height dependence of D are responsible 
for the appearance of terms in 3 n/3z and N, in addition to the second 
derivative <3 ^n/c)z^ which is characteristic of diffusion formulas. 

Since D increases exponentially upward, whereas the other coeffi- 
cients q and (3 in the continuity equation decrease upward, at some 
level diffusion dominates so completely that the continuity equation 
reduces to ]9 n = 0. The solution of this equation is 


XT . -1/2 z A -z 

N = A^e + e 


(in-56) 


The first term corresponds to diffusive equilibrium, with jv = 0, such 
that the ionization assumes a scale height twice that of the ionizable 
gas. For the second term, w^ ^ 0 ; and this represents a boundary con- 
dition of a finite flux of ionization at z -> + °o f which is upward if the 
coefficient > 0. If ionization is gained or lost by diffusion at the 
top of the ionosphere (i.e., to or from the magnetosphere), the N(z) 
distribution should contain a component of this type. 

Finally, we note that in the F region the ionization can diffuse 
only along the lines of geomagnetic force. The usual way of treating 
this problem is to use the components of w^ and (S/Sh) parallel to 
the field. This introduces a factor of sin^ I into the coefficient 
D, where I is the magnetic dip angle (Ferraro, 1945). However, this 
approach is incorrect if the neutral air is in motion, as may be the 
case (Dougherty, 1961 ). Near the magnetic equator the "sin I" factor 
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fails and a more complicated form of the operator J9 is required which 
includes horizontal as well as vertical diffusion. This has been derived 
by Kendall (1962) and Lyon (1963). 

The equation for the diffusion of ions (lII-50) serves as a convenient 
starting-point for discussing the equilibrium distribution of charged 
particles in the topside F region. At these heights, photochemical 
processes can be neglected, and the collision frequency v in is so small 
that the continuity equation reduces to the partial pressure, gravita- 
tional and electric field terms of Eq. (III-50). The situation is interest 
ing because several species of positive ion may be present. We suppose all 

ions to be singly charged and to possess the same temperature T , though 

T. .th 

this need not be the same as the electron temperature T^. For the j 
species of ion, and for the electrons, the equations are 

d(N kT. )/dh = -N .m .g + N .eE (III-57) 

v J i J J J 

d(N kT )/dh = -N m g - N eE (III-58) 

v e e e e e 


Since ZN = N , the equations can be added to eliminate E. We then 

j e 

denote the mean positive ion mass by m , so that m N = ZN m and set 

“I" +0 J J 

T /t. = T obtaining 
e i 



N dh 


e 


gm + 

kT . ( 1 + T ) 
1 ' 


(ill “59 ) 


when the small gravitational term (^ e m e &) neglected. The electric 

field E may now be found from Eq . (ill -58) and inserted in the ion 
equation to give 



m . 
J 


Tm 

4 - 

1 + T 




JL. 

kT. 


(III-60) 


We have neglected temperature gradients for simplicity, although they 
could be readily included. A further simplification is to assume thermal 
equilibrium, T = 1, which reduces the last two "scale height" equations to 
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H e = 2kT/(m + g) 


and 


(ill -61 ) 


H. = kT/(n,.g 


I \ g ) 


We see that if only one ionic species is present, its ’’effective scale 

height” is twice that of the neutral species of atomic mass m (as 

before). However, a light ion for which m. < — m actually has a 

J 2 + 

negative scale height, so that N^ increases upward, as discussed by 
Mange (i960). The Eqs. (III-59) and (lII-60) can be solved by a method 
due to Hanson (1962). It is convenient to use the "geopotential height’ 
defined below Eq. (l-9) in Sec. 1-3, in order to take account of the 
vertical variation of g. A specimen equilibrium distribution computed 

4. -i_ 4: : i ■ 

for a mixture of 0 , He and H ions, is shown in Fig. 18. 



FIG. 18. IDEALIZED DISTRIBUTIONS OF ELECTRONS (e) 
AND OF 0+, He + and H+ IONS, COMPUTED BY SOLVING 
EQS. (HI-59), (lII-60). Electron and ion con- 
centrations are given in terms of the electron 
density N ea at the level z = 0, at which height 
the ionic composition is taken to be [o + ] = [He 4- ] = 
49 $>, [h + ] = 2 $ 0 , The level at which the He + and H + 
concentrations are equal is near z = 17. The unit 
of reduced height z is the scale height of 
neutral atomic oxygen. 
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6 . Solving the Continuity Equation 


Much of the progress in ionospheric theory has been achieved by 
obtaining solutions of the continuity equation and comparing them with 
observation. Of course, the full equation is so complicated, and its 
coefficients so poorly known, that drastic simplifications have to be 
made. In this section, we will start with the simplified equation 

Sltf/St = q - l(n ) - S(Nw)/Sh (lTI-62) 


Only the vertical velocity component w is included in the movement term 
and all losses are included in l(n). When discussing the variation of 
these quantities with reduced height, we should remember that the unit 
of z is the scale height of the ionizable gas. 

In earlier sections, we have already discussed the "Chapman, " 
"Bradbury, " and "transition" layers which apply to the situation of 
photochemical equilibrium, when the equation reduces to q = L(n). If 
time variations are included, but movements neglected, a first-order 
total differential equation is obtained, which can be solved analytically 
in some cases, otherwise numerically (Millington, 1932). 

If the movement term is reinstated, the equation becomes a partial 
differential equation in h and t. It may be simplified by following 
the motion of a "cell" of ionization, using the total derivative appro- 
priate to vertical motion, 

d/dt = (S/St) + w(S/Sh) (HI-63) 

in which we have specified (dh/dt) = w. Insertion of this operator into 
the continuity equation yields 

dN/dt = q - L(n) - NSw/Sh (lII-64) 

The usefulness of these equations lies in the fact that if Sw/Sh can be 
neglected, as may be the case for electromagnetic drifts, Eq. (III-64) 
reduces to a purely "photochemical" equation. However, as the altitude 
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of the cell changes, q and L will vary, so that the equation must 
in general be solved numerically . Use of this technique has been limited 
by the unsatisfactory state of the theory of electromagnetic drifts. 

The principle of M f ollowing-the-cell ” may also be applied to thermal 
motions. This is especially pertinent to the F region where large 
diurnal temperature changes occur. Suppose that the neutral atmosphere 
expands and contracts, so that the air moves vertically. The gas con- 
centration, which largely determines q, and the diffusion coefficient 
D, changes in a complicated way at any fixed height h, but within any 
particular cell, the concentration just varies inversely with temperature. 
This is obvious from the perfect gas equation, if we remember that the 
pressure within a given cell is constant, being the weight of the over- 

lying air. Then q, p and D are best written as functions of pressure, 
or of reduced height z. 

An approach of this type was used in a discussion of temperature 
effects on the equilibrium electron distribution by Garriott and Rishbeth 
(1963). It is found that, within certain assumptions, the shape of the 
electron density profile in terms of reduced height, e.g. n(z), is 
unchanged with temperature variations, but that the magnitude is altered 
in proportion to T ^ , These results have been extended in an approxi- 
mate manner to the time varying case by Rishbeth (^i964j. It is found 
that the temperature variations can so distort the N(t) curves that 
the maximum of N occurs before noon, at heights below the peak.. Shima- 
zaki (l957) has also used this technique when both electromagnetic and 
thermal motions are considered. The implicit assumption, that the moving 
air carries the ionization with it, is probably valid except near the 
magnetic equator. There will be complications if the temperature changes 
affect the composition at the turbopause of if large horizontal flows of 
air occur, but even in these circumstances a "f ollowing-the-cell " approach 
may be useful. 

When diffusion is introduced into the continuity equation, as must be 
done if the F2 peak is to be adequately treated, the continuity equation 
contains the diffusion operator jQ which involves ci/dh and d 2 /dh 2 . 
Equilibrium solutions of this equation are of some interest, even for the 
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F2 layer because it is true during much of the day that dN/dt is 
smaller than other terms in the continuity equation. At least they 
provide an insight into the relative importance of different processes. 

Equilibrium solutions of the diffusion equation, with electromagnetic 

drift, were obtained by Yonezawa (1956, 1958). Rishbeth and Barron (i960) 

solved the equation numerically in a number of cases, using Chapman-type 

Z' ■” kz 

production functions and exponential formulas, D oc e and e , for 

the diffusion and loss coefficients. Several values were used for the 
constant k. (We recall from Sec. III-4 that k = 2 and k = 7/4 if 
the loss reactions depend respectively on molecular oxygen and molecular 
nitrogen in a diffusively-separated atmosphere, and that k = 1 if all 
the atmospheric constituents are fully mixed.) From such calculations, 
the following generalizations about the behavior of a daytime equilibrium 
layer can be made: 

(a) The F2 maximum electron density occurs at a level where dif- 
fusion and loss are of comparable importance, i.e., where 

8 - D /h?, in which the subscript V refers to the maximum. 

*rn m' x 9 

(b) At the maximum, and below it, the electron density is approxi- 
mately given by N ~ q/j3, just as it would be in the absence of 
diffusion. 

(c) Well above the maximum, the electron density distribution is 
exponential, and takes the form N ~ e“ z /^ [as we obtained in 
Eq. (III-60)]. 

(d) These conclusions are substantially unaltered if gradients of 
scale height are present. 

(e) Vertical electromagnetic drift, magnitude w, alters the level 
of the maximum by an amount of order wH^/D m scale heights; 
is still given, very roughly, by (b). 

Points (a), (b), (c) are illustrated by the sketch in Fig. 19, More 

detailed discussions of "equilibrium layers" have been given by Bowhill 

(1962) and Nisbet (l963). The theory of the decay of the nocturnal F2 

layer has been investigated by Martyn (1956) and Duncan (1956), who used 

_ 2z 

the assumption {3 oc e Z , and by Dungey (1956) who used the form (3 <x e 
All these authors assume D« e Z , They investigate solutions which decay 
without change of shape, with an "effective" decay coefficient (3 5 ; i.e., 
of the form 
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FIG. 19. EQUILIBRIUM ELECTRON DENSITY DIS- 
TRIBUTION N( z ) FOR THE F2 LAYER. Below 
the peak N » q/p (photochemical equili- 
brium) and above the peak N oc exp (-1/2 z) 
(diffusive equilibrium). At the peak, the 
relations (3 = d/h? and q = PN are 
approximately valid. 31098 



N( z , t ) = N g (z ) • e ^ * 


(ill -65 ) 


The results of these investigations may be summarized as follows. 

There exist a number of ”shape-preservin g rr or ’’stationary ” distributions 

N (z) which are such that their peaks lie at heights h . At these 
s s 

heights, the loss and diffusion coefficients are connected by equations 

2 

of the type p^ = u^D^/h^, and the effective decay coefficient is given 
by equations p' = up where u and u are numbers of order unity. 
Usually we are only interested in the most slowly-decaying stationary 
layer which has the smallest value of P’. The special case investigated 
by Duncan and by Martyn, namely P « e Z , is of interest because it happens 
that u = u^ = 1, and also that the ’’shape preserving” function is the 
’’Chapman alpha layer” which arose previously in an entirely different 
context, namely, the equilibrium between photoionization and a square- 
law loss process (Eqs. III-29 and III-21). If an upward (or downward) 
electromagnetic drift is applied, then the peak of the layer is raised 
(or lowered), but the effective decay coefficient is still equal to the 
value of p at the maximum. In other cases, such as that studied by 
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Dungey (with {3 <x e , which is probably closer to the situation in the 

actual F2 layer), the analysis is more complicated and the numbers u^, 

u take different values. For instance, in the absence of vertical 
2 

drift, p = 1.86 P so that the decay of the layer corresponds to the 
value of p at some height below the peak of the layer. 

Time-varying solutions of the full diffusion equation, giving the 
theoretical diurnal variation of electron density, are naturally more 
difficult to obtain. Gliddon and Kendall (i960) were able to obtain an 

analytic solution, in terms of Green's functions, but had to make the 

~ z z 

assumption p * e , They used the upper boundary condition N oc e , 

which does not correspond to diffusive equilibrium. In a comprehensive 

paper (Gliddon and Kendall, 1962) they compared results obtained with 

-z/2 

this boundary condition and with the more usual assumption, N ~ e , 
and also discussed eclipse phenomena and the effects of vertical electro- 
magnetic drift, which they had to assume independent of height and time, 
Briggs and Rishbeth (l96l) constructed an electromechanical analogue 
computer to solve the diffusion equation. This method, though approxi- 
mate, has the advantage that the functions q(h,t), p(h) and D(h) are 
not restricted to special mathematical forms but may be derived from any 
empirical atmospheric model. The assumption of complete mixing (cor- 
responding to p a e Z in a mathematical model) was used initially but 
in later work (Rishbeth, 1963) p was assumed proportioned to the 
concentration (corresponding to P <x e 

These calculations provide a rather simple description of the diurnal 
behavior of a layer controlled by photoionization, linear loss and dif- 
fusion. At sunrise, the electron density near the peak increases at a 
rate which depends primarily on the production rate q, and diffusion 
and loss play a secondary role. The height of the peak, which corresponds 
to the F2 peak in the actual ionosphere, falls because of the rapid 
production of ionization in the lower F region, and reaches a minimum 
before noon. By this time, diffusion and loss become important and for 
a few hours the layer is not far from equilibrium, in the sense that 
dN/dt is small compared to other terms in the continuity equation. 

Under these conditions, the peak remains near a level where P ~ d/h , 
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as in the case of the "equilibrium layer" discussed earlier. Also, at 

heights up to the peak, the "equilibrium approximation" N = q/|3 is 

valid, and the maximum value of N lags after noon by a time ^l/p(h). 

Above the maximum, however, diffusion is so rapid that the variation of 

h T closely follows the variation of N . Later in the day, N decreases 

m m 

and h rises as solar control weakens, 

m 

After sunset, h approaches the "night stationary" level h and 
m s 

the layer develops into the "shape preserving" form which decays exponen- 
tially with time throughout the night, in the manner of Eq . (lII-60), 

The pressure level (or z coordinate) of the maximum electron density 
is higher than in the daytime. Duncan (1956) showed that when (3 e~ Z , 

the maximum occurs for 0 = 0.25 d/h 2 , while Dungey (1956) showed that 

-2z 1 / 2 

for p °r e , the level was at 0 as 0.11 D/H^ . However, this increase 

is more than offset by the thermal contraction of the atmosphere at 

night and the real height of the maximum would be expected to decrease 

(Rishbeth, 1964). 
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IV. MORPHOLOGY OF THE IONOSPHERE 


1 . D Region 

Until the last few years, most of the information about the height 
variation of electron density and collision frequency in the lowest iono- 
sphere was derived from the indirect evidence obtained from radio propa- 
gation experiments. These consisted of phase and amplitude measurements 
of low frequency waves (in the range 10 - 200 kc/s) which are reflected 
in the D region; and of measurements of absorption, polarization and 
cross-modulation of waves at higher frequencies (up to a few Mc/s) which 
are reflected in the E or F regions but traverse the D region. 

Some of these techniques were discussed earlier in Sec. II, 

The long-wave experiments give information about reflection heights, 
which vary with the solar zenith angle. Frequencies around 16 kc/s are 
reflected at 70 - 75 km by day but from the E layer by night (Bracewell 
et al , 1951 ), the transition occurring rather suddenly at sunset and sun- 
rise. The data have been interpreted in terms of one or more discrete 
layers of ionization below the E layer (Bracewell and Bain, 1952), 

Absorption measurements, at a series of radio frequencies, have shown 
a frequency dependence consistent with simple magnetoionic theory except 
near the critical frequencies f^E and f^Fl, when complications are to 
be expected. The absorption is found also to depend on the solar zenith 
angle X anci on t ^ ie mean sunspot number R. Appleton and Piggott (1954) 
found that at 4 Mc/s 


[log p| oc (Ch x) 

|log p| oc (l + 0.01 R) 


(IV-1) 


where p is the TT reflection coefficient" of the ionosphere and is a 
measure of the total absorption along the path. These measurements sug- 
gest a strong solar control of the D region, at least at the mid-latitude 
station (Slough, 52°) where the most extensive studies have been made. 

In winter, however, the absorption is much greater on certain days and 
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the long-wave and reflection coefficient measurements suggest that addi- 
tional ionization is present in the D region. The occurrence of this 
"winter anomaly" has been discussed by Thomas ( 1962a) and Dieminger 
(1952). Thomas finds that the excessive absorption occurs over geographi 
cal areas of about 1000 km in extent, but there seems to be no correla- 
tion with magnetic disturbance, except perhaps on a long-term basis. It 
may be that the ultimate cause of the anomaly lies at lower levels in 
the atmosphere. 

The only method so far devised that makes extensive use of the long- 
wave phase and amplitude data depends on "full wave" solutions of the 
wave propagation equations (Barron and Budden, I960) for a variety of 
models, combined with trial -and -error fitting of the parameters to the 
data. 

According to Nicolet and Aikin (i960) the chief production process 
in the upper D region is the photoionization of NO by Lyman a radiation 
(A = 1216 %). This radiation penetrates into the D region because it 
cannot ionize any of the gases found at higher levels, and neither is it 
strongly absorbed by them. Below 100 km it is absorbed by molecular 
oxygen (which it dissociates) but a fraction is expended in the photo- 
ionization of N0 f which has a rather low ionization potential. Although 
neutral NO has not been detected spectroscopically in the D region, 
Nicolet (l96l) considers that various chemical reactions can supply the 
trace (only 10 ^ of the total concentration) that is required. 

The model of Nicolet and Aikin, Fig. 20, shows the component of ioni- 
zation due to Lyman a is greatest between 70 and 80 km (for X = o). It 
is possible that a peak of electron density exists in this region, but 
in any case a positive gradient of electron density exists between 70 
and 75 km and this could account for the long-wave reflections observed 
at this level. However, it may be that the main contribution to D 
region ionization is made, not by Lyman a, but by hard X rays in the 
range 1 - 8 8. Above 85 km, the bottom of the E region production 
curve (due to softer X rays of 10 - 100 8) gives rise to another large 
electron density gradient, but this should really be considered as part 
of the E layer. 
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FIG. 20. CONCENTRATIONS OF POSITIVE 
IONS, NEGATIVE IONS AND ELECTRONS 
VS ALTITUDE FOR A VERY QUIET SUN 
(after Nicolet and Aikin, I960). 

During solar flares, the D region ionization is greatly enhanced, 
severe absorption is produced, and the reflection height of long waves 
drops by several kilometers. These events are called "sudden ionospheric 
disturbances," or SID's and will be mentioned again in Sec. VII. Since 
such an event starts at the same time as the visual flare and is observed 
at all latitudes, it is presumably due to electromagnetic radiation which 
must be more penetrating (i.e., possess a smaller cross-section) than 
that responsible for the normal D layer, since the height of peak 
photoionization is independent of its intensity but depends on the absorp 
tion cross-section a. 

Another type of disturbance is the "polar cap absorption" (PCA) which 
we mention here because of its interest in connection with D region 
photochemistry. It begins a few hours after solar flares and is thought 
to be initiated by the arrival of solar protons with energies on the 
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order of 10 Mev . It persists for several days with a pronounced diurnal 
variation, which is attributed to attachment and photodetachment processes. 
However, an alternative suggestion is that the disturbance creates extra 
neutral NO which can then be photoionized (Herzberg, I960), 

Other D region disturbances are probably associated with ionization 
by charged particles. This is especially true of the "blackouts" and 
"auroral absorption" observed at high latitudes. Ionization by electrons 
has been discussed by Rees (1963), and by solar protons by Webber (1962). 

Even in the normal D region, most of the production below 70 km is 
due to cosmic rays. This source differs from the solar radiation in that 
it acts by night as well as by day, but it is very latitude-dependent 
because of the effects of the geomagnetic field on the motions of the 
incoming charged particles. According to the curves given by Webber 
(1962), the rate of electron production by galactic cosmic rays increases 
downward throughout the D region, and indeed throughout the whole 
mesosphere. But the N(h) profile is quite different from the produc- 
tion rate q(h), because of the altitude dependence of the various loss 
coefficients, and in most published models N increases monotonically 
with height up to (at least) 70 km. However, Pierce - ( 1963) has suggested 

that rates of collisional detachment are such as to produce a nighttime 

2 -3 

"C region" peak of electron density, perhaps 10 cm , at about 50 km. 

The existence of this peak had also been suggested on the basis of long- 
wave radio propagation (Moler, I960). 

In the lower D region, ionization has little effect on radio propa- 
gation, partly because the free electron density is small (N « N « N ) 

v e - + ' 

and partly because the collision frequency greatly exceeds the plasma 
frequency. Therefore, the ionization below about 50 km is generally 
considered to belong to the field of atmospheric electricity rather than 
ionospheric physics. 

Many of the papers already cited show values of the negative-ion 
ratio A = N_/n^ derived from theoretical calculations. Most authors 
find that the level where A = 1 occurs at about 70 to 81 km by day, 
and about 90 to 100 km by night. Some comparisons have been made by 
Hultqvist (1963) but there is still considerable disagreement on this 
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topic. The behavior of the D region at sunrise and sunset is very 
relevant to the question of ion composition. Long-wave reflection and 
absorption experiments show that the electron density starts to increase 
about an hour before ground sunrise (Bracewell et al , 1951 ). Although 
the height at which this ionization appears may not be well known (espe- 
cially in the case of absorption experiments), the fact that it is pro- 
duced when X > 90 deg implies that the effective radiation has traversed 
lower levels in the atmosphere. Aikin (l96l) showed that Ly CC and X 
rays do not penetrate the D region until near ground sunrise, so that 
the initial increases are attributed to photodetachment from negative 
ions. If the "screening height "- -the lowest height traversed by the 
radiation - -was well known, it would be possible to deduce whether visible 
or ultraviolet light is involved. One expects visible light to be cut 
off only by the solid earth (or by clouds in the lower atmosphere), whereas 
UV light is screened by the ozonosphere. 

This problem is acute in connection with the "polar-cap absorption" 
previously mentioned, because the values of X an ^ the related screening 
heights are inconsistent with the simple theory that electrons are pro- 
duced by photodetachment from 0 or 0 ions, which can be accomplished 
by visible light. Because of this difficulty, it has been suggested that 
the ionization responsible for polar-cap absorption lies below 55 km, in 
which case the observations would be consistent with screening by the 
solid earth (Hultqvist and Ortner, 1959). An alternative suggestion is 
that the negative ions are not 0~ or 0 as originally supposed, but some 
other ion (perhaps 0~ or NO") which has a high electron affinity, such 
that ultraviolet light is required for photodetachment. The screening 
height for this radiation would then be consistent with the data (Reid, 

1961 ) but this subject is still controversial. A recent review of D 
region processes is that of Reid (1964). 

2. E and FI Regions 

The E and FI layers of the ionosphere are generally regarded as 
being fair approximations to the idealized Chapman layer. Although the 
Chapman theory does provide a first-order description of the behavior of 
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the critical frequencies, more detailed study indicates that the layers 
are more complex. There is little doubt that the layers are within the 
"photochemical regime" described in Sec. III-4, although transport pro- 
cesses can produce appreciable perturbations. 

In applying the theory of photochemical processes to the E and FI 
layers, it is generally satisfactory to assume dN/dt = 0. For a simple 
Chapman layer in an isothermal atmosphere, we found in Sec. Ill -2 that 
the critical frequency is related to solar zenith angle X by the 
equation 


f Q = 9000 [(q 0 /a) cos x ]°- 25 (c/s) ( IV-2 ) 

where is the peak production rate for overhead sun (x = 0) and 

a the recombination coefficient, expressed in cgs units. Many investi- 
gators assume a relation of the type f oc (cos x) n and find the value 
of n as X varies with time of day or w r ith season, at any one station; 
or as X varies with latitude, for a number of stations. A study of 
diurnal variations (Tremellen and Cox, 1947) gave n = 0.3 (or somewhat 
more) for diurnal variations, but seasonal variations give a result 
closer to the theoretical value 0.25. A possible reason why n > 0.25 
in the E layer was mentioned earlier; that is, the existence of a 
positive scale height gradient, T = dH/dh. In Eq . (lll-3l), we deduced 

that n = 0.25(l + T), so that the diurnal variation of f E implies 

o 

that r zz 0.2, reasonably consistent with models of the thermosphere. 
However, there seem to be seasonal and geographical variations which 
cannot yet be explained in detail. 

The dependence of f E on the mean sunspot number R is studied by 
assuming Eq. (lV-2) to hold; it is then found that the square of the 
"noon” maximum electron density is given by 

q Q (E)/a(E) ft* 180 (l + 0.010 R) x 10 8 cm” 6 (lV-3) 

(Allen, 1948). The height of maximum production for x = 0 for the E 
layer has been found by Robinson (1959) to be about 108 km. 
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The recombination coefficient a has been determined to be greater 

than 10 -8 cm 3 /sec from eclipse observations (Ratcliffe, 1956b); and by 

measuring the "sluggishness," i.e., the interval by which the time of 

greatest N E lags behind local noon, which should correspond to 
in 

l/2 aN and is of order 10 min (Appleton, 1937, 1953). From the decay 

/x -8 3/ 

of the E layer at night, Titheridge (1959) finds a = 2 X 10 cm /sec. 
However, in all these analyses, the E region was assumed to contain 
only a single ionic constituent. It may be necessary to consider the 
presence of two ions with different recombination coefficients to satis- 
factorily explain the observations (Bowhill, 1961 ). 

As regards the latitudinal variations of f Q E, there are departures 
from the normal Chapman layer behavior which may be caused by electro- 
magnetic movements associated with the quiet-day (S^) currents flowing 
in the E region. This explanation is favored by Beynon and Brown 
(1959); moreover, Appleton, Lyon and Pritchard (1955) find a systematic 
latitude variation of the post -noon lag of N^E which seems to be due to 

movements rather than to changes in the "sluggishness time l/2 QN. 

We must point out that the vertical structure of the E layer is 

less simple than the previous discussion might imply. On ionograms, two 
or more cusps often appear, and so it may be difficult to* assign a 
definite E layer critical frequency. If a sequence of ionograms is 
studied, however, it is usually found that one particular cusp displays 
a more consistent diurnal variation than the others, and this is usually 
tabulated as the "actual" critical frequency. Occasionally a second 
cusp, at a higher frequency, is found to vary in a consistent way, and 

it is then denoted as f E2. Additional stratifications of a more or 

o 

less temporary nature are often observed (Becker and Dieminger, 1950, 
Dieminger, 1959; Robinson, 1959). 

The term "sporadic E, " or Es, is applied to E region ionization 
which does not behave in any regular manner. It is generally considered 
that this phenomenon is not closely related to the ordinary E layer 
produced by ultraviolet radiation, and we shall deal with it in a later 
section (lV-7). 

Since the E layer is produced, not by monochromatic radiation but 
by a variety of wavelengths (including X rays and ultraviolet lines 
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such as Lyman p) f it is not surprising that a complex structure exists. 

So far, however, the data are insufficient for any correspondence 
between individual spectral lines and ionogram features to be established. 
M though conventional ionograms do not give very detailed information 
about ionization between the E and FI layers, rocket N(h) profiles 
have shown that there is little (if any) decrease of N above the E 
layer, that is, there is no pronounced "valley. " 

For the FI layer, the solar cycle variation is quite well described 
by the formula 


Q Q /ct' = 500 (l + 0.016 R) x 10^ cm ^ (lV-4) 

(Ratcliffe and Weekes, I960). Eclipse and other observations, summarized 

by Ratcliffe (l956b), give value of a in the range (l/2 to l) lo’ 8 cm 3 sec 

The height of maximum production h^ for the FI layer is very difficult 

to determine from ionogram data. For summer noon at Slough, Thomas, 

Haselgrove and Robbins (1958) found h FI to lie at about 185 km at 

m 

sunspot minimum (1953) and 190 km at moderate solar activity (1950). 

These heights would be expected to lie near h , or a fraction of a scale 
height above it. 

Although the FI ledge is not always observed on ionograms, its 
behavior approximates that of a Chapman layer when it does appear. 
According to Allen (l948), the critical frequency f FI varies diurnally 
as (cos x) * • This is different from the behavior of f E f mentioned 

earlier, but is consistent with the theory that the FI ledge lies at 
the transition between regions of "quadratic" and "linear" loss laws. 

The index n would be expected to be less than the value 0.25, because 
the "transitional" recombination coefficient ap/(aN + (3) decreases 
upward. There is evidence of some geomagnetic control of the FI layer 
(Cummack, 1961 ), which might possibly 7 result from electromagnetic move- 
ments such as we discussed in connection with the E layer. 

3 . F2 Region Problems 

Early investigations of the F2 layer critical frequency revealed 
that it does not behave at all in the manner of a simple Chapman layer. 
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Although we now have available perhaps 10 values of together with 

many computed electron density distributions for selected stations, the 
physics and chemistry of the F region are only partially understood, and 
the rates of the photochemical processes are known to within only a factor 
of two, at best. Some of the outstanding problems of F2 layer morphology 
are: 

(a) behavior of the F2 peak (N m F2, ' h m F2) and the diurnal behavior 
of electron densities at fixed heights- 

(b) the topside of the F layer 

(c) the seasonal and annual anomalies 

(d) the equatorial anomaly 

(e) high latitude behavior 

Each of these problems will be discussed in turn and some of the 
proposed explanations advanced, even though they may be only tentative 
suggestions. Little mention of electromagnetic drifts will be made, 
although these may be responsible for some of the anomalies. Mainly, 
this omission is made because of the great difficulty in obtaining quanti- 
tative measurements of electromagnetic drifts. The theory of drifts is 
described in more detail in Sec. VI. 

We should state also that all our remarks apply to the average 
behavior of the undisturbed F2 layer. More than the other layers, the 
F2 layer is subject to day-to-day variations of perhaps 20 percent, and 
may even change appreciably from one hour to the next. Some discussion 
of this variability is given in a review of F2 layer morphology by 
Wright (1962). 

(a) The F2 peak and N(h,t) curves . In Figs. 21a and 21b, we 

reproduce the variations of N F2 and h F2 for quiet days at Slough 

r mm 

for a number of months of low solar activity (1953) and moderately high 
activity (1950), given by Thomas, Haselgrove and Robbins (1958). These 
display some of the anomalies we shall discuss. Figure 21a shows the 
seasonal anomaly (for Jan 1950) and also reveals that N^F2 f a i- r ly 

sharply peaked about noon in winter but is rather irregular in other 
seasons . 
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a. N m F2. The six curves show the seasonal and 
solar cycle variations. 



31277 

b. h F2 (after Thomas et al) 
m 7 

FIG, 21. AVERAGE DIURNAL VARIATION FOR THE 
10 INTERNATIONAL QUIET DAYS IN EACH MONTH 
AT SLOUGH. 


The N(h,t) curves of Fig. 22 show many features typical of the 
mid-latitude F2 layer. Up to 200 km, the curves are fairly regular, 
consistent with control by photochemical processes at these heights. 
Higher up, where transport processes are important, the curves are less 
regular and not symmetrical about noon. Although we should expect the 
maximum of N to occur at a time about l/{3 later than noon, it seems 
that movements so affect the layer that the maximum N may be attained 
before noon. This may well result from thermal expansion of the F 
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FIG. 22, n( t ) CURVES, SHOWING THE 
AVERAGE DIURNAL VARIATION OF ELEC- 
TRON DENSITY AT FIXED ALTITUDES 
(in kilometers) FOR MAGNETICALLY 
QUIET DAYS DURING ONE MONTH AT 
SLOUGH. (After Croom, Robbins 
and Thomas, 1959a.) 


region. As the F region cools and contracts in the evening, there is 
a tendency for N to increase and on occasions both morning and after- 
noon maxima of N at fixed heights, or of N , are observed. 

The heating and cooling of the atmosphere alters the height at which 

2 

d/h « 8, near which h F2 should occur. Although the drop of h F2 
at sunrise, and its subsequent rise (Fig, 21b), are broadly consistent 
with the theoretical behavior outlined in Sec. III-6, the difference 
between noon and midnight h^F2 (generally 100 km) is much greater than 
would be expected from theory, especially when the day-night temperature 
changes are allowed for (Rishbeth, 1964). At lower latitudes, the dif- 
ference is smaller, and changes sign near the equator where h^F2 is 
higher by day than by night. These facts are unexplained, though some 
kinds of movements might be responsible. 

At night, N F2 decreases rather irregularly. This decrease may not 
’ m 

be sustained throughout the night, especially on long winter nights when 

N F2 appears to reach a "base level." The N(t) curves of Fig. 22 
m 

illustrate this. N may then fluctuate irregularly but does not show 

m 

a steady decrease. In fact, there is sometimes a slight increase well 
before F layer sunrise. 
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The base level is also evident in the behavior of the total electron 
content of the whole ionosphere, as determined by lunar-echo Faraday 
rotation observations (Evans and Taylor, 196l). Again, the explanation 
may be an upward drift of the ionosphere due to electromagnetic forces, 
or it might be due to corpuscular ionization, such as the flux of soft 
electrons postulated by Antonova and Ivanov-Kholodnyy (l96l). An alter- 
native source might be provided by slow downward diffusion of ionization 
from the exosphere, as investigated by Hanson and Ortenburger (l96l), but 
theoretically this appears to be inadequate to maintain the F layer 
(Hanson and Patterson, 1963). Also, it seems unlikely that photoioniza- 
tion by EUV could provide enought ionization at night, since an attempt 
to detect He II 304 8 radiation at night yielded negative results (Byram 
et al , 1961). 

(b) The topside of the F layer has been investigated in many 
rocket experiments. Most of these yield an exponential N(h) profile 
roughly consistent with diffusive theory, and give useful information 
about the variations of temperature in the F region, though the results 
are sometimes anomalous (Bowhill, 1963). 

One of the interesting discoveries from the UK-1 satellite ’"Ariel" 
is that f i eld -al igned "ledges" exist in the topside F region (Sayers, 
Rothwell and Wager, 1963). One of these seems to be associated with the 
"troughs" of the F region equatorial anomaly, while the others are 
found at higher latitudes. Further discussion of these ledges, based on 
"Alouette" data, has been given by Lockwood and Nelms (1964). 

The "Alouette" topside sounder has provided a powerful tool for 
investigat ing the upper F region. So far, only preliminary results 
have been published, including a number of topside N(h) profiles and 
some latitudinal distributions (Warren, 1963; King, 1963; Knecht and 
Van Zandt, 1963). Plots of log N(h) sometimes show a change of slope 
between 600 and 800 km, which is thought to represent the transition 
between 0 + ions and He + ions (Calvert, Rishbeth and Van Zandt, 1964). 

Not all the profiles, however, can be fitted by a simple model of this 
type (Bauer, 1964). 
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c- Seasonal and annual anomalies . These anomalies have been known 

for a long time. The most detailed analysis of seasonal and annual 

variations of N F2 has been made by Yonezawa and Arima (1959) and 
m 

Yonezawa (1959) , Over the world as a whole , there exists an annual 

variation of N F2, which is about 20 percent greater in December than 
m 

in June. This may be compared with the annual 3 percent variation of 
the earth-sun distance, which is a minimum in January. In addition, 
there is a "seasonal" anomaly in that daytime NJF2 tends to be greater 
in winter than in summer. In the northern hemisphere these two effects 
act in phase to give a resultant anomaly, which is very large at high 
solar activity (see data for 1950 in Fig. 21a) but scarcely evident at 
sunspot minimum (as for 1953). Croom, Robbins and Thomas (i960) have 
shown that this effect is most pronounced within a certain range of 
latitudes (probably magnetic latitude) and the solar-cycle dependence 
has been further studied by Thomas (1963). 

Rishbeth and Setty (l96l) found that, just after layer sunrise, 
dN/dt in the F2 layer is much greater in winter than in summer; and 
that the increase of N begins at a larger value of X (about 97 deg) 
in winter than in summer (about 93 deg). This seasonal sunrise anomaly 
can be seen in the slopes of the N^t) curves of Fig. 21a and is 
particularly interesting because at sunrise the continuity equation is 
dominated by the photochemical terms, especially q. These authors sug- 
gested that seasonal changes of atomspheric composition--in particular, 
the ratio nfo]/n[N ]--could affect q and p. 

Such changes could occur in the F region if there occurred a change 
in the o/n^ ratio at the turbopause near 100 km, where diffusive sepa- 
ration of neutral gases begins (Sec. 1-4); or if the altitude of the 
turbopause itself varied. These ideas have been developed by G. A. M. 
King (1961), who suggested that travelling disturbances might affect the 
height of the turbopause; and by Wright (1963, 1964). Recently, it has 
been suggested that large-scale meridional circulations of air in the 
mesosphere and lower thermosphere give rise to changes of the o/Ng ratio 
at the turbopause, and hence also in the F region (King, 1964). This 
can be related to seasonal changes in mesospheric temperature (Young and 


SEL-64-111 


106 


liiiiiJililliiiii 




Epstein, 1962); here we may have an important connection between the 
ionosphere and lower atmosphere. We note that no satisfactory explanation 
of the F region seasonal anomalies in terms of ionospheric movements, 
or in terms of temperature changes, has ever been devised, 

(d) The equatorial F2 layer . The F2 layer at low latitudes is 
very peculiar. Sometimes, electron densities are greater at midnight 
than at noon. Norton and Van Zandt (1964) succeeded in explaining day- 
time N(h,t) curves in terms of photochemical processes only, taking 
into account the diurnal temperature changes in the atmosphere. Vertical 
diffusion is neglected because ionization cannot diffuse across the 
geomagnetic field lines. It can, however, diffuse along the field (i.e., 
horizontally) and this may affect the latitude distribution of ionization. 

A plot of noon values of N F2 as a function of latitude, shows a 
pronounced "trough" centered on the magnetic dip equator, with "crests" 
about 15 deg or 20 deg to the north and south (Appleton, 1946), This 
trough is also found in plots of N at fixed heights below (Croom, Robbins 
and Thomas, 1959b) and above (King, 1963) the F2 peak. Studies by 
many authors, most recently Lyon and Thomas (1963) show that the anomaly 
exists during most of the day, being most pronounced around sunset, but 
disappears after midnight. The anomaly shows rather different features 
in different longitudes. 

-^It was suggested by Mitra (1946) that the anomaly might be due to 
the diffusion of ionization away from the equator, causing an accumulation 
of electrons to the north and south. To test this suggestion theoreti- 
cally, it is necessary to use a form of the diffusion equation which 
takes account of the geometry of the geomagnetic field. Solutions of 
this equation, for equilibrium conditions, lead to the conclusion that 
diffusion is not sufficient to produce a trough as large as is actually 
observed and that other processes are responsible (Rishbeth, Lyon and 
Peart, 1963), though this subject is very controversial (Goldberg, Kendall 
and Schmerling, 1964). It is also quite possible that electromagnetic 
drifts contribute to the equatorial anomaly (Duncan, I960). 

(e) High latitudes . The morphology of the F2 layer at high lati- 
tudes is very complicated. Probably transport processes play a larger 
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part than at lower latitudes, especially during the polar night when the 

p layer persists despite the absence of solar photoionization. Observa 

tions of f F2 at the South Pole (Knecht, 1959) and Halley Bay (Bell- 
o 

chambers and Piggott, I960) show unusual variations of F2 layer param- 
eters, for which no complete explanations have been given. 

4 . F Region Rates 

In principle, it should be possible to determine the rates of pro- 

duction and loss in the F region from the large amount of available 

N(h,t) data. But in practice it has proved difficult to assign reliable 

numerical values to these rates. The reason for this seems to be that 

the F2 layer is so complex that it is difficult to isolate the effects 

of any one process. In the lower F region, photochemical processes 

dominate and the situation should be simpler, but it turns out that the 

best data give only the ratios of quantities such as q, a, (3 and D, 

not their absolute values. This situation is now being resolved by 

rocket observations of the solar spectrum, which give absolute values of 

q. Even these, however, are influenced (to some extent) by uncertainties 

about atmospheric composition. Another approach, of course, is the 

laboratory determination of rate coefficients. It cannot be yet said 

that these are sufficiently accurate to provide definitive values. In 

some cases, there is conflict with the ionospheric evidence and it may 

be that laboratory and ionospheric conditions are not really comparable. 

Ionospheric data relating to various process rates are considered below. 

Noon electron density. We saw previously that observations of the 

FI layer lead to the relation (N^Fl) sec X = ( c l o / 0C ) = 500 0 + 0- 16 R ) 

X 10 8 cm” 6 (Eq. IV-4). For the F2 layer, at heights up to the peak, 

we have suggested that N s q/p, but owing to the variability of the F2 

layer, no general formula for N^F2 can be given. Allen (1948) finds 

that at several stations N m F2 « (l + 0.02 r), and this linear dependence 

incidentally supports the theory that the rate of loss in the F2 layer 

2 

is given by (3N and not by 0(N . 

Night decay. Values of P for the F2 layer have been deduced 
from nighttime N(t) variations. After sunset, the decay of N at 
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a ny height h should at first depend on the local loss coefficient p(h). 
Once the nighttime 'stationary" layer is established, the decay would 
depend on the value of £ at the peak (Sec. III-6), but electromagnetic 
and thermal motions apparently complicate the analysis. Ratcliffe et al 
( 1956) suggested the empirical formula (with h in km) 

0(h) = 10' 4 exp [(300 - h)/50] sec ^ (lV-5) 

from a study of N(h,t) data from several stations. Owing to the dif- 
ficulty of obtaining acceptable nighttime N(t) curves for the lower 
F region, values of a have not been deduced for the FI layer by 
this method. As we have discussed, it now appears that the loss coef- 
ficient is quite sensitive to temperature variations, since the neutral 
concentrations at any fixed height vary with temperature. This effect 
has been studied by Nisbet and Quinn (1963) and others even more recently. 

Eclipses . Since this topic is discussed in Sec. IV-5, we only mention 
here that the observations do not lead to unambiguous values of a for 
the FI layer. Similar remarks apply to methods based on the classical 
Chapman layer theory (Appleton, 1937). In the F2 layer, the situation 
is usually hopeless because of movements, although in some cases this 
difficulty appeared not to be serious (Van Zandt et al , I960). 

Sunr i se . Just after F layer sunrise, the peak of q lies in the 
F2 layer and we expect to have dtf/dt « q(x) * q Q /ch X- Rishbeth and 
Setty ( 1961 ) found values of q(x) but, owing to uncertainties in Ch x 
and in other terms, this method does not lead to accurate values of q 

o' 

From these methods, we can obtain values for production and loss 
rates, but these are probably only accurate to within a factor of two. 
Values of the ambipolar diffusion coefficient D have to be deduced 
indirectly, since direct measurements have not been made and the theory 

f_ _ _ 

is not entirely satisfactory. Shimazaki (1964) has obtained values of 
D and £ from detailed analysis of night N(h,t) data. His values of 
D are considerably smaller than those needed to explain the behavior 
of the daytime F2 layer. 
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We must emphasize that the coefficients q, p and D depend very 
much on neutral air densities, which in turn depend on scale heights. 
Hence, they must vary considerably with solar activity, and there is no 
such thing as the value of P at any given height. Some illustrative 
calculations have been made by Rishbeth (1964), whose suggested values 
of P at 300 km, as a function of solar activity and F region tempera- 
ture, are compared with other determinations in Figs. 23a, b. The noon 
value of q at 300 km increases very rapidly with increasing solar 
activity, because the incident flux of radiation and the concentration 
of ionizable gas, atomic oxygen, are both increased. Values of q, P 
and D are tabulated in Table III. Theoretical N(h,t) curves calcu- 
lated with these coefficients agree with many daytime F2 data, but in 
order to make them consistent with the survival of the F2 layer through- 
out the night, it has to be assumed that some nocturnal source of ioni- 
zation exists. 

5 , Eclipse Effects 

On those infrequent occasions when the sun is partially or totally 
occulted by the moon, first order perturbations are observed in the 
electron density distribution in the ionosphere. This is particularly 
true in the D, E and FI regions, but sometimes they also extend up 
to the F region peak. These events provide useful methods for the 
calculation of the production and loss terms in the continuity equation 
and provided early evidence that the E and FI layers are produced by 
elec tromagnetic- -not corpuscular- -radiation . Although absorption and 
therefore electron density varations in the D region have been observed 
during eclipses, most studies relate to the higher layers. 

When considering the E and FI regions, we may start with the 

2 

appropriate continuity relation c)N/dt = q - CiN , which we found in an 
earlier section. The variation of q^ with time on a non-eclipse or 
"control" day is sketched in Fig. 24 as a solid line. From Chapman's 
theory, q = q (cos x) for an isothermal atmosphere, but on the eclipse 
day we must multiply this by a factor <>(t) equal to the uneclipsed 
fraction of the sun's disk, assuming the ionizing radiation is emitted 
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O: Norton and Van Zandt (1964), equatorial F2 layer model, t = 12 hr 

©: Norton et al ( 1963), model of E and F layers, t = 10 hr (below) and 12 hr (above) 

#: Van Zandt et al (i960), eclipse observation, t = 09 hr (below) and 12 hr (above) 

(in the last two cases, the values for 12 hr are rough estimates, obtained by connecting the values 
given for the earlier hours for changes in atmospheric density between those hours and noon.) 




TABLE III. ROUNDED VALUES OF F2 LAYER PARAMETERS 
(from Rishbeth, 1964) 


Parameters 

Sunspot 

Minimum 

Sunspot 

Maximum 

Unit 

Zurich sunspot no. R* 

0 

200 


Solar 10.7 cm flux S* 

70 

200 

-22 - 2 / / v -1 

10 wm (c/sj 

Values at 300 km: 
Temperature T 
(noon ) 

840 

1610 

°K 

(Harris and Priester, 




midnight ) 

580 

1180 

°K 

Photoionization q 
(noon) 

50 

750 

-3 -1 

cm sec 

Loss coefficient p 
(noon) 

0.4 

7 

_ -4 -1 

10 sec 

(midnight ) 

0.06 

4 

_ -4 -1 

10 sec 

Diffusion coefficient D 
(noon) 

7 

2 

in 10 2 - 1 
10 cm sec 

(midnight) 

14 

2 

- n 10 2 - 1 
10 cm sec 


Both these parameters are used to describe solar activity. For 
present purposes, a linear relation between them is assumed, and is 
defined by the values quoted. 


uniformly across the disk. Thus q = $q^(cos x) t an< ^ this function, as 

well as the maximum electron density N , are shown in Fig. 24 for a 

m 

hypothetical total eclipse. 

At second contact t^, the production function is zero and we may 
calculate the recombination rate as 


5n /St 

nr 



m 


(IV-6) 
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in which the implicit assumption has been made that the height of N 
does not change during the eclipse. At the E region peak, Minnis 
(1955) has found d ^ 1.5 x 10 cm^ sec * using this method. 



ON A CONTROL DAY AND THE PRODUCTION RATE AND MAXIMUM 
E REGION ELECTRON DENSITY ON THE DAY OF A HYPOTHETI- 
CAL TOTAL ECLIPSE. From curves similar to these, it 
is possible to calculate the recombination coefficient. 

It is also possible to assume a value for a and predict the shape 
^m VS ^ me ^ or the ec lipse, looking especially for the time sepa- 
ration between third contact and t , the electron density minimum. 

The time difference is expected to be a few minutes, but the observations 
do not show any very consistent pattern ( Ratcliff e, 1956). 

It now appears likely that a has been somewhat underestimated in 
the eclipse analyses because some of the X ray radiation may originate 
in the solar corona and thus fail to be occulted by the moon. Elwert 
(1958) has calculated that 10 to 20 percent of the X ray radiation 
may come from above the photosphere. Since q is usually a much larger 
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term than Sn/cH in the E region continuity equation, the electron 
production remaining at the time of total optical eclipse would generate 
a much larger proportional error in OL . Furthermore, the radiation is 
probably more intense from various active regions on the sun and the 
production function would not be the smoothly varying geometrical func- 
tion originally assumed. Values for a between 10 and 10 cm /sec 
are usually adopted for the E region, and are somewhat smaller in the 
FI region. Thus eclipse results fall in the lower portion of this range. 

The various eclipse measurements of a become more consistent with 
each other and also with the results from entirely different methods 
(e.g., nocturnal decay and f, sluggi shness” ) when a two-ion mixture is 
considered (Bowhill, 196l). If the ions have distinctly different values 
of a, one ion will decay more rapidly, resulting in a diminishing recom- 
bination coefficient with time. Bowhill suggests that N0 + and 0^ are 
these ions and that the latter may decay most rapidly. After sunset, the 
NO + /o^ ratio should increase and this seems to be borne out by the rocket 
ion spectrometer results of Fig. 17. 

Eclipse results in the F2 region are more difficult to interpret 
than those at lower altitudes. Near the peak, little or no effect may 
be observed and movements resulting from diffusion or electromagnetic 
drift may also be important. The most successful analysis has been made 
by Van Zandt, Norton and Stonehocker (i960) for the October 1958 eclipse 
at Danger Island, They made several simplifying assumptions which appear 
to be rather well satisfied by the results. They have first assumed that 
temperature changes and all movement terms are negligible. Since Danger 
Island is at a magnetic latitude of only 11 deg, vertical diffusion is 
certainly slow. Also, their calculations have been made only at alti- 
tudes greater than 280 km, where the loss of electrons is believed to be 
linearly related to the electron concentration (£N). Furthermore, these 
altitudes are all well above the production peak and have small optical 
depths. Thus, the production rate at any fixed height is approximately 
constant without dependence on cos X f° r control days; on the eclipse 
day, q need only be multiplied by the uneclipsed portion of the sun’s 
disk. The radiation principally responsible for F2 ionization is in 
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the UV, and it is believed to originate more uniformly across the disk, 
rather than partly in active regions and in the corona. This eliminates 
one of the principal difficulties encountered at lower altitudes. Within 
these assumptions, the continuity equation for the F2 region can be 
written 

= q(h) 

Since q(h) is relatively constant well above the production peak 
as discussed earlier and 0(h) should be nearly constant if temperature 
changes are negligible, Eq. (lV-7) is of the form y = qx - 0. The 
values of y and x can be determined experimentally as a function of 
time at any fixed altitude, and we may expect a graph of y vs x to be 
a PProximat ely linear, with a slope of q(h) and y intercept equal to 
p(h). Van Zandt et al (i960) have made these calculations at each 10 km 
knight interval between 290 km and 450 km and obtained rather good con- 
sistency throughout the range. Their results may be expressed mathemati- 
cally as 

q(h) = 880 exp [(300 - h)/l86] elec/cm^/sec (lV-8) 

0(h) = 6.8 X 10 4 exp [(300 - h)/l03] sec" 1 (lV-9) 

where h is measured in kilometers. 

Although this value of 0 is much greater than that obtained by 
Ratcliffe et al (1956) from nighttime N(t) curves at times of moderate 
sunspot numbers, Fig. 23 shows that these two results are reasonably 
consistent. At the time of the Danger Island eclipse, the sunspot number 
was almost as high as has ever been recorded; hence, as discussed in 

.. . ..... . _ . g 

Sec. IV-4, large values of 0 are to be expected. If the scale height 
in Eq. (lV-8) is associated with atomic oxygen, and that in (lV-9) with 
molecular nitrogen or molecular oxygen, a temperature of over 3000 °K 
is deduced, but this would seem to be rather excessive, since other 
evidence favors values of 2000 °K to 2500 °K at sunspot maximum. 
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6, Ionospheric Irregularities 


The preceding sections dealt with only the large-scale or "background" 
distribution of ionization. But, superimposed on this background, small- 
scale irregularities of ionization seem to exist at every level in the 
ionosphere. The basic properties to be measured are their size and shape, 
intensity (i.e., fractional deviation of electron density) and their 
direction and speed of motion. The measurements are, however, liable to 
be influenced by instrumental selection. For Instance, a spectrum of 
irregularities may exist at a particular level, of which only certain 
components will be recorded by any particular equipment. Some informa- 
tion about irregularities is given by standard ionograms or h T (f) curves, 
and by continuous recordings of phase heights at fixed frequencies, or 
h r (t) curves. These are particularly useful for the detection of large 
"travelling ionospheric disturbances," to be discussed later in this 
section . 

Irregularities are also detected because of their ability to scatter 
radio waves. Signals at frequencies of 50 to 100 Mc/s, too high to be 
reflected from a smooth ionosphere even at oblique incidence, may be 
received at distances of a thousand kilometers from the transmitter, 
having been scattered from irregularities in the D region at 80 to 90 
km (Bailey, Bateman and Kirby, 1955). This forward scatter is of prac- 
tical importance, but most of the scientific study of the irregularities 
themselves has been conducted with more specialized techniques. One of 
these is the reflection method of Mitra (1949), which uses a pulse 
transmitter at a frequency of a few Mc/s, and three receivers in a tri- 
angle, spaced roughly one radio wavelength apart. Another method makes 
use of a radio star as the source and amplitude fluctuations produced 
by the ionosphere are again observed with spaced receivers. This latter 
method detects mainly irregularities in the F region, whereas the 
reflection method can be used to study the normal E layer by day, and 

has been applied to E ionization at night. If a suitable radio f re- 

fa 

quency is chosen, it can also be used to study the F layer. If the 
fading is caused by the steady drift of irregularities in the ionosphere, 
then the records from each receiver are similar, but there are systematic 
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time delays between the appearance of particular features (such as maxima 
and minima of signal strength) on the different records, from which the 
drift velocity of the irregularities can be determined. 

This simple "time-delay " approach is only valid if the ionospheric 
irregularities possess certain statistical properties. It has largely 
been superseded by a more sophisticated analysis which treats fading sta- 
tistically (Briggs, Phillips and Shinn, 1950; Phillips and Spencer, 1955; 
Bowhill, 1961). First, it should be noted that the wavefield at the 
ground is formed by diffraction from an irregular "screen" in the iono- 
sphere. This "screen" consists of irregularities which impose phase 
variations on the wavefront of the downcoming radio wave. If the varia- 
tions are "deep, " in the sense that the phase deviations exceed one 
radian in the width of one Fresnel zone, the pattern at the ground is 
smaller in scale than the screen, but it may retain the general shape 
and orientation of the irregularities (Hewish, 1951; Fooks and Jones, 

1961 ). For a "shallow" screen, in which the phase changes are less 
than a radian in a horizontal distance of one Fresnel zone, the ampli- 
tude pattern on the ground should be very similar in size and shape to 
the structure of the irregularities in the ionosphere. The distinction 
between these two cases therefore depends on both the wavelength and 
the altitude of the irregular! ties . The velocity of the wave pattern 
at the ground is equal to the drift velocity in the ionosphere in the 
radio-star scintillation method and twice as great for the reflection 
method . 

To describe the shape and size of the irregularities in the screen, 
use is made of the autocorrelation function of the electron density dis- 
tribution. In the two-dimensional case, this function may depend on 
direction and it is generally assumed that the irregularities can be 
described in terms of a "characteristic ellipse, " whose radius in any 
direction is the distance in which the correlation falls to the value 
l/2. The shape, size and orientation of this ellipse provides the infor- 
mation about the ionospheric irregularities. 

Since the pattern of the irregularities may change as it moves, another 
parameter V c has been defined which is zero if the pattern drifts without 
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change. If V c is comparable to the drift velocity V of the pattern, 
it implies that random changes in the pattern contribute to the observed 
fading as much as does the bodily drift, and this appears to be the situ- 
ation usually encountered. 

All these quantities can be determined from sets of fading records by 
a complicated procedure known as "full correlation analysis." This has 
been applied by Fooks and Jones (l96l) who compare the results with those 
of a simpler analysis developed by Yerg (1956) and the straightforward 
"time-delay " approximation described earlier. They further discuss the 
circumstances under which the simpler methods can give meaningful results. 
The time-delay method seems to lead to overestimates of drift velocity. 
Moreover, if the irregularities are elongated, the simple time-delay 
method tends to emphasize motions normal to the axis of elongation, and 
thus give incorrect results concerning the direction of drift. 

The reflection method of Mitra has been used to study drifts in many 
parts of the world. Most of the data have been analyzed only by simple 
procedures, such as the time-delay method. Deduced drift velocities are 
typically a few decameters per second, but the details of their magnitudes 
and directions are extremely complicated and consistent patterns do not 
always emerge. Survey articles, with references, have been given by 
Briggs and Spencer (1954), Shimazaki (1959), and Rawer (1963). 

The irregularities in the E and lower F region are often quoted 
as being a hundred meters or so in dimension. Probably there exists a 
distribution of excess electron density which contains a broad spectrum 
of spatial periods. When "correlation ellipses" are computed, the ratios 
of their major axes are generally found to lie in the range from 1.5 to 
2. The irregularities are thus not greatly elongated, and the directions 
of elongation do not seem to follow a very clear-cut pattern. The ratio 
V /v which, as previously stated, measures the relative importance of 
random change and bodily drift in producing fading is sometimes small 
and sometimes large, but its median value is near unity. Typical irregu- 
larities are thought to contain a fractional deviation of electron 
density, An/n, of a few percent. Fooks (1962) has quoted a figure of 
2 percent for irregularities in the E layer and lower F layer by day. 
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The irregularities principally responsible for radio-star scintil- 
lations are located in the F region, at heights of 250 km and above. 

On occasion they are above the level of the F2 peak and are therefore 
inaccessible to conventional ground-based sounding, although they can be 
studied by satellite techniques such as beacon transmissions and the 
topside sounder. 

Early observations of radio-star scintillations at 3.7 m wavelength 
(81 Mc/s) showed that records obtained from spaced receivers showed good 
correlation for spacings up to a few kilometers, but for spacings of 
200 km, the occurrences of scintillation at the two sites are related but 
the scintillations themselves are quite uncorrelated (Smith, Little and 
Lovell, 1950). This indicated a terrestrial (presumably ionospheric) 
origin of the scintillations. By using diffracting screen theory, Hewish 
(1952) was able to estimate an altitude of 400 km and dimensions of 
several kilometers for irregularities producing scintillations at 6.7 m 
wavelength (45 Mc/s). Steady drift velocities of 100 m/sec or more, 
were found. Later work on nighttime scintillation has tended to confirm 
Hewish 's conclusions in general, but some scintillation at low elevation 
angles has been associated with E g ionization (Wild and Roberts, 1956) 
and E region irregularities (Chivers and Greenhow, 1959). 

There seems to be general agreement that irregularities in the F 
region are considerably elongated in the direction of the geomagnetic 
field lines. Axial ratios as great as 10:1 have been reported, and under 
these circumstances it is difficult to deduce reliable north-south com- 
ponents of velocity. 

It is sometimes possible to make continuous observations of a cir- 
cumpolar radio star. The amount of scintillation may then depend on the 
elevation of the source above the horizon. At low elevations, scintil- 
lations become more severe because of the greater length of the slant 
ray path in the ionosphere and because irregularities are further from 
the observer. Since the ray path from a source at low elevation inter- 
sects the F region a thousand kilometers or so from the observing 
point, latitudinal variations of the severity and frequency of scintil- 
lation may be studied. 
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In order to make a statistical study of the occurrence of scintilla- 
tion, it is necessary to assign numerical indices to describe the amount 
of scintillation or the degree of spreading. This tends to be somewhat 
subjective, though it is possible for individual observers to achieve 
consistency in the analysis of a long series of records. Prolonged series 
of observations, extending over several years, have been made in Britain 
(Briggs, 1964; Chivers, I960) and in West Africa (Koster and Wright, 

I960). Scintillations occur mainly at night, but are sometimes observed 
by day. It is not clear whether the apparent seasonal variations observed 
in Britain are real, or result from the dependence of scintillation on 
the zenith angle of the radio source. But it is clear that scintillations 
are more common at sunspot maximum than at sunspot minimum, although this 
is the opposite of the solar-cycle dependence of spread F. It has also 
been found from satellite observations that the regions producing scin 
tillations are quite dependent on magnetic latitude (Yeh and Swenson, 

1959; Kent, 196l). A comprehensive survey of the effects of the iono- 
sphere on radio-wave propagation, including absorption, refraction and 
scintillation is given by Lawrence et al (1964). 

Large irregularities may be detected as distortions of ionogram 
traces or of high-frequency backscatter patterns. A particular type, 
known as "travelling ionospheric disturbances," produces "kinks" in the 
”o" and "x" traces near the F2 cusps on an ionogram, and these 
gradually move down toward lower frequencies (Munro and Heisler, 1956; 
Munro, 1957). These are believed to be wavefronts of distrubances , many 
hundreds of kilometers in extent, which may travel for horizontal dis- 
tances of 3000 km or more, with speeds of up to 10 km/min. The apparent 
vertical motion seen on ionograms can be accounted for if the wavefronts 
are supposed to be inclined to the vertical. Hines (i960) interprets 
these disturbances as gravity waves in the atmosphere. Because they 
travel so far and last for so long, it seems that their energy must be 
trapped in some kind of horizontal duct or "waveguide." If this ducting 
results from large scale wind circulations, as discussed by Hines (1963), 
then it may be possible to account for the preferred directions of motion 
noted by Munro. The energy is probably trapped at heights around 100 km, 
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but some upward leakage occurs and causes the observed effects in the 
F region. 

Spreading of F region traces on ionograms has been divided into 
two main types. One is "range spreading," in which multiple traces are 
seen at frequencies well below f^F2, and which seems to be prevalent at 
lower latitudes. The other is "frequency spreading," in which the high- 
frequency ends of the traces are branched or blurred. This is the type 
most commonly seen at higher latitudes, and the range of frequencies 
covered by the "spreading" gives an indication of the deviation of elec- 
tron density, AN, present in the irregularities. Spread F indices 
have been defined empirically and used for statistical studies, including 
comparisons with the empirical indices used for radio-star scintillations. 
Published bulletins of ionospheric data contain information about the 
occurrence of spread F, but scaling practices are known to vary and care 
must be exercised in the use of such data for detailed statistical studies. 
On account of their availability, they have been widely used for studies 
of worldwide occurrence of spread F. Among these may be mentioned the 
work of Shimazaki (1959) and Singleton (i960, 1962). 

From such studies it has been found that spread F occurs most fre- 
quently in equatorial and auroral latitudes, less frequently in temperate 
and polar regions. Spread F occurrence is positively correlated with 
magnetic disturbance at high latitudes but negatively at low latitudes. 

At low latitudes the phenomenon appears to be closely connected with 
the increase of height of the F layer after sunset. It has been sug- 
gested by Martyn (1959) that there is a causal connection between these 
phenomena. Singleton further finds that probability of occurrence is 
negatively correlated with F2 layer critical frequency. This appears 
to hold also for the solar-cycle variation at temperature latitudes. 

Thus spread F is less common at sunspot maximum than at sunspot minimum, 
opposite to the behavior of radio-star scintillation. This has been used 
as evidence against the view that the phenomena are closely related 
(Chivers, I960), but Briggs (1964) thinks that at sunspot maximum the 
greatest degree of irregularity occurs above the F2 peak, where it 
cannot produce effects on bottomside ionograms. 
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Ionograms obtained from the Alouette topside sounder have also given 
considerable information about the distribution of spread F. Calvert 
and Schmid (1964) found that it occurred most frequently at low latitudes, 
during the night and mid-morning, and at high latitudes, at any hour. 

It was found to be less common at middle latitudes. 

Much study has been devoted to the detailed form of spread F echoes 
on ionograms, both bottomside and topside, and the way in which they are 
related to the geometry of the irregularities. This naturally involves 
wave propagation calculations in an irregular medium. Considerable 
success has been achieved in relating field-aligned irregularities asso- 
ciated with equatorial spread F to the "ducts" which guide the radio 
waves transmitted by Alouette. In circumstances when multiple reflections 
occur, waveguide calculations have proved fruitful in explaining features 
observed on ionograms (Calvert and Cohen, 1961; Pitteway and Cohen, 1961; 
Muldrew, 1963), 

7 . Sporadic E 

Ionosondes often detect dense layers or patches of ionization in the 
E region, at heights of 100 km to 120 km, which do not seem to be related 
to the normal daytime E layer. This phenomenon is known as "sporadic E" 
or "Es" because it does not show any regular behavior. Sometimes Es 
appears in sheets which completely hide the overlying F layer. At other 
times it may be patchy and partially transparent to waves reflected from 
higher layers. 

In their comprehensive survey article, Thomas and Smith (1959) adopt 
the I GY definition of a sporadic E reflection as an E region echo 
characterized by one or more of the following: 

random time of occurrence, 

partial transparency (echoes also obtained from higher layers), 

variation of penetration frequency with transmitter power, 

virtual height independent of frequency. 

Much of the work carried out on Es has necessarily consisted of 
classification and tabulation of data. For this purpose, several types 
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of Es have been defined, which are mostly distinguished by the appear- 
ance of the Es traces on ionograms; sketches of these are shown in the 
Thomas and Smith review cited above, A more recent review, containing 
many original papers, has been published by Smith and Matsushita (1962). 

Sporadic E is studied with the aid of vertical incidence ionograms, 
oblique incidence ionograms and fixed frequency recordings. Backseat ter 
methods are also useful, particularly if a steerable antenna is used so 
that the motions of individual Es patches can be followed (Peterson, 

Egan and Pratt, 1959). Rockets have been flown through Es layers, and 
data obtained from propagation and probe techniques have indicated the 
presence of very thin layers of dense ionization at heights of 100 to 
130 km. There is some evidence, that these occur at a series of "pref erred " 
heights, at intervals of about 6 km (Pfister and Ulwick, 1958; Gringauz, 
1958; Smith, 1962). 

The parameters used to describe Es include the critical frequencies 
f Q Es and f Es . Es traces often do not show the usual "cusps" and 
sometimes a "top frequency" at which reflections are obtained, f Es, is 
used (Bibl, Busch, Rawer and Suchy, 1955). In many cases it may be 
assumed that f Es = f Es. The "blanketing frequency" f Es is the 
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lowest frequency at which echoes from higher layers are received through 
the Es layer. Generally f^Es < f^Es, so that the Es layer is 
partially transmitting over a (small) range of frequency. Sometimes the 
Es layer completely blankets all layers above it, and thus prevents 
observation of the F layer. Another parameter is the virtual height 
h'Es which is usually independent of frequency over most of the trace, 
in which case it should be nearly equal to the actual height of the base 
of the Es layer. 

Many surveys of Es occurrence do not distinguish between different 
types, but merely assess the frequency with which some arbitrary criterion, 
such as f^Es ^ 5 Me /s , is met. In temperate latitudes, occurrence of 
this condition has a probability of about 50 percent around noon in sum- 
mer, and there are lesser peaks of occurrence near noon in winter and 
midnight in summer. Diurnal minima of occurrence are found near sunrise 
and sunset, and seasonal minima at the equinoxes. In addition, variations 
with longitude may exist. 
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"Patches" of Es may be tens or hundreds of kilometers in extent, 
and drift with velocities of order 50 m/s. There seems to be little 
correlation between mid-latitude Es and magnetic activity. Some con- 
nection between Es and thunderstorms has been suggested (Rastogi, 1957), 
and certain E region echoes (not necessarily within the major Es 
categories) have been attributed to meteoric ionization (Naismith, 1954). 

It is difficult to see how sources of ionization could be sufficiently 
localized to produce Es, except perhaps for energetic charged particles, 
which seem more likely to be associated with high latitude Es, rather 
than the mid-latitude types. Any such source would have to operate for 
the observed lifetime of Es patches, which may persist for some hours, 
since the mean lifetime of ions in the E region is only a few minutes. 
Consequently, most attention has been given to explanations involving a 
redistribution of existing ionization rather than the production of extra 
ionization . 

It does not seem that ordinary turbulence could produce the strong 
nonuniformities of ionization present in Es. Dungey (1956) showed that 
a vertical gradient of the horizontal wind velocity (that is, a wind 
shear) could cause irregularities, and Whitehead (l96l) has shown that 
shears of east -west wind velocity can be very effective in producing 
thin layers of ionization. This theory seems promising, especially since 
rocket experiments reveal the existence of strong shears (see also Hines, 
1964). Moreover, the theory can explain the connection between Es 
occurrence and the strength of the horizontal geomagnetic field component 
(Heisler and Whitehead, I960). It has been shown that gravity waves in 
the neutral atmosphere may be the source of the wind shears necessary to 
this theory of Es (Axford, 1963), 

In high latitudes, one type of Es is associated with aurora (Knecht, 
1956). Other types are strongly associated with magnetic disturbance, 
and their geographical and temporal distribution shows interesting pat- 
terns (Thomas, 1962b) which suggest that direct particle bombardment is 
responsible for production of Es ionization. 

Near the magnetic equator, a distinctive type of "equatorial" Es 
is observed, which is patchy and transparent to waves reflected from 
higher layers. It is strongly associated with the "electrojet " current 
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(Sec. VI ) which flows along the magnetic equator by day (Skinner and 
Wright, 1957) and, like the electrojet, it is subject to some lunar con- 
trol (Matsuchita, 1957). The equatorial Es irregularities appear to 
be aligned with the geomagnetic field (Bowles, Cohen, Ochs and Balsley, 
I960), and it has been suggested that they are caused by plasma instabili- 
ties arising from the flow of the large electrojet current (Farley, 1963). 

8 . Artificial Disturbances 

The availability of rockets has led to several methods in which 
modifications to the naturally-occurring ionosphere are possible. A 
number of rocket flights have released various chemicals including sodium, 
aluminum oxide, and cesium at altitudes within the D and E regions. 

The objectives of a series of flights in 1960 (Project Firefly) were the 
study of the motions, luminescence, and decay of regions of high electron 
density (Gallagher, 1963) which describe these results, An interesting 
collection of ionograms relating to these experiments has been presented 
by Wright (1964), The wind velocities are found to range up to 70 m/s 
in the E region on many occasions and at times to exceed 150 m/s. 
Usually, there is also a variation of both speed and direction with 
altitude; the resulting wind shear may well lead to the sporadic E 
layers as we discussed in the preceding section. 

The radial growth of the artificial ion-electron clouds which are 
produced by these chemical releases permits the appropriate diffusion 
coefficients to be calculated. Zimmerman and Champion (1963) find the 
measured values for the diffusion of lithium, sodium or cesium through 
molecular nitrogen to be a factor of two to four smaller than theoretical 
values at altitudes between 100 and 150 km, 

The explosion of nuclear bombs at high altitudes also generates a 
variety of geophysical effects which have been used to study the upper 
atmosphere. The great release of energy will obviously produce greatly 
enhanced ionization over a considerable volume and a visual aurora may 
be observed at the magnetic conjugate point ow ? ing to the dumping of high 
energy electrons trapped in the earth’s field. Other electrons may 
"mirror” at higher altitudes and continue to precess about the earth in 


125 


SEL-64 -111 



"shells" for periods of months or even years. The synchrotron radio 
emission from these electrons has been measured and is the source of 
considerable concern to some radio astronomers (Peterson and Hower, 1963; 
Ochs et al , 1963). VLF and hydromagnetic waves are launched at the time 
of the explosion and the propagation of these waves may be studied all 
around the globe. These and other phenomena are described in the 
"Collected Papers on the Artificial Radiation Belt from the July 9, 1962, 
Nuclear Detonation" (1963), 
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V. RELEVANT ASPECTS OF GEOMAGNETISM 


1 . Introduction 

Any discussion of the ionosphere would be quite incomplete without 
consideration of the terrestrial magnetic field. In Sec. III-5 we found 
that in the transport regime” above the F2 peak, plasma diffusion was 
constrained to follow the earth’s magnetic field lines, and this fact was 
included in evaluation of the diffusion coefficient in the continuity 
equation. 

A number of the experimental techniques described earlier rely on- 
magnetic field effects. The polarization rotation produced by the Fara- 
day effect was measured in propagation experiments with satellites and 
moon echoes. The received spectrum from incoherent scatter was found to 
be greatly modified by the presence of a magnetic field. 

In addition to these reasons for studying the earth’s magnetic field, 
there are several phenomena which have not yet been considered. It 
appears that tides and thermal expansion of the neutral atmosphere gen- 
erate a sufficient emf by dynamo action to drive currents and produce 
additional polarization fields which can move the plasma in the ionosphere 
The height at which the electron-neutral collision frequency equals the 
gyromagnetic frequency determines the altitude of maximum conductivity, 
and thereby the currents which flow in response to the induced emf’s. 

We shall also describe ’’magnetic storms” which are usually accompanied 
by auroras and ionospheric deviations as large as ±50 percent from the 
normal electron density values. 

One of the earliest experiments related to geomagnetism was published 
in 1600 by Queen Elizabeth’s physician, William Gilbert, and it provided 
considerable insight into the nature of the earth’s magnetic field. As 
described by Chapman (l95l), Gilbert ”cut a spherical piece of the 
naturally magnetized mineral called lodestone, and examined the distri- 
bution of direction of the magnetic force over its surface by means of 
tiny magnetized needles freely pivotted. He saw that the distribution 
of dip agreed with what was known of the earth’s field. Hence, he con- 
cluded that the earth itself is a great magnet, similar to his magnetized 
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sphere except in size.” We can now show ma thematically that the external 
field of a uniformly magnetized sphere is identical to that of a dipole 
of appropriate magnitude placed at the center of the sphere. Since we 
find, as did Gilbert, that this does provide a reasonable first-order 
fit to the observed field near the earth’s surface, we will want to con- 
sider the dipole field in more detail a little later. 

By 1635 Gellibrand had discovered the secular variation in declination 
at London, and observations since then suggest a cyclic "wandering" of the 
magnetic declination and inclination with a period of about 500 years. 

The secular variation in the magnitude of the field has only been studied 
for a little over a century, since instruments to measure the absolute 
magnitude were not devised until the nineteenth century. 

Another major advance was made in 1722 by Graham, who was a clock- 
maker in London. He discovered the regular diurnal variations of the 
compass needle and also the larger perturbations which occur at times of 
magnetic storms. Further work by Celsius and Wilcke showed that magnetic 
disturbances are quite widespread and well correlated with the visible 
auroras. Surprisingly, it was still another century before Schwabe in 
1843 noted the periodicity of the sunspot number. It was soon established 
that this periodicity is reproduced in the frequency of magnetic distur- 
bance. 

In 1882 Balfour Stewart connected these magnetic variations with the 
region we now identify as the ionosphere. He argued that the diurnal 
variations were produced by currents flowing in the upper atmosphere in 
regions of appreciable conductivity. He also indicated that these currents 
could be driven by electric fields induced by motion of the air across the 
earth’s magnetic field lines and this is essentially the position main- 
tained today. The presence of the ionosphere was further implied by 
Marconi ! s high-frequency trans-Atlantic communication in 1901, but "direct" 
measurements were not made until about 1925 by Breit and Tuve and by 
Appleton, as we discussed earlier. 

The strength of the earth’s field is now measured at numerous mag- 
netic observatories located all over the globe. The standard components 
which are recorded on magnetograms are usually: (a) northward X; 
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eastward Y; downward Z; or (b) easterly declination D = arc tan (y/x); 

2 2 1/2 

horizontal intensity H = (X + Y ) 7 ; vertical intensity V = Jzj. 
Sometimes the dip angle or inclination, I = arc tan (z/h), is also given. 
Further details of a historical nature are given by Chapman in 
"Solar Plasma, Geomagnetism and Aurora" (l963). Probably the best dis- 
cussion of the entire area is found in the two-volume edition of Chapman 
and Bartels, "Geomagnetism" (1940), 

2 , The Geomagnetic Field 

In many cases, a very convenient way to describe a field is with a 
"potential function," Since the currents which produce the main field 
of the earth do not flow across the earth’s surface (measured atmosphere 
to ground currents average about 10 amps/m , which are too small to 
produce significant magnetic fields), the field is "curl-free" at the 
surface, V x B = 0, and it therefore may be obtained from a scalar 
potential, B = -W. Since V • B = 0 also, we see that the potential 
function V satisfies Laplace's equation. The solution to this equa- 
tion in spherical coordinates is well known to be a product of three 
functions, each containing only one of the variables, colatitude (e) 
increasing southward, longitude (<t>) increasing eastward and radial 
distance (r) increasing outward. It is written as 

. ^ °o n 

V = a zl (r) ^ P n m ( cos cos + h n sin (V-l) 

n-1 m=0 

In this equation P (cos 0) is the associated Legendre polynomial of 

n , m 

degree n and order m, and a is the earth radius. We have ignored 
a radial solution proportional to (r/a) n , which restricts our results 
to the region r > a and assumes that the currents producing the fields 
are inside the earth's surface. This assumption has been verified to 
within 0,1 percent, which is as close to zero as may be determined from 
presently available magnetic data. The index n = 0 is missing because 
this would correspond to an unrealistic magnetic monopole. The index m 
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increases only to n, since P (©) = 0 for m > n. We see that the 

n , m 

index n determines the order of the magnetic multipole, e.g., n = 1 
gives a dipole potential, n = 2 gives a quadrupole potential, etc. 

The orthogonality relation for associated Legendre polynomials gives 
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2 (n + m) 1 
( 2n + 1 ) (n - m ) l 


for n 1 ^ n 

(V-2) 

for n ! = n 


in which x s cos 9. We can see that when n = n’, the value of the 
integral is very strongly dependent on the choice of m. This would 
eventually be reflected in the values of g^ and h^ when these coef- 
ficients are evaluated by comparison with magnetic observations. To 
overcome this defect, Schmidt is credited with introducing a partially 
normalized Legendre polynomial (see Chapman and Bartels, 1940) defined 
by 


P m (x) = P (x) for m = 0 

n' ' n,m' 
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for m > 0 


(V-3) 


Equation (v-l) may be rewritten in terms of P™(x) instead of P r m (x), 

and the coefficients g m and h m will have different numerical values 

n n 

owing to the normalization. It should also be noted that these poly- 
nomials can be obtained from the ordinary Legendre polynomials P n (x), 
since 
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(V-4) 


When we consider a single magnetic dipole at the earth’s center, 
oriented along the line of the geographic poles, these equations become 
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very simple. In this case n = 1, m = 0 and P°(x) = P^ Q (x) = P^x) 
cos e. Thus 

V = 3 g l (r) COS 9 (V-5) 


Since B = -W = -l^bv/br) - i Q (l/r)(dv/de) - l^l/r sin e)(dv/d*), we 
may obtain the magnetic field components by differentiation. Usually 
the values of Z = (dv/dr), radially downward, X = ( l/r) (dv/d0 ) , north- 
ward, Y = -(l/r sin e)(dv/d<t>), eastward are those obtained. For the 
potential (V-5) we obtain 


v 0 /a\ 3 . „ 0 

X = " g l (r) Sin 0 = X 1 

Y = 0 (V-6) 

Z = " 2g l (f) COS 6 = Z i 


and 


B 



+ Y 


♦ z 2 ) lh 



+ 3 cos 


2 



Another useful parameter is the magnetic dip angle I, defined by 


tan I = 


(x 2 + Y 2 ) 1 / 2 


= 2 cot 0 


(V-7) 


for the simple dipole field. The value of g^ may be identified with 
the southward component of the field strength (more precisely, the mag- 
netic induction) at the equator and, for the earth, g^ = -0.305 gauss = 
-3.05 x 10 5 webers/m 2 . 
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It is impossible to fit the earth's field very well with a single 
dipole term as we have considered above. As the next better approximation, 
we can let n = 1 and m = 1, which generate two additional terms in 
the potential expression. We find P*(x) = X (x) = sin 9, and 

2 

V = a ^ jg° cos 6 + ^g* cos <t> + h* sin <t^ sin ej (V-8) 


The first term in the bracket corresponds to (V-5) while the second term 
provides two new contributions to each of X, Y and Z, Physically, 
these may be considered to be generated by two additional orthogonal 
dipoles placed at the center of the earth but with their axes in the 
plane of the equator. From (-W), we obtain 


3 

xj = cos 0 cos + h* sin 


= (l) («i 


sin - h cos 


(V-9) 


- -(|) 2 sin 9 (gj 


cos 4> + h^ sin 


)J 


The terms in Eqs. (v-9) and (V-6) must be added as appropriate to obtain 
the fields produced by potential (V-8). The effect of g* and h* is 
to incline the total dipole term to the geographic pole by an amount 



(V-10) 


Thus, the potential (V-8) could be produced by a single dipole inclined 
at an angle 0! to the geographic pole and with an equatorial field 
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strength 


g l/ + ( g l) + ( h l) J ' For the eart *i, the single dipole 

equivalent would be tipped by about 11.5° toward 70° W longitude and 

have an equatorial field strength of about 0.312 gauss. 

For n > 1, additional multipoles contribute to V. Finch and Leaton 
(1957) have computed the first 48 terms in Eq . (v-l) (through m = 6, 
n_=_6) for epoch 1955 and their first eight terms are given in Table IV 
below. The dipole term clearly predominates and the higher order multi- 
poles are not only smaller in magnitude, but they also diminish more 
rapidly with altitude. However, even this model cannot give any repre- 
sentation of the local variations of the geomagnetic field, as only 
features whose scale is of order one-sixth of the earth's radius can be 
accurately described by formulas involving only m ^ n ^ 6. 


TABLE IV. FINCH AND LEATON COEFFICIENTS 
(values in gauss) 


0 




g l 

— 

0.3055 


1 

g l 

= ^ 

0.0227 

= + 0.0590 

0 

g 2 

= - 

0.0152 


1 

g 2 

= + 

0.0303 

h* = - 0.0190 

2 

g 2 

= + 

0.0158 

h o= + 0.0024 


A significant improvement can be obtained, while retaining a simple 
mathematical form, by using "eccentric dipole coordinates" (Cole, 1963). 
In this model, the dipole is moved away from the earth center (about 
400 km toward the western Pacific) to obtain a better world-wide fit to 
the magnetic field. 

3. Regular Variations 

We have already noted that the magnetic elements (X, Y, z) undergo 
a regular diurnal oscillation about their mean values. Figure 25 
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shows the solar daily magnetic variation, or Sq, for quiet days in the 
1922-1933 equinoctial periods at various latitudes (Vestine et al , 1947) 

The more obvious features include the phase reversal in AY and AZ 
(but not AX) across the equator and the large increase in AX near 
the equator. The latter feature is evidence of the equatorial electro- 
jet, previously discussed in connection with ionospheric conductivities. 


Local mean hours 



FIG. 25. THE Sq VARIATION IN AX, AY AND AZ AT VARIOUS 
LATITUDES (after Vestine, I960). 
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It is again convenient to find a scalar magnetic potential from 
which these variations may be calculated. The starting point is usually 
with data such as shown in Fig. 25, from which a Fourier series is cal- 
culated to represent the diurnal variation at each station. Then it is 
assumed that the same series is applicable at all points along the cir- 
cle of latitude through the station under consideration. This permits 
the variable in the Fourier series (time) to be related to the variable 
in the spherical harmonics (longitude). The proce dur e is unfortunately 
rather complicated, but it is explained very carefully in Chapter XX of 
Chapman and Bartels (1940). The analysis is further complicated by the 
necessity of separating the potential into internal and external parts. 
Chapman found the external part of the potential to exceed the internal 
part by about 2.5 to 1.0. The latter is presumably produced by currents 
within the earth induced by currents in the ionosphere. 

Another very convenient way to represent the source of the regular 
magnetic variations is by a current sheet, flowing on the surface of a 
sphere above the earth. This surface is usually assumed to be at an 
altitude of about 100 km where the maximum conductivity is found. The 
current density can be expressed directly in terms of the magnetic 
potential (or vice versa), as shown by Chapman and Bartels. The current 
systems which they have found required to produce the external part of 
the Sq variation for sunspot minimum, equinoctial conditions, is shown 
in Fig. 26. Current flows parellel to the streamlines with 10,000 amperes 
between each line. 

The main features of Fig. 25 can be related to the current system of 
Fig. 26 by inspection. We remember that an eastward sheet current above 
the ground will produce a northward magnetic field at the surface, etc. 
Therefore, between the latitudes of about ±40 deg, the eastward daytime 
current produces a positive AX component, as can be seen at Alibag and 
Honolulu. At higher latitudes, however, the phase reverses to a nega- 
tive AX (southward) by day. The transition occurs near the current 
focus at the latitudes of Tucson and Watheroo. The current sheet does 
not delude the equatorial electrojet required to explain the large 
variation at Huancayo. A narrow current belt along the magnetic dip 
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FIG, 26, THE OVERHEAD CURRENT SYSTEM CORRESPONDING TO 
THE EXTERNAL PART OF Sq FOR SUNSPOT MINIMUM, EQUI- 
NOCTIAL CONDITIONS (after Chapman and Bartels, 1949), 


equator could be added to produce this feature. The AY variation 
implied by Fig. 26 should be eastward prior to about 1100 local time and 
westward in the afternoon for northern latitudes and reversed in the 
southern hemisphere and this is also observed in Fig. 25. The AZ 
variations are not understood quite so easily, since they must be pro- 
duced by currents removed from the point in question since overhead 
currents produce horizontal variations only. The vertical variations 
are most pronounced near the current focus because the counter-clockwise 
current flow produces an upward magnetic field near the focus in the 
northern hemisphere while a positive AZ is produced in the southern 
hemisphere near Watheroo, At night, the current densities are much 
smaller owing to the reduced conductivity in the E region, 

A similar magnetic variation can be found with the fundamental period 
of a lunar day, about 24.8 hours. The magnitude of the lunar (l) varia- 
tion is much smaller than the solar (Sq) variation, being on the order 
of ±4 7 . Also, the predominant mode is a semidiurnal variation as is 
found in ocean tides, rather than the predominant diurnal period in Sq. 
This is particularly important since the L variation is presumably 
induced by tides alone while the Sq variation may be produced by both 
solar gravitational tides and by thermal expansion. 
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4. 


Storm Variations 


In addition to the regular Sq and L variations discussed above, 
uhere are sudden and sometimes rather large disturbances of the magnetic 
field which are observed all over the world. As these transient ampli- 
tude fluctuations grow larger, the "magnetic activity" is said to increase 
and the period may then be identified as a "magnetic storm. " 

The degree of magnetic activity is determined by an examination of 
the magnetograms obtained from the various observatories. For each 
station individually, according to its own average activity, a K index 
scale is assigned which varies between 0 and 9 . The index is not linear 
but varies quasi -logarithmically wdth magnetic disturbance. The quali- 
fication for K = 9 might be a disturbance of about 3007 for a low 
latitude (but not equatorial) station, 5007 for a mid-latitude station 

and 20007 for a station in the auroral zone. We remember that 1 gauss 

5 -4 2 

is equivalent to 10 7 or 10 webers/m . The K indices from a set of 
observatories are then combined - to yield a "planetary" figure, Kp, for 
every 3-hour interval of each day. The Kp scale is determined to an 
accuracy of l/3 unit, running as: Oo, 0+, lo, 1+, 2-, 8 +, 9-, 9o. 

The maximum value 9o has been recorded only on a few occasions. 

The range of variation in gammas for each of the elements (d, H, v) 

or (X, Y, z) at each station may also be scaled from the magnetograms . 

The greatest range in each three-hour interval is called (a) and it is 

from (a) that the K index was assigned. The arithmetic sum of the 

(a f s) at 12 representative stations yields a planetary value (a ) and 

P 

its daily average is published as A . The planetary figures Kp and 
^p are computed in De Bilt, Holland and are published in several places, 
such as the National Bureau of Standards Publication, CRPL-f(b) series. 

The days on which magnetic storms occur may follow the observation 
of a solar flare on the sun or at times of high solar activity. These 
flares or emitting regions are believed to eject streams of charged 
particles with velocities on the order of 300 to 1000 km/sec which then 
require two to five days to reach the earth. The interaction of the 
charged particles with the earth’s atmosphere and magnetosphere then 
produces the magnetic storm, but the manner in which this is done is 
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still an area of active debate and research. The arrival of the solar 

plasma is many times revealed by a "sudden commencement” observed nearly 

simultaneously all over the world. 

The disturbed value of a magnetic element (we will use X as an 

example, but any of the other elements may be treated similarly) is 

defined by AX = X , - X , where X . is the observed value and X 

J D obs q obs q 

is the quiet day value, including the regular Sq and L variations. 

The values of AX^ depend not only on "storm time" T, which is the 
number of hours elapsed after the commencement of the storm, but also 
on the colatitude 0 and longitude A of the observing station. The 
disturbed value is usually split into two parts: (a) the Dst or 

storm time variation, which is the average value of AXp around a circle 
of constant latitude and (b) the SD solar daily disturbance or the DS 
disturbance local-time inequality . All of these three quantities vary 
with storm time, and are defined by 

= Dst + DS = Dst + SD (V-11) 


The SD description was first developed when it was thought that 
the longitude variation of AXp maintained a fixed orientation in solar 
time. It is now known that the phase and amplitude of the longitude 
(or local time) variation changes during the course of a storm, and the 
quantity DS has therefore been defined (Akasofu and Chapman, 1964). 
These authors now consider that "SD is merely the mean of DS over 
a day," but it would appear that SD could be the average of DS over 
any convenient interval. These relations have been expressed as 


where 


X_ = a + > c sin (nA + £ ) 

D o n , n' 

n 


(V-12) 


Dst = a 

o 


DS = 


I 


c sin (nA + £ ) 
n v n' 
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The evaluation of these coefficients at any latitude requires either the 
observation of many storms of comparable intensity at a single station 
or the observations of several stations at the same latitude for a single 
storm. Ideally, both of these approaches will lead to the same result, 
but some differences are to be expected in practice owing to local 
variations, different storm characteristics, etc. 

These quantities have been defined with a particular model of a 
magnetic storm in mind, although there is no assurance that this model 
is correct in any detail. The Dst variation, which is the same all 
around a circle of latitude, must be produced by a current loop encircling 
the globe. Observations suggest that the main phase is produced by such 
a ring current flowing from east to west and perhaps at an altitude of 
several earth radii (Akasofu and Chapman, 1961 ). The DS variation 
would be produced by currents flowing in the lower ionosphere, particu- 
larly intense near the auroral zone and presumably related to the influx 
of particles in these regions. 
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VI. DYNAMO THEORY 


1 . Atmospheric Oscillations 

The dynamo theory was first suggested by Stewart (1882) to account 
for the daily variations of the geomagnetic field and developed quanti- 
tatively by Schuster (1908). Although there are serious difficulties 
with some aspects of it, the theory may be partially correct and we 
shall review it here, A good account has been' given by Hines (1963). 
The basic stages in the discussion are as follows: 

(1) The sun and moon produce tidal forces in the atmosphere, the 
periods being fractions of the solar day (24 hr) and lunar 
day (24.8 hr ) . 

( 2 ) These forces set up standing waves in the atmosphere below 
100 km, which result in (primarily horizontal) air motions. 

( 3 ) There may be a natural resonance of the atmosphere of about 
12 hr period, which selectively amplifies the solar semi- 
diurnal component. 

(4) The motion of the air across the geomagnetic field induces 
electromotive forces, which drive currents at levels where 
the electrical conductivity is appreciable (principally in 
the E region), thus causing the solar quiet-day and lunar 
magnetic variations. 

(5) Because of the vertical and horizontal variations of conduc- 
tivity, currents cannot flow freely in all directions; and 
polarization charges are thereby set up, modifying the flow 
of current, 

(6) The electrostatic fields associated with these charges are 
transmitted to the F region, via the highly-conducting geo- 
magnetic field lines, where they cause electromagnetic drifts. 

Stages (l) to (3) come within the field of atmospheric dynamics and 
may be thought of as the ’'mechanical 11 side of the dynamo theory, which 
is shown on the left of Fig. 27. The ’’electrodynamic" side comprises 
stages (4) to (6). 

Let us now review the mechanical side in greater detail. First, 
we note that the 12-hr oscillation is observed to be the dominant com- 
ponent of the pressure oscillations at ground level. At low latitudes, 
the 12-hr oscillation is about 1 mb in amplitude and is visible on 
ordinary barograms, though in middle latitudes it is concealed by the 
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MECHANICAL 


ELECTROOYNAMIC 


3 l 4-62 



FIG. 27. BLOCK DIAGRAM OF THE DYNAMO THEORY, SHOWING 

(LEFT) ITS MECHANICAL AND (RIGHT) ITS ELECTRODYNAMIC 

ASPECTS. The predominant period of the variation is 

indicated in parenthesis, but other components may 

also be of appreciable size. Only a small fraction 

of the total current J is required to establish 

and maintain the electrostatic field, E . 

— s 

irregular meteorological variations. Some typical amplitudes of pressure 
oscillations are as follows: 


Solar diurnal 
Solar semidiurnal 
Lunar semidiurnal 


Tropics Mid-latitudes 


o 

00 

mb 

0,3 

mb 

1.3 

mb 

0.4 

mb 

0.08 

mb 

0.03 

mb 


For gravitational tides which have a fundamental semidiurnal period, 
the lunar driving force (period 12.4 hr) exceeds the solar force (period 
12 hr). This is consistent with the observed domination of ocean tides 
by the moon, and we might expect the same situation to apply in the 
atmosphere. However, there is also a solar driving force arising from 
the periodic input of thermal energy. This has a fundamental 24-hr period 
but its harmonics (l2-hr, 8-hr, 6-hr, etc.) might have appreciable ampli- 
tude. Even so, it is not obvious why the 12-hr period should predominate. 

To overcome this difficulty, it was suggested that the atmosphere has 
a natural period of oscillation close to 12 hr, which amplifies the solar 
semidiurnal oscillation but is too selective to respond to the lunar 
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12. 4 -hr oscillation (Thomson, 1882). Although this coincidence is pos- 
sible, it might seem too good to be true since resonance is rare in 
nature. But before we can discuss it further, we must consider how the 
atmosphere responds to these periodic driving forces. The energy, both 
thermal and gravitational, is introduced primarily into the lower atmo- 
sphere, where the density is greatest. The pressure oscillations are 
propagated upwards, and increase in amplitude because of the upward 

decrease of density. Provided no energy is dissipated, the energy density 

2 

(proportional to pW where W is the amplitude of the oscillatory wind 

_l/2 

velocity) is independent of height; thus, W cc p . Now p decreases by 

6 ' 
a factor of about 10 between the surface and 100 km, so there is a con- 
siderable upward increase of W, known as "tidal amplification,” A 
consequence of this amplification is that the fractional oscillations 
of pressure and density, Ap/p and Ap/p, increase upward. At the 
ground Ap/p ~ 10~ but at 100 km, this ratio would approach unity, 
but for the limitations set by other processes. 

The upward progress of the tidal wave is interrupted by partial 
reflection, which occurs at heights where the temperature gradient is 
negative. This reflection is expected to be so pronounced in the meso- 
sphere at about 60 to 80 km that a system of standing waves is established 
and most of the oscillatory energy is trapped below the mesopause. There 
is, however, some upward leakage of energy which is dissipated at greater 
heights through various processes. One reason for the dissipation is the 
tidal amplification discussed above, because the simple linear wave theory 
holds only where Ap/p « 1. The tidal wave then "breaks" and its energy 
is degraded into smaller scale motion, ultimately being lost in viscous 
and electrical resistive damping. Figure 28 shows the amplitude and 
phase of Ap/p for the solar semidiurnal oscillation, as a function of 
height (Weeks and Weekes, 1947). There is a node at about 30 km, and 
the phase gradient above 80 km is connected with the leakage of energy 
from the top of the system. 

The trapping of tidal energy below the mesopause gives rise to the 
possibility of resonant oscillations, which we mentioned earlier. The 
work in this field has been reviewed by Wilkes (1949). The resonant 
periods are very sensitive to the temperature profile and so should 
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FIG. 28. AMPLITUDE (SOLID LINE) AND PHASE ANGLE 
(BROKEN LINE) OF THE SEMIDIURNAL PRESSURE 
VARIATION, AS A FUNCTION OF HEIGHT, AFTER 
WEEKES AND WILKES (l947). The dotted curve 
shows a growth amplitude given by exp (z/ 2 ). 

The departure of the solid line from this 
curve near 30 km represents a node in the 
standing wave pattern, produced by internal 
reflections from the ground and the mesosphere. 

vary with position and time. In addition, the reflection of energy is 
affected by winds in the stratosphere and mesosphere (Charney and Drazin, 
1961); thus the calculations of resonant periods are very uncertain. 

There may well be a broad resonance at a period of about ten hours, which 
would have some effect in increasing the ratio of the 12-hr to the 24-hr 
oscillation, but the theory that a sharp resonance exists near the 12-hr 
period now appear untenable. Small and Butler (1961, 1963) have computed 
the response of the atmosphere to a thermal tide arising from absorption 
of solar ultraviolet radiation in the ozonosphere, and find the semi- 
diurnal component to exceed the diurnal. 

As with ocean tides, the tidal velocities in the atmosphere are 
principally horizontal. At the ground, the tidal wind velocities are 
only of order 0.05 m/s and are thus very difficult to detect in the 
presence of much larger meteorological winds, but at 100 km the velocities 
have amplitudes of about 50 m/s, as measured by drifts of meteor trails. 
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2. Conductivities 


Before we can discuss the flow of current resulting from dynamo action, 
we must deal with the electrical conductivity of the ionosphere, which is 
anisotropic because of the geomagnetic field. We shall derive the results 
obtained by Baker and Martyn (1952, 1953), Chapman (1956) and others. 

We will start by investigating the motions of the charged particles, 
and use the following symbols: mass m, charge e(±), electronic charge 

-e, angular gyrof requency go = Be/m, collision frequency with neutral 
particles v. Suffixes i and e are used for positive ions and elec- 
trons. We ignore negative ions, since they probably do not play an 
important role at the heights where conductivities are large, at any rate 
by day. 

Consider first a particle subject only to collisions and a force F. 

We are primarily concerned with electrical forces for which F = eE. If 
each particle has a constant probability of collision, v per unit time, 
and has (on average) zero velocity after each collision, we may write 
the drift velocity as v = kF, in which k is the mobility of the parti- 
cle. The effect of collisions is then represented as a constant retarding 
force mvv. Then the equation of motion in electric and magnetic fields 
is 

m dv/dt = eE + ev X B - mvv (VI-l) 

If e j| B, the motion is unaffected by the presence of the magnetic 
field. For the more interesting case of E 1 B, the motions are as indi- 
cated in Fig. 29, where for convenience, we have omitted the last term in 
(VI-l) and assumed that the particle is brought to rest at equal intervals, 
l/v. The trajectories are then cycloidal, and are shown for different 
values of v/co. Both ions and electrons move in the direction E X B 
if v/cjd = 0, and in general, their motion is inclined at an angle 
[arc tan (co/v)] to E. 

To obtain the actual mobilities, we find the average drift velocities, 
over periods long compared with both l/v and l/co, by setting dv/dt = 0 
in (VI-l). Hence, in a cartesian coordinate system in which B is 
directed along the z axis, 
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FIG. 29. IDEALIZED TRAJECTORIES FOR ELECTRONS AND IONS SUBJECT 
TO AN ELECTRIC FIELD IN' THE PLANE OF THE DIAGRAM AND A MAG- 
NETIC FIELD DIRECTLY OUT OF THE PLANE. Charged particles are 
assumed to collide with neutral particles at regular inter- 
vals l/v, and to possess zero velocity after each collision. 
In order to show both electronic and ionic gyrations, the 
diagrams are drawn as though a) e /a> ± = 10 (instead of ^ 10 4 
as in reality). All the trajectories refer to equal inter- 
vals of time, namely 5 ionic (or 50 electronic) gyroperiods. 
Numbers in brackets refer to approximate heights, in kilo- 
meters, at which the conditions occur. In general, the head 

of the vector representing v, or v lies on a semicircle. 

— l — e 

0 = eE + ev B - mw 
x y x 


0 = eE -evB-mvv 
y x y 


( VI -2 ) 


0 = eE - mvv 
z z 


By solving these equations for v , v , v we find three different com- 

x y z 

ponents of the mobility. The solution can be compactly written in terms 
of a tensor mobility k, which satisfies the equation v = ek • E. It is 
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where the suffixes "0," "t," "h" denote "direct" (or "longitudinal"), 
"transverse," and "Hall." Note that the charge € is not included in 
K, so that according to our definition K is the ratio of velocity to 
force. This is applicable whether the force is electrical or mechanical, 
and is found more convenient here than a quantity defined as a ratio of 
velocity to electric field . 

If we define oo = Be/m (always positive) we can write each component 

in two forms, one of which includes B explicitly. Because the Hall 

current flows in the direction of B X E in the ionosphere, whereas the 

ions and electrons actually drift in the direction of E X B, the quantity 

k in the preceding equation is negative for ions and positive for elec- 
h 

trons . 


, u „ v 1 1 CD 

"Direct" (||B, ||e) , k 0 = mv “ B e v 

. ,i \ , 1 v 2 L vcu 

"Transverse" (-UB, ||E) , k t = mv 2 2 “ Be 2 2 

V + CD V + CD 


"Hall" (- for 
ion, + for 
electrons ) 

At any given height in the ionosphere, oo/ v is not the same for ions 
as for electrons. The height variation of the mobilities for one kind of 
carrier is shown in Fig. 30, in which altitude is measured from the level 
where v = co and the logarithmic mobility scale is referred to the value 

l/Be. 


(iB, lE),±k h = jjk-2 


CDV 


CD 


+ CD 


Be 2 2 

V + CD 
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FIG. 30. MOBILITIES FOR A SINGLE CHARGED SPECIES, 
AS A FUNCTION OF REDUCED HEIGHT z. It is 
assumed that the collision frequency v oc e" z 
and that the gyrof requency a) is height inde- 
pendent, The level z= 0 is taken where 
v = ch. The longitudinal, transverse and Hall 
mobilities (kg, k^, ^h) are plotted in units 
of l/Be. The combination of two such diagrams, 
one each for electrons and ions, placed at their 
appropriate heights, gives the conductivity per 
ion pair (Fig. 31 ) . 


To find the electrical direct -current conductivity cr, we combine 
the ionic and electronic motions, using the equations 

i = or • E = Ne(v. - v^) (VI-4) 


where is the current density and N the electron and ion concentra- 

tion. It is customary to define four components of conductivity as, 


"Longitudinal” or "Direct": cr = Ne (k + k ) 

0 ' Oe Oi ' 


2 

"Transverse" or "Pedersen": cr = Ne (k + k ) 

1 v te ti 


*2 = Ne Ke + k hi^ 
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"Hall"; 


"Cowling": 


= (4 * 4)^1 

These conductivities are all positive quantities. The electron and ion 

contributions to cr^ are of opposite sign, but as we shall see, 

[ k I 1 Ik I at all heights in the ionosphere, 
he hi 

Three of these appear in the tensor expression for a, which follows 
from these definitions and the expression for k: 

(VI -5) 


whereas the Cowling conductivity appears at a later stage in the theory. 

The vertical variations of these quantities are sketched in Fig. 31, in 
which the contributions of the ions and electrons are shown in units of 
Ne/B. The resulting sums are thus conductivities per ion pair , and should 
be multiplied by the electron density profile to obtain the actual vertical 
variation of each cr. 

The important levels in the diagram are defined by (about 

80 km in the actual ionosphere), cd.cjd = v. v (about 105 km) and cd. = v, 

r ' x e i e v ii 

(about 140 km). We assume ( v /a> ) = 1000 ( v e / CD e )» this ratio being 

relatively independent of height. Peaks of ct^/n occur at the levels 

03 = v and o> = v , but the former does not give an actual peak of cr_ 

e e l i i 

because of the very low value of N at 80 km. Numerical values of the 
conductivities have been given by Chapman (1956). 

Baker and Martyn (1953) have discussed the "effective” or "layer" 
conductivities in the actual ionosphere. These arise because of the 
limited vertical extent of the conducting layer in the E region. If 
an electric field is generated in the ionosphere by dynamo action, the 
resulting current cr * E may contain a vertical component. If so, charge 
will accumulate at the boundaries of the conducting layer because the 
current cannot flow into the region of low conductivity. These "polari- 
zation" charges will modify the electric field E until the resulting 
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Ne/B 





FIG. 31. CONDUCTIVITIES PER ION PAIR (a Q , Uj , a 2 , cr 3 ) PLOTTED 
ON A LOGARITHMIC SCALE RELATIVE TO THE VALUE Ne/B, AS A 
FUNCTION OF REDUCED HEIGHT z FOR AN IDEALIZED ISOTHERMAL 
MODEL ATMOSPHERE. It is assumed that collision frequency 
v cc e” z and gyrof requency a) is independent of z, for both 
positive ions and electrons; also that Vj_o) e / v e a^ £ 1000 
(independent of height), this being a reasonable value for 
the actual ionosphere. The zero of z is the level where 
v i v e = ^i^e' ancJ the im Portant levels where and 

v e ~ U) e are therefore situated at z = ±l/2 ln(l000); also 
shown are the approximate altitudes at which they occur in 
the actual ionosphere. 


flow is horizontal. Therefore, when we impose the condition of zero 
vertical current, we can eliminate the vertical electric field from the 
equations, and we find that the 3x3 tensor a can be replaced by a 
2x2 tensor a 1 , representing the "layer conductivity” whose components 


*■ -» ■’ 



- IMSC l 



SEL-64-111 


Hill KJII.Ii !. I :i 111 m : lil :|l!il KlipliH I f|i! « 


depend on the magnetic dip angle I. From this point on, we use coor- 
dinates x, y for the magnetic southward and eastward directions, and 
denote the dip angle by I. Then we can write the layer conductivity as 


cr’ - 


a cr 

xx xy 

-<J cr 

xy yy 


(VI-6) 


where 


°o a i 


xx . 2 2 2 

0q sin I + o*^ cos I sin I 


cr 


°0 a 2 Sin 1 


xy . 2 T 2 _ sin I 

ct 0 sxn I + cos I 


(VI-7) 


cr 


yy 


cr 2 cos 2 I 

2 2 + °1 
cr^ sin I + cr^ cos I 


The approximations given at the right arise because generally cr^ » cr 
or a , but they are not valid near the magnetic equator where 1=0, 

A 

and where we have 


tr = cr n , cr = 0, cr = (cr 2 + cr 2 )/cr = cr (VI-8) 

xx 0’ xy yy V 1 2/' 1 3 

Two consequences follow from these latter results. First, the very 
high conductivity cr = <t along the field lines insures that these 
lines are approximately electric equipotentials . Second, the east-west 
conductivity u at the equator is extremely large, being the "Cowling" 

yy 

conductivity cr^ which is comparable to This highly conducting 

strip along the magnetic equator carries a large current known as the 

"equatorial electrojet," which is confined to the region a few degrees in 

2 2 

width, where cr^ sin I « cr^ cos I. 
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The simple model of the dynamo region is therefore a relatively thin, 
horizontally stratified layer. The effects of the currents as observed 
at the ground can most conveniently be calculated with the use of "inte- 
grated layer conductivities," 2^ - /a^dh and = /<y dh, where the 
integration is made vertically through the conducting region. Since the 

major part of the conductivity is found in the E region, 2^ and 2 

•§T.®.^ a PP rox i ma tely proportional to N^E. At night, they become small 
because the E layer almost disappears. 

We should also consider the effects of a force due to a wind of veloc- 
ity W in the neutral air, as Martyn (l953z) and Weekes (1957) have done, 
Wg may take the force exerted as F = Wmv, as in Eq. (VI-l). If W is 

parallel to the magnetic field, the plasma moves with velocity v = kF = 

( 1 /m * ) W = W. However, all the other components of the tensor mobility 
k ...depend on ( //a>) or ( v'/ao) and thus decrease rapidly upwards. At 
heights above the E layer a wind is ineffective in moving ionization 
across the field lines, whereas an electric field E is much more effec- 
tive since it causes the ionization to drift with the velocity E X b/b^. 
In the D region, where v » cjd for ions and electrons, winds rather 
than electric fields are effective in moving ionization. In the E 
region, both winds and fields produce motion of ions and electrons in 
different directions, therebj^ generating electric currents. 


Motions in the Ionosphere and Magnetosphere 

The discussion up to this point may be summarized w ? ith the two 
equations 


J = Z 1 


(VI-9) 


E t = E ± + E _ = W x B - 


(VI-10) 


The horizontal sheet current J f lows in response to the total electric 
field E^ , which is composed of an induced part (w X b) and an electro- 
static part, -V<t , which was established in order to require the current 
flow to be horizontal. Obs ervations of the ma gnetic field strength at 
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the ground permit the current flow to be calculated from Ampere's circuital 
relations. If horizontal variations can be neglected, the southward (x) 
and eastward (y) components of magnetic field H (or induction B) and 
layer current density J are related by 


- - J = H = B/p 
2 y x x' "o 



X 


= H 


y 



(VI-ll) 


The numerical factor in these equations should be about 3/4, instead 
of l/2, because of the contribution of currents induced m the ground by 
the overhead current systems. Detailed analysis shows that ground cur- 
rents cause about one-third of the observed magnetic effect, though the 
proportion is not the same for the different magnetic components (Chapman 
and Bartels, 1940). If the ground were a perfect conductor, and compli- 
cations due to horizontal variations and the earth’s curvature were 
neglected, the ground and ionospheric currents would be equal and oppo- 
site, and would contribute equally to the observed magnetic variations. 

The coefficients in Eqs. (VI-ll) would then be 1. 

If there are electrostatic fields in the E region, they must exist 
also in the F region owing to the very high conductivity parallel to 
B. Hence, the E and F regions are electrically "coupled," as sug- 
gested by Martyn (1947) and investigated by Farley (1959) and others. It 
is believed that E region fields are accurately "mapped" into the F 
region, even for scale sizes down to a few kilometers. Above the E 
region, only the electrostatic component of the electric field is present, 
because the tidal winds are thought to be damped out by viscosity. Hence, 
the electric field acts on the ionization in the F region, causing it 
to drift with the velocity 


E X B 
— s — 



(VI- 12 ) 
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This has been likened to an electric "motor" running without load, as 
the "back emf" v x B set up by the drift of ionization across the 
field lines, just balances the applied field E g . Actually, a small cur- 
rent flows, because the electron and ion drift velocities are not pre- 
cisely equal, implying that some energy is expended in collisions and in 
overcoming gravity. The vertical component of the F region drift veloc- 
ity is produced by the eastward electric field component E and is given 

by v = (E /b) cos I. 
z y ; 

Dynamo theory, as condensed into Eqs. (VI-9) to (vi-12), is perhaps 

deceptively simple. In practice, there are numerous difficulties in the 

application of this theory, such as the interaction of the drifting plasma 

with the neutral atmosphere (Cowling, 1945; Hirono and Kitamura, 1956; 

Dougherty, 196l), which can reduce the vertical velocity v of the 

z 

plasma . 

We now wish to extend our discussion to higher heights, and introduce 
some concepts used in some modern theories of magnetic storms, which 
provide an alternative way of looking at the motions of ionizations. 

Since the geomagnetic field lines above the E region are (almost) 
electric equipotent ials , the electrostatic fields which cause F layer 
"motor" drifts must also permeate the magnetosphere. At these levels 
the fields. B and E and the ionization drift velocity y satisfy the 
so-called "f reezing-in " equations 

O = E + v X B 

( VI-13) 

v = E X b/b 2 


The second of these has already appeared as the F region motor equation 
(VI-12). 

The "freezing-in theorem" is of very general application to the 
magnetosphere, stellar atmospheres, and ionized interstellar matter. It 
is valid whether or not the electric field E is purely electrostatic 
(i.e., derivable from a potential $), providing its time variation is 
slow enough. Basically, the theorem states that the ionized matter and 
magnetic field moves in such a way that 
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(i) a set of charged particles, if connected by a field line, 
remain so connected at all times, 

(ii) any given tube of ionization bounded by a set of field lines, 
always encloses the same magnetic flux. 

The proof of these statements is somewhat involved, and has been given 
by Dungey (1958) and by others. Since a "line of force” is an abstraction 
and not a physical object, its ”motion” is somewhat arbitrary. In ele- 
mentary magnetostatics, "lines of force” are visualized as being fixed 
in space; yet the result (i) just quoted entitles us to identify a field 
line by the material situated on it, so that its motion is taken to be 
that of the material. We then think of the effect of a transverse 
electrostatic field as a translation of the field line itself, with 
velocity E X b/b , in which all charged particles participate. This 
applies strictly only to thermal particles, but holds to some extent for 
energetic particles. 

In the atmosphere, the translational motion of magnetic field lines 
can take place in such a way that B does not change with time. In the 
case in which the field is electrostatic, this must always be the case, 
because then 


<^B/c)t = curl E = curl grad 0 = 0 (VI-14) 

There is then a "circulation" or "convection” in three dimensions (Gold, 

1959), which is possible because the presence of the nonconducting lower 

atmosphere enables the field lines in the upper atmosphere to move 

independently of those in the solid earth. 

The motions of field lines are opposed by frictional collision forces 

acting on the ions in the E region, since v > ox below about 140 km 

and the ions are not "frozen in" to the magnetic field. But since 

v « cd in the dynamo region, the electrons are "frozen in, " and do 
e e 

take part in the motion of the field lines. Thus, the ions and electrons 
in the F region drift with the same velocity as the electrons in the 
E region, apart from small differences due to the vertical changes of 
B. All have the cycloidal types of motion shown in Fig. 30, for the 
case v « cd and move along electric equipotentials . In the E region, 
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the relative motions of ions and electrons represent an electric current 
and the ion-neutral collisions represent joule loss. This is just another 
way of looking at the motions from which the conductivities were derived. 

This description implies that any current flowing in the E layer 
is accompanied by motions of ionization in the F region and magneto- 
sphere. Conversely, if any motion of ionization normal to B is set 
up in the magnetosphere by external causes, it will produce corresponding 
motions of the F region ions and electrons and will set up electron 
currents in the E region, which will dissipate energy. The circulatory 
motion in the magnetosphere, however, is frictionless because of the 
absence of collisions, but has to take place in such a way that the field 
lines are merely interchanged and the field is not distorted in shape. 

The magnetospheric circulation does not directly produce any magnetic 
variations at the ground; these result indirectly from the E region 
currents and induced earth currents, in which the energy is eventually 
dissipated . 

4 . Observational Data 

It is useful to keep in mind the orders of magnitude of the various 
parameters of importance to the dynamo theory. Some of these are listed 
in Table V, and a more detailed discussion of the observational results 
is given below. 

Barometric oscillations have been analyzed in detail, but the com- 
parison between ground-level and dynamo region conditions involves cal- 
culation of the tidal amplification factor and phase variations, using 
an atmospheric model such as done for Fig. 28. As we have discussed 
earlier, this has not been a particularly fruitful approach to wind 
determinations . 

Drifts of meteor trails can be observed by radar technqiues, and 
give winds in the neutral air at 80 - 110 km height, but the data are 
limited to a few locations. A long series of measurements has been made 
at Jodrell Bank (53° N) where the amplitude of the 12-hr oscillation is 
somewhat greater than that of the 24-hr wind (Greenhow and Neufeld, 1961 ); 
whereas at Adelaide (35" s) the 24-hr component is more pronounced 
( Elf ord , 1959). 
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E layer drifts have been measured at various stations. Early 
results were reviewed by Briggs and Spencer (1954) and an analysis of 
I GY data given by Shimazaki (1959). A more recent survey has been given 
by Rawer (1963). Strong semidiurnal components are found, whose vectors 
generally rotate clockwise with advancing time in the northern hemisphere, 
and counterclockwise in the southern hemisphere. Diurnal components com- 
parable to the semidiurnal ones, have been reported but their detection 
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necessitates inclusion of data from nighttime layers, which may not 

be strictly comparable with the normal E layer data. 

F layer drifts are also reviewed in the papers quoted above. When 
the velocities arc subjected to harmonic analysis, the diurnal component 
is found to predominate. Often there arc fairly abrupt changes near 
sunrise and sunset, which perhaps arise from the changes of E region 
conductivity accompanying the formation and disappearance of the E 
layer . 

Magnetic variations of the D, H, Z elements on quiet days follow 
a regular latitudinal pattern and can be represented by overhead current 
systems (see Figs. 25 and 26). The dominant period of variation is 24 
hr, doubtless because of the variation of E region conductivity. These 
daily variations contain an unknown "zero level" which unfortunately 
cannot be separated from the earth’s main field, so that variations of 
each component are conventionally measured from its average level through- 
out 24 hr. However, it would seem much more probable that the real zero 
levels are near the nighttime values, since the nighttime conductivity is 
so small. This is one of the major problems in relating magnetic varia- 
tions to dynamo currents, although it could be resolved by determining 
the nighttime currents by means of a sensitive rocket-borne magnetometer. 

The calculations required in dynamo analysis can be approached from 
two standpoints. One is to solve the equations on a local scale, insert- 
ing observed values of the local magnetic variations and comparing them 
with other data. The most that can be said in this regard is that 
limited agreement is sometimes obtained between the E region data, 
meteor winds, and ground-level pressure oscillations; but that it is 
difficult to relate F layer drifts to the other data. 

To some extent, this may be due to the difficulty of separating the 
electric field into its "induced" and "electrostatic” components. 

Another (and more basic) difficulty is that the "thin sheet" approach 
to dynamo theory fails if the wind velocity W varies considerably with 
height through the conducting regions. There are strong indications 
that it does, and that the height variations may be too irregular to 
enable any meaningful "average" velocity to be defined. 
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The other approach is to assume that a steady state exists and to 
solve the equations on a global scale to obtain the layer current density 
J, electrostatic potential $ and wind velocity W. Not only must these 
quantities satisfy the dynamo equations at all points of the conducting 
layer, but in addition the current continuity equation, div J = 0, must 
be satisfied everywhere. Fortunately, it then becomes possible to sepa- 
rate the "induced" and "electrostatic" components by procedures which 
make use of the identity curl grad 0=0. In addition, some dynamical 
constraint on the wind velocity must be assumed. It is generally 
assumed that the wind velocity W can be derived either from a scalar 
potential (the terms "scaloidal" and "irrotational " being used a describe 
this case) or from a vector potential ("toroidal" or "solenoidal") . The 
latter approach, however, seems to lead to difficulties at the poles of 
rotation. In addition, the Coriolis forces arising from the earth s 
rotation introduce considerable complications, so are often neglected. 

All calculations are based on models for the conductivity g. Since 
E region electron density changes by a factor of 30 (or more) between 
day and night, E undergoes a variation which is predominantly diurnal, 
but contains higher harmonics. Lunar tidal calculations are simpler 
because the moon does not affect the conductivity, provided the data are 
averaged in such a way as to remove the solar effects. 

The calculations of Baker (1953) and Fejer (1953) proceed in the 
"forward" direction, left to right in Fig. 27. They neglect the diurnal 
variation of conductivity and therefore obtain only semidiurnal magnetic 
effects and drifts from their assumed semidiurnal pressure variations. 
Several of the Japanese workers start with the magnetic Sq data and 
proceed in the reverse direction (e.g., H. Maeda, 1955; Hirono and 
Kitamura, 1956; Kato, 1956). The computed E region winds turned out 
to be predominantly diurnal, whereas the observational data we quoted 
earlier (meteor trails, E region drifts, barometric oscillations) sug- 
gest that semidiurnal variations dominate. 

Considering the uncertainties that arise in discussing just the 
dynamo region itself, it is not too surprising that the extension of the 
theory to F region drifts gives few firm conclusions. A number of F2 
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layer phenomena are attributed to electromagnetic drifts, sometimes 
because no other explanation is available. We can see that drifts might 
be important by a simple order -of -magni tude argument, by which we expect 
a vertical drift velocity to cause noticeable effects at the F2 peak 
if it moves ionization through one scale height H in a period equal 
to its average lifetime 1 fa. Thus v ~ Hg = 50 km X 10" 4 sec" 1 = 5 m/sec, 
say. Velocities deduced from dynamo theory are generally larger than this, 
so should be of importance. 

To determine the actual effect of the drift on the electron density 
N(h,t), one must solve the continuity equation with drift and diffusion 
included. This is complicated but has been attempted for some idealized 
conditions. Some conclusions relating to "day equilibrium 11 and "night 
stationary" layers were mentioned in Sec. III-6, but so far little suc- 
cess has been achieved in relating calculated and observed drifts to 
detailed time variations of the F2 layer. 

It may have been noted that the principal unknown in all of the cal- 
culations is the electrostatic field, E . If this could be computed or 
measured with confidence, the applicability of dynamo theory would be 
much better established. Unfortunately, it is extremely difficult to 
measure this quantity in situ , because its magnitude is on the order of 
1 mv/m. Rocket and satellite potentials are typically 100 to 1000 times 
larger and quite variable due to photoemission, induced potential and 
other perturbing factors. 
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VII. STORMS AND THEIR IONOSPHERIC EFFECTS 


1 . Synopsis of Storm Effects 

The events known as magnetic storms comprise a very complicated 
sequence of phenomena. Our account is necessarily condensed, and is 
largely devoted to ionospheric effects. We must caution that most of 
our statements are subject to reservations of one sort or another. Most 
of the theories are tentative, and individual storms vary so much that 
exceptions may be found to almost any description of the observed 
phenomena; 

Magnetic storms are initiated by solar disturbances. Most severe 
storms have a sudden commencement (SC) which takes place 1 to 2 days 
after a solar chromospheric flare, though often there is difficulty in 
associating storms with particular flares and it may be that some SC 
storms originate differently. Another type of storm, generally less 
severe, is not associated with flares. These storms start gradually and 
tend to recur at intervals of 27 days, the period of solar rotation as 
viewed from the earth. They have been attributed to periodic immersion 
of the earth in particle streams emitted, more or less continuously, from 
active areas on the sun which became known as (magnetic) M regions. 

Since M regions could not be positively identified with visual phenomena 
on the solar disk, even though they must persist for some months, they 
remained somewhat hypothetical. It has been suggested (Piddington, 1964) 
that M region storms are associated with certain features in the inter- 
planetary magnetic field. These rotate with the sun, hence causing a 
27-day recurrence, but there need be no association between the storms 
and visible features on the solar disk. M region storms are most common 
during the declining part of the solar cycle. SC-type storms are most 
frequent at sunspot maximum, when large solar flares are commonest. 

Table VI shows some storm phenomena, with notes on their durations 
and possible causes; some of these are reasonably well established but 
others are quite speculative. Any numerical values quoted are intended 
to be mere orders of magnitude, characteristic of strong (but not excep- 
tional) storms. Timewise, the phenomena fall into two categories. Those 
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listed under (l) and ( 2 ) in Table VI, accompany the visual flare, and 
involve electromagnetic radiation travelling at the speed of light. 

Effects (3) to ( 7 ), which comprise the storm proper, depend on the much 
slower propagation of clouds of solar plasma. A typical delay time of 
l-l/2 days between the flare and the SC implies a speed of -10 3 km/s 
from the sun to the earth. But solar cosmic-ray increases and PCA events 
listed under ( 2 ), begin not long after the flare and seem to be due to 
solar protons propagating at velocities of up to 10 5 km/s. 

The storm sudden sudden commencement (SC) is thought be be caused 
by the impact of a shock wave which precedes the plasma cloud on the 
magnetosphere of the earth. The disturbances are then propagated to 
the earth in the form of hydromagnetic waves. There are differences of 
order 1 min in the time of SC at different places, presumably due to the 
travel time of the waves. Sometimes an impulse resembling an SC is 
observed, but no storm ensues, in which case the term "sudden impulse" 

(Si) is used. 

In a typical magnetic storm, the SC is followed by an "initial" or 
"positive" phase lasting a few hours. During this time the geomagnetic 
field intensity is increased, the "stormtime" component Dst ( h ) being 
positive. This is probably due to the compression of the geomagnetic 
field by the solar plasma. After a few hours the main phase sets in, in 
which Dst(H) becomes negative and the field is depressed' for a day or 
so. This reduction can be represented as the field of a "ring current" 
which is oppositely directed to the main geomagnetic field, the location 
and nature of which constitutes one of the central problems of storm 
theory, Dst ( h) reaches its greatest negative values about 18 - 36 hr 
after the SC (soonest in great storms). There are often irregular fluc- 
tuations of fields during the initial and main phases, but after I- 1/2 - 
2 days the storm passes to the recovery phase and the field returns 
smoothly to normal, with an exponential time constant of a day or so. 

The amplitude of the Dst magnetic field is greatest in low latitudes, 
and decreases polewards. This is opposite to the field of the current 
systems associated with the DS disturbance local-time inequality. This 
field is most intense at auroral latitudes, where it exceeds Dst, but 
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decreases towards low latitudes, where Dst is the dominant storm 
variation. In general, the variation observed at any station is a com- 
plicated combination of DS and Dst, which depends on the local time at 

which the SC occurs. 

In early stages of a storm, the DS system changes rapidly, but after 
some hours it settles into a pattern which remains nearly fixed with 
respect to the sun, and gradually decays as the storm progresses. Once 
these patterns are established, the harmonic components of the DS varia- 
tion at any station maintain almost constant phases in local time, though 
the amplitudes decay gradually. A comprehensive study of weak, moderate 
and great magnetic storms has been made by Sugiura and Chapman (1958), 

Much of the interest in geomagnetic storms is centered on the aurora, 
cosmic rays and the radiation belts. These are included under (4) and 
(5) in Table VI, but we do not propose to discuss them further here. The 
mid-latitude red arcs are a fairly recently discovered phenomenon, which 
promise to be of considerable interest in relation to upper atmosphere 
processes. An account of many of their known properties has been given 
by Roach and Roach (1963). 

The increase of satellite drag during storms is almost certainly due 
to an increase of atmospheric temperature. A detailed analysis of this 
effect has been given by Jacchia (1961, 1963) who find that the tempera- 
ture in the upper thermosphere increases by 1.0 to 1.5 °K per unit of A^ 
It has been suggested that joule heating by DS currents is responsible 
for these temperature increases in the neutral atmosphere (Cole, 1962). 
The connection between current flows and the electric fields which might 
cause mid-latitude "red arcs" has been discussed by Carleton and Megill 
(1964). So far, it has not been possible to assess the relative impor- 
tance of this mechanism, corpuscular heating, and other possible causes 
of heating such as absorption of hydromagnetic waves (Dessler, 1959). 
Waves of periods around 1 sec are strongly absorbed in the ionosphere 
( Karplus , Francis and Dragt, 1962); there seems to be no evidence that 
the intensity of such waves incident on the ionosphere from above is 
sufficient to account for the heating. 
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2 . Storm Effects in the Lower Ionosphere 


Disturbances of the lower ionosphere during storms are most pro- 
nounced in high latitudes. This scarcely surprising, because of the 
known precipitation of energetic particles in the auroral zone. Some 
effects, however, occur at other latitudes. 

It is well known that radio communications are affected during 
magnetic disturbance, the terms "fadeout" and "blackout" being in common 
use. Though part of the difficulties in ionospheric propagation arise 
from decreased F2 critical frequencies, to be discussed in the next 
section, they are largely due to D region absorption. Reid and Collins 
(1958) showed that high latitude absorption comprises three distinct 
types of phenomena which we list here. Their "type numbers" are not 
now in general use but are sometimes found in the literature. 

IyP e 1 ls the "short wave fadeout" (SWF), part of the "sudden iono- 
spheric disturbance" (SID) which accompanies solar flares (Dellinger, 
1937). This is observed all over the sunlit hemisphere, so is not par- 
ticularly a "high-latitude" phenomenon. 

T -yP e 11 is the "auroral absorption" which accompanies the main phase 
of geomagnetic storms. It is associated with visual auroras or those 
detected by radar, and its positional distribution has been found to 
follow well-defined "spiral" patterns, as we shall see later. 

TyP e 111 is the "polar cap absorption" (PCA) which begins a few hours 
after solar flares--and thus, well before the storm SC--and persists for 
some days, with pronounced diurnal variations. Maps showing the develop- 
ment of PCA have been published by Obayashi (1959, I960). 

The "sudden ionospheric disturbance," SID, which accompanies solar 
flares, is associated with a number of other phenomena which begin 
simultaneously, including: 

SWF. Short wave fadeout or "Type I" absorption., Waves reflected from 
the E and F layers are strongly attentuated by enhanced ioni- 
zation in the D region, 

SCNA: Sudden cosmic-noise absorption. A fairly large flare may cause 
a 2 db increase in absorption or even more at a frequency of 
25 Mc/s ( Bhonsle, I960), 
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SPA- Sudden phase anomaly. The change of phase of long waves reflected 
at oblique incidence indicates a drop of several kilometers in the 
height of reflection (Bracewell et al , 1951; Bracewell and Straker, 

1949) . 

SEA: Sudden enhancement of atmospherics recorded at frequencies around 

20 kc/s (Ellison, 1953; Bureau, 1937). Like the SPA, this indi- 
cates a change of VLF propagation conditions, 

SFE: or crochet: Magnetic solar flare effect. A disturbance of the 

geomagnetic field, attributed to increased ionospheric conduc- 
tivity, which leads to increased current flow. 

All these effects indicate an increase of electron density in the 
D region, perhaps by a factor of ten in a major flare. The fact that 
they begin at the time of visual sighting of the flare implies that they 
are due to electromagnetic, rather than corpuscular, radiation. The 
fact that ionization is produced at lower heights than usual provides 
strong evidence that the radiation is harder (i.e., possesses a smaller 
absorption cross-section) than that responsible for the normal D layer. 
It is known that solar X ray emission is greatly enhanced during a 
flare, whereas Ly Ct emission is only slightly increased. 

With regard to the magnetic crochet, we recall that the conductivity 
per ion pair is so small below 80 km that huge changes of electron density 
would be required to enable sizable electric currents to flow. Thus, it 
may be that the currents producing the crochet flow in the E region, 
and this might be connected with the increase of f Q E which has been 
reported on occasions (Taubenheim, 1957). But information on electron 
density changes above the D region, during solar flares, is rather 
incomplete, largely because of absorption in the D region. 

Solar flare phenomena were observed in an extreme form during the 
famous event of February 23, 1956. A few minutes after the start of the 
flare, which was accompanied by an intense SID on the sunlit side of the 
earth, a worldwide increase of cosmic rays was observed and strong D 
region disturbances occurred in the dark hemisphere. Many papers on this 
event have appeared, and a summary of the observations has been given by 

Bailey (1959). 

During the main phase of storms, the lower ionosphere is not drasti- 
cally affected at middle latitudes. There is some evidence that the 
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critical frequency f E is depressed by a few percent (Sato, 1957), 

D region absorption is almost unaffected, but there are abnormal diurnal 
variations of amplitude and phase of low frequency waves (<300 kc/s) 
reflected from the ionosphere (Bracewell et al , 1953). Belrose has 
found that these persist for two weeks or more after the storm (Ratcliffe 
and Weekes, I960), and similar long-lived phenomena have been reported 
after a nuclear explosion (Brady et al , 1964). Thus, the D region may 
possess a long time constant for recovery from storm effects. 

An early description of high-latitude storm effects was based on 
observations at Tromso during the Second Polar Year of 1932/33 (Appleton 
et 1937). In recent years, much study has been devoted to the dis- 

tribution of "polar blackouts" and storm sporadic E, using routine 

tabulations of ionospheric parameters such as f Es and f (the 

o min ' 

minimum frequency at which ionospheric echoes are received, this being 
determined principally by absorption). The occurrence of these phenomena 
shows marked maxima at certain times of day, which vary from station to 
station. On polar maps plotted in magnetic coordinates, "isochrons" can 
be drawn to join points at which the maxima occur simultaneously (in 
Universal Time), and these are found to be spiral in form (Agy, 1954; 
Thomas, 1962). A similar spiral pattern has been found for magnetic 
activity at high latitudes (Nikol T skii, 196l), and for auroral activity. 

Originally, these were interpreted as Stormer spirals, which are the 
locus of precipi tat ion of charged particles originating in a distant 
source and deflected by the geomagnetic field (stormer, 1917). If this 
were correct, then the sense of the spiral would indicate that blackouts 
and storm Es are produced by positively and negatively charged par- 
ticles, respectively. However, it is now clear that Stormer T s theory 
is not applicable to particles in the required range of energy, because 
their density must be high enough for interactions between particles to 
be important. The Stormer theory neglects such interactions and is prob- 
ably valid only for high-energy cosmic ray particles whose motion is 
influenced only by the geomagnetic field. Instead, the explanation of 
the observed spirals has become a part of modern theories of geomagnetic 
storms . 
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3. Storm Effects in the F Region 


Solar flare effects in the F layer have proved difficult to detect, 
and were considered to be absent by Berkner and Wells (1937), This prob - 
ably implies that the effects produced by a flare lasting (say) half an 
hour are no larger than the random variations of N m F2 in such a period, 
which might be 10 percent. Knecht and Davies (l96l) have reviewed a few 
cases in which increases in N^F2 have been observed. These increases 
range from 20 to 50 percent and generally disappear within about one 
hour. The Orbiting Solar Observatory recorded a 14 percent increase in 
the intensity of solar Hell 304 % radiation during a fairly large (Class 
3) flare (Lindsay, 1963), so the F2 layer effects may well be due to 
increased photoionization. 

In contrast to these rather elusive flare effects, changes in the 
F layer during a storm are profound. The critical frequency f Q Fl is 
not greatly affected; Sato (1957) found a slight decrease of about 0.3 
Mc/s (which implies a reduction of about 10 percent in N^l). We may 
take this to indicate that the production rate q is not greatly changed 
during a storm, as we might indeed expect. But the reduction of F2 
layer electron density sometimes results in the FI cusp becoming 
visible on ionograms at times when it would not normally be seen. This 
is probably due to the reduction of electron density in the overlying 
F2 layer, which would tend to make the FI ledge more conspicuous, as 
we discussed in Sec. III-4. 

During storms, the h'(f) curves for the F2 layer are greatly 
changed. The critical frequencies f Q F2 are usually reduced and the 
virtual heights h' greatly increased, sometimes attaining values of 
600 km or more. These were originally thought to indicate real changes 
in the height of the F2 layer, but this interpretation was shown to 
be false when the N(h) profiles were computed (Thomas and Robbins, 
1958). The real height of the F2 layer is actually the same, to within 
about ±20 km, on storm days and quiet days. The great virtual heights 
occur because the disturbed electron density distributions are such that 
the virtual height integral, h' = / p'(f,N) dh, takes very large values. 
We say that the waves reflected from the F2 layer suffer very large 


SEL-64-111 


168 



"group retardation" below the height of reflection. If follows that 

increase of h r F2 should be regarded only as a sensitive indicator of 

F layer disturbance. Many authors have used semi empirical measures of 

F2 layer height (deduced from ionograms by assuming parabolic N(h) 

profiles, or from published maximum-usable frequencies), but it is not 

clear that these parameters give reliable information about h F2 . 

m 

The storm changes in f Q F2 can be analyzed in terms of the SD and 
Dst variations originally applied to magnetic parameters. This was 
done by Martyn (l953x; 1953y), and a more extensive analysis was carried 
out by Matsushita (1959), who used data for 109 storms at 38 observa- 
tories situated in eight latitude zones. Matsushita’s Dst curves are 
reproduced in Fig. 32. 

In the Zones 1-3 (geomagnetic latitudes 45 deg to 60 deg) the 

Dst curves are positive for the first few hours and then become negative, 

reaching a greatest depression of 30 percent in N F2 (for strong 

m 

storms) about 24 hr after the SC. In geomagnetic latitudes 20 deg to 

45- deg, the "positive" phase lasts longer and the "negative" phase is 

weaker (-10 percent depression of N F2). In low latitudes, there is 

m 

no negative phase and Dst is weakly positive (+5 percent to +10 percent 
in N^F2) throughout the storm. 

It seems therefore that ionospheric storms can be either positive or 

negative, according to whether N F2 is increased or decreased from 

m 

mean quiet day values. Positive storms appear to be more common in winter 
(Sato, 1957). Negative storms generally possess an initial positive 
phase, so that the curves of Dst (N^F2) possess a superficial resem- 
blance to magnetic Dst(H) curves. However, we may note that the posi- 
tive phase lasts longer in the N F2 variation than in the magnetic 

m 

field variation. Some detailed studies have been presented by Appleton 
and Piggott (1952), by J. W. King (l96l) and by Benkova (l96l). Gen- 
erally, the greatest depressions of N F2 are observed on the day after 

m 

the greatest disturbance of the magnetic field. 

In addition to the Dst components, there are SD (local-time) com- 
ponents of the variation, which display fairly complicated changes of 
amplitude and phase, especially during the earlier parts of the storm, 
but gradually die away in the recovery phase. 
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STRONG STORMS WEAK STORMS 



FIG. 32. Dst VARIATIONS OF N m F2 IN EACH 
OF EIGHT LATITUDE ZONES FOR STRONG AND 
WEAK MAGNETIC STORMS WHICH EXHIBIT SUDDEN 
COMMENCEMENTS , The ordinate is the 
approximate percentage deviation from the 
quiet day behavior vs storm time. The 
zone number is shown in parenthesis 
between the applicable geomagnetic lati- 
tudes: 60° ( 1 )55° ( 2)50° ( 3)45° (4)40° (5 ) 

30° ( 6)20° ( 7 )l0° (8 ) -10° (Matsushita, 
1959). 
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Satellite and moon-echo observations show that the total electron 

content / Ndh is generally reduced during storms, though probably to 

a lesser extent than the peak electron density (Lawrence et al , 1963), 

This would imply that the F2 layer is slightly thicker during storms, 

although Titheridge (1964) finds little change. The reduction of total 

electron content was extremely pronounced during a severe storm in 

November 1960 (Taylor, 1961 ). In addition, the electron content increased 

much more slowly after sunrise on the disturbed day than on quiet days. 

It had previously been known that dN(h)/dt and dN /dt were reduced 

m 

during storms, and the total content observations show that the reduction 
cannot be interpreted in terms of a mere redistribution of ionization. 
Therefore, the observations suggest that the storm interfered with the 
F2 layer photochemical processes in some manner. However, the prime 
cause of these effects must be explained by more comprehensive theories 
of geomagnetic storms (Sec. VII-4). 

These F region storm effects present many unsolved problems. We 
can do no more than review briefly some possible explanations, which 
include "thermal expansion," "photochemical changes" and "electromagnetic 
drift" theories. Very possibly, all of these contribute to F layer 
storm phenomena. 

Early papers (Berkner, Wells and Seaton, 1939) attributed the depres- 
sion of electron density to thermal expansion of the F region, this 
being largely suggested by the increase of F2 layer virtual height, 
which is now known to be misleading. Although satellite drag studies 
have shown that the atmosphere is indeed heated during storms, this 
explanation is probably inadequate, and could not account for the occur- 
rence of positive storms in which N F2 increases. A more recent dis- 

m 

cussion has been given by Matuura (1963). 

Thermal expansion of the F region would be expected to cause some 

increase in h F2 and a modest decrease in N F2. If it were the main 
m m 

process, there should then be a correlation between the storm changes 
N^F2 and h ^F2 w k ic h has not, in fact, been shown to exist. In 
order to explain large decreases of N , it has been suggested that the 
loss coefficient is enhanced, either because of a postulated temperature 
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dependence of the ion-atom interchange coefficient (Yonezawa, 1963) or 
by changes of atmospheric composition which increase the proportion of 
molecular gases in the F region (Seaton, 1956; Kamiyama, I960). Such 
a change might conceivably occur if dynamic processes in the vicinity 
of 100 km were enhanced such as to increase turbulence and raise the 
upper limit of mixing of molecular gases (0 2 and N g ) with atomic oxygen, 
as discussed by King (1962). 

The drift theory of Martyn (l953y) attributes storm effects to 
electromagnetic movements, produced by the electric fields associated 
with the storm current systems. Although this seems quite a reasonable 
concept, the detailed calculations have not led to a really satisfactory 
model. This is because they resemble a complicated chain which contains 
weak links in the form of unproven assumptions. It has very much in 
common with the "dynamo theory scheme" of Fig. 27. The drift calculations 
can be divided into three main stages: 

(a) Computation of ionospheric storm current systems from the mag- 
netic data. This requires assumptions about the contribution of 
induced ground currents to the magnetic effects. 

(b) Calculation of electric fields driving these currents, by 
inserting conductivities, and then the derivation of F region 
drift velocities by assuming transmission of the electric field 
up the geomagnetic field lines. On the simplest theory, these 
velocities are = E g X b/b 2 . 

(c) Insertion of the drift velocity Vj) into the continuity equation 
and calculation of the resulting change of electron density. In 
general, the vertical component of Vj) would be expected to 
produce the greatest effects, but these are likely to be modified 
by the effects of ion-drag (Baker and Martyn, 1953). 

We do not further discuss (a) and (b) here, because the problems arising 
are very similar to those discussed earlier. One simplification exists, 
however, because the storm electric fields are unlikely to be produced 
by dynamo action, but are probably entirely electrostatic in nature, so 
that no separation of the field is necessary. Regarding (c), there have 
been considerable doubts about even the sense of changes of N produced 
by given drifts, because of the complexity of the F layer continuity 
equation. Maeda and Sato (1959) attributed the occurrence of positive 
and negative storms at different latitudes to a latitude-dependent com- 
bination of quiet day and storm (Sq and DS) drift velocities. 
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Complete solutions of the F2 layer continuity equation with time- 
varying drift, diffusion and photochemical processes, will doubtless be 
developed. One can obtain some idea of the effects of drifts quite 
simply. An upward drift transports ionization into a region of lower 
loss, thus prolonging the average lifetime of ionization. Thus, during 
the day, upward movements would certainly increase the total electron 
content / Ndh, and would probably increase N . A downward movement 

TTl 

would have opposite effects. At night, once a "stationary" layer were 
established, upward or downward movement would alter the effective rate 
of decay of the layer. 

Some authors attribute ionospheric DS effects to DS currents, and 
ionospheric Dst effects to Dst ring-currents, the latter being sug- 
gested by the superficial similarity of typical Dst ( h) and Dst (n F2) 

m 

curves. In our opinion, the similarity may be accidental. Most of the 
Ost ring current is thought to flow at geocentric distances of several 
earth radii, and the presence of Dst current in the ionosphere has not 
been established. Because conductivities and production rates are dif- 
ferent by day and by night, a periodically-varying drift velocity with 
average value of zero can produce an average non-zero change of electron 
density. In this way a periodic (DS) current system in the ionosphere 

can give rise to both Dst and DS effects in N F2. This is consistent 
^ m 

with the argument that there are differences in the magnetic SD varia- 
tions for ionospherically positive or negative storms (Rishbeth, 1963). 

4 . Theories of Geomagnetic Storms 

In the last section, we saw that no detailed explanation yet exists 
of how the storm effects in the ionosphere arise. But the problem of 
accounting for the geomagnetic phenomena is perhaps more basic. We can 
only attempt to summarize some of the theories which have been advanced. 

The most obvious features of the magnetic disturbance are the storm- 
time Dst decrease of the field, which can be represented by a ring 
current, and the auroral zone electrojets which belong to the DS current 
system. There is a case for saying that the electrojets are associated 
with the prime cause of the DS currents, and that DS currents outside 
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the auroral zone largely represent return flows of current which complete 
the circuit; though Akasofu and Chapman (1964) have presented evidence 
for a different view. 

If the DS auroral electrojets are to be produced by wind-driven 
dynamo action within the atmosphere, there must exist very localized 
wind systems (with velocities ~1 km/s) or strong enhancements of iono- 
spheric conductivity (Nagata and Fukushima, 1952; Cole, I960). The 
latter seems more plausible, because the electron density in the auroral 
zone is certainly enhanced by ionization produced by incoming charged 
particles, especially during storms. In the absence of any other driving 
mechanism, the currents would have to be sustained by abstracting more 
power from the atmospheric tides than is absorbed under quiet conditions. 
Since the height-integrated current density in the jets is of order 
3 amp/m for an intense storm, as compared to a quiet -day current density 
of 0.1 amp/m in the mid-latitude ionosphere, the increase in conductivity 
would have to be roughly 30:1 if no other driving mechanism existed. It 
is not clear whether the DS currents could arise purely by this means, 
without any active driving mechanism being present, so it seems reason- 
able to look for driving mechanisms within the auroral zones, or in the 
parts of the magnetosphere linked by field lines to the auroral zones. 

The review by Nagata (1963) discusses these possibilities. 

The pioneer development of geomagnetic storm theory was carried out 
by Chapman and Ferraro (1931 and later papers). They solved idealized 
mathematical problems in relation to the interaction of neutral ionized 
streams projected from the sun with the geomagnetic dipole field. They 
found that oppositely charged particles tend to be deviated towards 
opposite sides of the earth; some particles might enter captured orbits 
and would give rise to a westward ring current far outside the earth. 

This might produce the Dst decrease of magnetic field, but it is not 
so clear how the DS current systems arise. 

Though historically important, the Chapman-Ferraro theory is inade- 
quate because it treats interplanetary space as a vacuum, ignoring both 
the ionized interplanetary matter, which may be regarded as an extension 
of the solar corona, and the conducting plasma in the earth's magneto- 
sphere. Disturbances can be transmitted rapidly through the magnetosphere 
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in the form of hydromagnetic waves, and most modern theories use this 
idea in one form or another. There is little doubt that the magnetic 
SC is caused by the arrival of a shock wave preceding the solar plasma 
cloud, and is transmitted to the earth’s surface as a hydromagnetic 
wave. Then the initial phase Dst increase of the geomagnetic field 
is ascribed to a sustained pressure exerted by the impinging solar 
stream. 

In some theories, it is supposed that ■ solar plasma streams interact 
with the magnetosphere as they sweep past it, and that the main phase 
Dst decrease and the DS current systems arise from the interaction. 

Thus, Piddington (i960) attributes the jnain phase Dst decrease to an 
outward force exerted on the geomagnetic field by the solar stream, which 
causes field lines to be drawn out into a geomagnetic "tail," He regards 
the DS fields as arising initially from transmission to the earth of 
hydromagnetic "twist" waves generated in the vicinity of the magneto- 
pause. 

Axford and Hines (l96l) envisage a gigantic circulation in the mag- 
netosphere, driven by frictional interaction with the solar wind at the 
magnetopause. Because the plasma is frozen to the magnetic field lines, 
the circulation takes the form of a "convection” of field lines, as 
described in Sec. VI -5. Figure 33 shows the proposed form of the elec- 
tric equipotentials , which are also streamlines of the convective motion, 
in the earth ’ ^equatorial plane . The pattern at ionospheric heights is 
then obtained by mapping the streamlines along geomagnetic field lines. 

The solar wind interaction generates turbulence in the outer mag- 
netosphere, and the turbulent material is conveyed by the circulation 
into the interior of the magnetosphere. The regions in which turbulence 
is expected to be present, when mapped along magnetic field lines, cor- 
respond to a band which resembles the locus of peak magnetic and auroral 
activity, and the high-latitude ionospheric phenomena described in Sec. 

VII -2 . 

We have devoted some space to a description of the Axford and Hines 
theory, as its concepts overlap to some extent those of other theories,.. 
Nevertheless, there are essential differences between different theories 
which we have not attempted to bring out. 
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SOLAR WIND 



FIG. 33. STREAMLINES OF THE MAGNETO SPHERIC CONVECTION 
IN THE EARTH’S EQUATORIAL PLANE. The curves also 
represent electric equipotentials, and are obtained 
by combining a symmetrical pattern caused by inter- 
action between the solar wind and the outer magneto- 
sphere, and an asymmetric pattern caused by rotation 
of the magnetosphere. An additional circulation may 
exist in the inner magnetosphere (after Axford and 
Hines, 1961 ). 


In reality, the magnetosphere contains not only the thermal plasma 
but also a wide spectrum of energetic particles. Low energy particles 
take part in the circulation, but high energy particles trapped in the 
geomagnetic field describe "mirror" orbits, and produce a "ring current" 
because they precess around the earth. Of course, no sharp distinction 
exists between "low" and "high" energies, though for rough calculations 
1 kev might be taken as the dividing line. In principle, the precession 
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of energetic ions and electrons can give rise to a ring current which 
could produce the Dst field in the main phase of a storm. The mechanisms 
by which particles enter and escape from the belts form a very active field 
of experimental and theoretical research. 

Kern (1962) has suggested that, should asymmetries exist in the 
distribution of charged particles, a separation of charge would occur 
which might be capable of producing the DS current system. The mechanism 
proposed by Fejer (1963) does not require any asymmetry, but requires 
the trapped radiation to contain a slight preponderance of particles of 
one sign (e.g., protons). Electrical neutrality of the magnetosphere 
as a whole is then preserved by a slight excess of oppositely charged 
particles (e.g., electrons) in the low energy plasma which takes part 
in the large-scale circulation of the magnetosphere. Because the paths 
of high-energy and thermal particles, when projected into the equatorial 
plane, are different, electric fields are set up, and are neutralized by 
flow of current along geomagnetic field lines and through the polar 
ionosphere. If the sign of the excess charges is as indicated above, 
the theory can qualitatively account for the flow of the polar electro- 
jet current. More recently, Fejer (1964) has extended his analysis to 
the calculation of worldwide DS current systems resulting from this 
mechanism. Promising results have been obtained with a relatively crude 
model, in which diurnal variations of ionospheric conductivity are 
neglected. 

We may safely say that, in spite of the great progress achieved in 
recent years, the problems posed by geomagnetic storms will remain of 
interest for many years to come. Indeed, the same can be said for upper 
atmosphere physics as a whole, as we have tried to show in this account. 
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